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Experiments on the Production of Very Low Temperatures by 
Two-Stage Demagnetization 


Braye OARDY, J2 HATTON, Bb. V. ROLLIN, E. F. W. SEYMOUR 
AND H. B. SILSBEE 
Clarendon Laboratory, Oxford 


MS. received 5th Fune 1951 


ABSTRACT. A two-stage demagnetization apparatus employing a superconducting wire 
as the thermal switch between the stages is described. "Temperatures in the region of 
0-003° k. have been obtained using a magnetic field of only 4,200 gauss. The problem of 
thermal insulation was investigated and the residual heat inflow after demagnetization was 
reduced to one erg/min. The possibility of employing a superconducting switch in cyclical 
magnetic cooling systems and in nuclear demagnetization is discussed. 


§1. INTRODUCTION 

HE success of the method of cooling by adiabatic demagnetization depends 
| to a large extent on the amount of entropy which is removed in the initial 
isothermal magnetization. In order to remove an appreciable fraction 
of the electron spin entropy in a magnetic field H at temperature T we require 
a value of H/T greater than about 15 kilogauss per degree. The magnetization 
is normally carried out at about 1°k., which is the lowest temperature obtainable 
by pumping liquid helium, and large water-cooled magnets are necessary to 
produce the required field. ‘The use of such large magnetic fields could be 
avoided by employing a two-stage process in which the magnetization of the 
second stage is carried out at a temperature lower than 1°. obtained by 

demagnetization of the first stage. 

The practical realization of a two-stage process has been impeded by the 
technical difficulty of obtaining a satisfactory way of making and breaking 
thermal contact between the stages. In experiments by Mendoza (1948) contact 
between two magnetic stages was made by copper foil which was mechanically 
broken after magnetization of the second stage. ‘These experiments demonstrated 
the possibility of the method, but the subsequent heating of the second stage was - 
exceptionally rapid. 

Experiments on the thermal conductivity of superconductors at low 
temperatures (de Haas and Rademakers 1940, Rademakers 1949, Hulm 1949, 
Heer and Daunt 1949, Mendelssohn and Olsen 1950) have led to the development 
of an alternative method for obtaining a satisfactory thermal switch. This 
depends on the fact that for some pure superconductors there is a considerable 
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difference between the thermal conductivity in the Ue Ona CIE and 
non-superconducting states. 

We have made a number of experiments with a two-stage apparatus, based 
on the application of this property. The apparatus was designed to enable 
temperatures below a hundredth of a degree to be obtained, using a low-power 
magnet giving a field of only about four kilogauss. 


§2. THE THERMAL SWITCH 


In a normal metal, at temperatures sufficiently low for the mean free path 
of the conduction electrons to be independent of temperature, one would expect 
the thermal conductivity to be proportional to the temperature. For 4 
superconducting metal the electrons in the superconducting state make no 
contribution to the thermal conductivity, so that at temperatures appreciably 
below the transition temperature the conductivity should be of the same order 
as that of a non-metal and should vary as the cube of the temperature. According 
to Heisenberg’s theory of superconductivity (Heisenberg 1949), the ratio of the 
conductivities in the superconducting and normal states is given approximately 
by K,/K,=2(T/T.)?/[1+(T7/T,)*], where 7, is the superconducting transition 
temperature. 

In experiments by de Haas and Rademakers (1940) it was shown that the 
thermal conductivity of lead was much less when in the superconducting state 
than when the superconductivity was destroyed by application of a magnetic 
field in excess of the threshold value. At temperatures below 4°x. the 
conductivity in zero field was very small and, as the temperature was lowered, 
decreased approximately as T?. Measurements by Heer and Daunt (1949) on 
tin and tantalum have confirmed that the ratio of the conductivities in the 
superconducting and normal states is approximately proportional to (7/T,)? 
for T<T,.. It has, however, been found by Mendelssohn and Olsen (1950) 
that quite different results are obtained for columbium and also for 
superconducting alloys, so that it appears that the relation is valid only for certain 
pure superconductors. 

The application of a magnetic field to a superconducting wire as a method 
for varying the thermal contact between specimens in demagnetization experiments 
was suggested by Gorter (1948), and experiments demonstrating the use in a 
single stage apparatus have been made by Heer and Daunt (1949). In the 
experiments of Heer and Daunt a paramagnetic salt was connected to a helium 
cryostat at 1:1°K. by a tantalum wire. When the salt was at a temperature in 
the region of 0-6°K., the rate of heat inflow when the wire was superconducting 
was about twenty times smaller than when the superconductivity was destroyed 
by a magnetic field. 

For a satisfactory thermal switch it is clearly desirable to have the highest 
possible ratio between the normal and superconducting thermal conductivities. 
If 7, and T, are the temperatures of the ends of the wire, and if the conductivities 
K, and K,, in the two states vary as T° and T respectively, then if 7,> 7, the 
rate of heat flow in the two cases is approximately proportional to T,4 and 7,?, 
so that the ratio varies as 7,?.. As a rough approximation this ratio can be 
equated to (7,/7,)”, so that in order to obtain a ratio of 1: 1,000 the temperature 
of the warmer end of the wire should not exceed about one-thirtieth of the 
transition temperature, 


Two-Stage Demagnetization Experiments 863 


§3. CHOICE OF PARAMAGNETIC SALTS 

The considerations affecting the choice of the paramagnetic salts for a 
two-stage apparatus are rather different for the two stages. For the first stage 
we require: (i) the entropy removed on magnetization to be large and, in any 
case, to be considerably greater than the entropy of magnetization of the second 
stage; (ii) the temperature after demagnetization to be low enough to enable 
an appreciable degree of magnetic saturation of the second stage to be obtained 
with the available field. For our purpose it appeared that the most suitable 
substance to use was iron ammonium alum. 

For the second stage we require a salt which will attain a very low temperature 
when demagnetized from conditions approaching magnetic saturation. In the 
ideal case, the temperature should be determined only by the magnetic interaction 
of the ions, which could, if necessary, be reduced by dilution. In practice, 
other effects, such as exchange interaction, Stark effect and interactions between 
electrons and nuclei are often the limiting factors. In the preliminary experiments 
(Darby, Hatton and Rollin 1950) we used diluted copper potassium sulphate, 
but there are considerable nuclear and exchange interactions in this salt (Benzie 
and Cooke 1949), and in subsequent experiments we decided it would be better 
to use dilute mixed crystals of potassium aluminium alum and potassium chrome 
alum, containing 5° of potassium chrome alum. Although this substance is 
not free from all defects it is reasonably satisfactory and, furthermore, its thermal 


and magnetic properties at very low temperatures have already been investigated 
(de Klerk, Steenland and Gorter 1950). 


S44 DESCRIPTION OF THE APPARATUS 
The two stages were cylinders of compressed salt, the first stage of diameter 
16mm. and length 5-7cm., and the second stage of diameter 9mm. and length 
2:3cm. ‘The cylinders were suspended in a framework by an axial cotton thread 


0 5 10 15 20 
Time (minutes) 


Figure 1, Figure 2. 


which was maintained in tension by springs (Figure 1). Each cylinder was 

provided with six copper wires of diameter 0-2mm. around which the salt was 

tightly compressed. ‘These wires enabled thermal contact to be made between 
56-2 
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the salt and the lead wire which served as the thermal switch. The lead wire 
was made from lead of the highest purity obtainable and was of length about 3 cm. 
and diameter 0-3 mm.; it was situated midway between the cylinders. 
The framework containing the two stages was inserted into the helium 
cryostat through a ground joint at the lower end which was sealed with glycerine. 
The magnetic susceptibility, and hence the Curie temperature, of each 
stage was obtained by measuring the mutual inductance between coils mounted 
outside the cryostat. The separation between the stages was such that the upper 
salt had no disturbing influence on the measurement of the susceptibility of the 
lower salt. The magnetic dilution of the second stage was great enough to make 
correction of the 7* values, due to the non-spherical shape of the specimen, 
negligible. The measurement of mutual inductance was made by means of an 
audio-frequency bridge of the Hartshorn type operated at a frequency of 25 c/s. 
The magnet used had pole faces 6 in. high and 4 in. wide with a pole gap 
of 2 in., and gave a field of 4,200 gauss with a power consumption of 600 watts. 
The height of the pole faces was sufficient to enable both stages to be magnetized 
simultaneously, which is a desirable, though not essential, feature. 


§5. EXPERIMENTAL PROCEDURE 

(i) The temperature of the cryostat was reduced to about 1°xK. and both 
stages were magnetized in a field of 4,200 gauss, the heat of magnetization being 
removed by helium exchange gas at a pressure of about 5 x 10-4mm. Hg. 

(ii) After about 10 minutes the exchange gas was pumped away and the 
pumping continued for a further 30 minutes in order to remove as much as 
possible of the gas adsorbed on the surface of the cryostat. 

(11) ‘The magnet was slowly lowered so that the upper stage was demagnetized 
while the lower stage remained in the field. During this time the lead wire was 
in a field greater than the threshold value and, as the cooling proceeded, heat was 
transferred from the lower stage to the upper. This operation was carried out 
very slowly so that the temperature difference between the two stages was as 
small as possible, and thus the minimum amount of entropy was transferred to 
the upper stage. 

(iv) When the lowering of the magnet was completed, the lower stage was 
still in a field of 4,200 gauss and the lead wire was in a stray field of about 
1,000 gauss, and so still in the non-superconducting state. The stray field of 
about 500 gauss on the upper stage was now compensated by means of a small 
permanent magnet. 

(v) After waiting about 30 minutes for equilibrium to be established, the 
magnetic field was reduced to zero. As soon as the stray field on the lead wire 
fell below 800 gauss, the wire became superconducting and so the thermal 
conduction between the stages was reduced. 

(vi) The magnet and compensator were removed and the Curie temperature 
of the two stages was measured. ‘The Curie temperature was obtained from the 
readings of the mutual inductance bridge using a calibration previously obtained 
in the temperature region between 4:2°K. and 1-5°K. 


So. RESULTS 
The lowest Curie temperature reached after demagnetization was about 
0-009°. From the measurements of de Klerk, Steenland and Gorter (1950) 
this corresponds to an absolute temperature of about 0:003°K. In order to 
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obtain this temperature in a single stage demagnetization from 1-:05°K. a magnetic 
field of 15,200 gauss would be required. The temperature of the upper stage 
was about 0:25°K. 

In the initial experiments the heating rate of the second stage after 
demagnetization was rather rapid, being sometimes as high as 50erg/min. 
However, after a number of modifications had been made to the apparatus, a 
heating rate of only about 1 erg/min. was obtained. In Figure 2 is shown a 
typical heating curve for the second stage. The temperature remained below a 
hundredth of a degree for about 3 minutes. The heating rate of the first stage 
was negligible, being only of the order of 1:5x10-* degree per minute, 
corresponding to a heat inflow of 50 erg/min. 


§7. EXPERIMENTS ON HEAT INFLOW 

In experiments at very low temperatures one of the most important 
considerations is the rate of heat inflow to the specimen. In previous experiments 
by Cooke and Hull (1942) it has been shown that with suitable experimental 
arrangements the heat inflow can be reduced to the order of one erg per minute, 
though other workers, in this and other laboratories, have usually observed much 
larger heat inflows. 

The high heat inflow which was observed in our initial experiments was first 
attributed to conduction along the superconducting wire. Experiments with 
different lengths of wire, however, disproved this hypothesis. In order to find 
whether the heat inflow was due to the presence of residual gas in the cryostat 
and to conduction along the supporting threads, the apparatus was modified so 
that the lower stage was surrounded by a copper cylinder which was soldered to 
wires embedded in the upper stage. In this arrangement the lower stage was 
completely shielded from any heat inflow from the cryostat. Although this 
arrangement made some improvement, the heat inflow was still large. After a 
number of further experiments, the origin of the anomalous heating was finally 
traced to the effects of small vibrations which produced heating due to frictional 
losses in the suspended system. In order to eliminate this effect it was found 
necessary to insulate the apparatus from all vibrations due to the vacuum pumps. 
It was found that even the vibration due to the boiling of the mercury of a small 
diffusion pump caused a heating of the upper stage of more than 300 erg/min. 
A similar effect has been observed in experiments by Hull, Wilkinson and Wilks 
(1951). In the experiments of Cooke and Hull (1942) we understand that the 
cryostat was, by chance, well insulated from vibrations, so that this may account 
for the discrepancy between their results and those of subsequent workers. 

It is of interest to note that the limiting heat inflow of 1 erg/min. which was 
finally obtained is about 60 times less than the conduction along the lead wire 
calculated by extrapolation of the measurements of de Haas and Rademakers (1940) 
according to a T? law. This discrepancy may be partly due to the fact that the 
specimen used by these authors was a single crystal of relatively large diameter 
(2-5 mm.). 

Some rough measurements were made of the heat inflow to the lower stage 
when a field greater than the threshold value was applied to the lead wire by means 
of a small permanent magnet. With the upper stage at 0-3°K. and the lower 
stage between 0-01° and 0-1° heating rates of the lower stage of more than a 
hundred ergs per minute were observed. This is consistent with the expectation 
that the heat flow should be increased approximately as (T,/T,)?. 
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$82 POSSIBILITY OF CYCLICAL PROCESS Es 


The use of superconducting wires for heat transfer in a cyclical magnetic 
cooling process has been proposed by Daunt and Heer (1949). It appears to us 
that some difficulty would be experienced in attempting to operate such a process 
with a starting temperature as high as 1°K. because the superconductor is not a 
very efficient thermal switch at this temperature. The main difficulty is that 
during the period of magnetization, when the salt is at a temperature above 1°K., 
a considerable amount of heat is transmitted along the second wire even though 
it is in the superconducting state. 

It would, however, seem practicable to operate a cyclical process starting 
from a lower temperature, say in the region of 0-2°K., using as a heat reservoir a 
large quantity of a suitable paramagnetic salt such as iron alum. In this way it 
should be possible, with a single stage cyclical process, to maintain a temperature 
of about a hundredth of a degree for several hours provided the heat influx from 
external sources were not more than a few ergs per minute. 

Although our apparatus was not designed for cyclical operation, we have 
found it possible to carry out several successive demagnetizations of the second 
stage without remagnetizing the first stage. It seems that an apparatus of the 
kind we have described could, with minor modifications, be adapted for operation 
as a cyclical system. 


§9. APPLICATION TO NUCLEAR DEMAGNETIZATION 


It has been suggested (Gorter 1934, Kurti and Simon 1935, Simon 1939) 
that it should be possible to extend the lower limit of temperature attainable 
by the use of nuclear paramagnetism instead of electron paramagnetism. Since 
nuclear magnetic moments are of the order of one thousandth of the electron 
moment, it is impossible to remove an appreciable amount of nuclear spin entropy 
at 1°K. by the application of any practically realizable magnetic field. The 
magnetization must be carried out at a temperature lower than 0-01°K., so that 
the use of an initial stage of electronic demagnetization is essential. 

It seems that the superconducting switch would be well adapted for controlling 
the flow of heat between the nuclear stage and the electronic stage. At such low 
temperatures an extremely high on-off ratio for the superconducting switch 
should be obtainable. One practical difficulty is that, in order to remove the 
heat of nuclear magnetization in a reasonable time, a large cross section of metal, 
of the order of 1cm?, would be required. Some difficulty might also be 
experienced in establishing adequate thermal contact between the metal of the 
switch and the paramagnetic salt of the first stage. 

Even at a temperature of 0-01°K. a magnetic field of the order of 50 kilogauss 
would be necessary to remove an appreciable fraction of the nuclear spin entropy. 
It does not seem practicable to perform a satisfactory nuclear demagnetization 
by using smaller magnetic fields and more stages of electronic demagnetization 
to obtain a lower temperature for the nuclear magnetization process. Owing to 
unavoidable heat inflows, the nuclear stage would need to be an appreciable 
fraction of a mole in order to remain at the low temperature long enough to 
enable observations to be made. The requirement that each stage must contain 
a greater molar fraction than the subsequent one in order to take up the necessary 
entropy without appreciable rise of temperature means that the first stage 
would have to be of impracticably large dimensions, 
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Note added in proof. In more recent experiments with an apparatus similar 
to that described above, we have used a water-cooled electromagnet giving 
a field of about 9,000 gauss. The Curie temperature reached by the second 
stage, which contained 8-4 x 10-4 mole of Cr ions, was 0:0045°. By extrapolation 
of the results of de Klerk, Steenland and Gorter (1950), this corresponds to an 
absolute temperature of approximately 10-°°x. 

The temperature remained below 0-01°K. for 40 minutes. 
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A Note on Frohlich’s Theory of Superconductivity 
By KUN HUANG 


Department of Theoretical Physics, University of Liverpool 
Communicated by H. Frohlich; MS. received 21st February 1951 


ABSTRACT. It is shown that the physical mechanisms underlying the two energy terms 
E, and E, in Fréhlich’s theory are very different: EF, is largely the energy change due to a 
modification of the lattice vibration frequencies caused by the adiabatic deformation of the 
electrons; E, is on the other hand a dynamical term, representing the energy change due to 
virtual collisions between lattice oscillations and electrons. A proof is given that certain 
inadmissible consequences of the theory, pointed out recently by Wentzel, derive from the 
energy term E£,; the superconductive behaviour discussed by Fréhlich arises however 
entirely from the term E,. Although both terms follow from the same perturbation treat- 
ment, the analysis given makes it seem not unlikely that whatever mechanism actually 
inhibits EZ, will leave the superconductive behaviour deduced from E, not substantially 
affected. 
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Frohlich’s theory (Frohlich 1950) which have been obtained in the course of 

the investigations undertaken by the author in active collaboration with 
Mr. M. J. Buckingham and Mr. B. V. Paranjape of this department; they are such 
as to give a better understanding of the physical mechanisms underlying the 
different energy terms in Fréhlich’s theory. A clarification of such issues is 
deemed particularly desirable, for they are closely related to certain points raised 
recently by Wentzel (1951)*. Wentzel’s conclusions can be summarized as follows: 
(a) the interaction term between the lattice and electrons required for super- 
conductivity is so strong that the convergence of the perturbation method used 
in Froéhlich’s theory is bad; (6) the interaction is of such a magnitude as to 
cause in some cases an apparent instability in the lattice, which must in practice be 
prevented by other stabilizing influences such as the Coulomb interactions; these 
influences should be taken into account in a realistic treatment. 

To save space we shall use the original notation in Frohlich’s paper (Frohlich 
1950) without further explanation ; hereafter his paper will be referred to as (F). 

Frohlich’s treatment refers throughout to the case that none of the lattice 
vibrations are excited. A significant difference between the energy terms denoted 
by Frohlich as £, and EF, emerges, however, when one considers an arbitrary 
excited state of the lattice. Let us consider the part of these energies contributed 
by a particular mode w, which is in the quantum state v. We have now to 
consider the energy terms due to the virtual absorption as well as emission of a 
phonon. Consider first the terms corresponding to F,. The emission terms can 
be taken directly from (2.14) of (F) by omitting the summation over w and 
multiplying by (v+1): 

16Flo,(v ar 1) | Sic tu-wt 
3n(27)?>_ J |k-—w|?--k? + a9w 

where we have expressed q in (F) explicitly as k—w. The absorption terms can 
be obtained from (1) by replacing (v+ 1) by v, k—w by k+w and apw in the 
denominator by —o pw: 


3n(2n)* J [k-+ wk? — oo dk (absorptionterms). ...... (2) 


Let us consider the integral in (2). If we express k in terms of q=k+w, the 
i lb 
integral becomes j wfawha 

g —|q—wl)?— aw 


Ip is the purpose of the present note to describe certain conclusions as regards 


dk (emissionterms) ...... (1) 


dq, 


which is clearly the negative of the integral in (1). Thus, adding (1) and (2),we 
obtain the total contribution to £, due to both virtual absorptions and emissions as 
16Floy (fifi 
3n(277)8 |k—w/|?—k? + oqw Ce OE acre ge (3) 
The point to note is that the expression is independent of the lattice quantum 
number w. 

E,, on the other hand, is not independent of the lattice quantum number. The 
lattice instability as discussed by Wentzel means simply that the energy of the 
system can be lowered indefinitely by increasing indefinitely the lattice quantum 
number. ‘This situation can only be caused by terms dependent on the lattice 
quantum number; Wentzel’s considerations thus concern only the E, term. 


* Tam grateful to Professor Frohlich for access to Professor Wentzel’s work before its publication. 
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The contribution to E, due to a particular mode w in state v, when worked out 
explicitly, is given by 
Flo, . 4 
3p saaae +1) {Eh —(w+t op) JIn\5 
Flogv A 2ky —og + 
3n(2m)e {4 = (w = Go) ] In pI re ae 4k,(w = Cees 
where the two lines correspond respectively to the emission and the absorption 
terms. : 
For lattice waves long compared with the lattice constant, we expand (4) with 
respect to w and ignore second and higher order terms: 
4FCk oa. = 8 FL Ryo 
— Faldnye ~ 3n(Ony? (Zo NW) g0 0 a hie eheroak as (5) 
The excitation energy of the lattice vibration in state v is given on the other hand by 


CIOS meek st ed eh LP) pesetiin (6) 
On combining (5) and (6), we find 
4FCRoo" . 16FCRooy 

PaaS ' ~ 3a(anyh } 
which, we note, includes all terms dependent on v. When the coefficient of the 
second term is positive, the term is equivalent to a reduction of the sound velocity 
from the value s to the value” 
16FCkyoy 
3n(27)?h 
If however (8) is negative, the energy of the system can be indefinitely reduced by 
increasingv. Thisis the instability discussed by Wentzel. 

To obtain a better understanding of this change of sound velocity let us 
suppose as a pure hypothesis that during a lattice vibration the electrons are merely 
deformed adiabatically. At every instant, we would thus imagine the electrons as 
being subject to a static perturbation corresponding to the instantaneous 
configuration of the lattice. According to the adiabatic approximation, the corre- 
sponding electronic energy operates as an additional potential function for the 
lattice motion. We wish to determine this potential function and find out how it 
modifies the lattice motion. 

Asin(F), we have so far been using the normal coordinates introduced originally 
by Peierls. In Peierls’ scheme (cf. Sommerfeld and Bethe 1933), the lattice 
displacement is given by 


+4k (00+ 00) | i 
ere. 


2Ro pete 
99 


= Cb aan Pee > bee (8) 


S— 


yy 2 ew ite exp (iw.x)+a,*exp(—iw.x)} (e€,=e _y) 

where e,, is the polarization vector and ay, ay * are respectively the annihilation and 
creation operators for the normal mode w. ‘These operators depend on the 
momenta as well as the normal coordinates. ‘The lattice—-electron interaction, 
being constructed from u(x), depends thus also on the momenta. ‘This fact 
makes it inconvenient to use Peierls’ coordinates to consider the adiabatic 
modification of lattice motion. Thus instead we shall make use of the more usual 
normal coordinates g,, which are a to the normalized Fourier coefficients 


O,, of u(x): 


u(x) = 


u(x) = aH ie Xi €@wOw eXp (7W. X) 
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by Oy =O_w* =(dwttq_w)/\/2, W in this relation being restricted to one side 
of a plane through w=0. 

Consider now a static deformation due to gy alone; the corresponding lattice 
displacement is given by 


1 1 
u(x) = VN ord 75 {exp (tw. x)+exp(—zw.x)}. Matte (9) 
Let us contrast it with the expression in Peierl’s scheme : 


u(x) Cw {ay exp (tw. x) +a,,* exp (—7w.x)} 4 ema: (10) 


el 
Say 
which underlies our earlier considerations of the contributions to E, and F, due to 
a single Peierls mode w. We note that (10) would reduce to (9) if we replace both 
Gy and Ay* by qw/1/2. This fact makes it possible to obtain the electronic energy 
(in the approximation of a second order perturbation calculation) due to the static 
deformation (9), from the expressions for £, and E, by suitable replacements. The 
second order terms due to the static perturbation have a strict correspondence with 
the terms in £,+E,. Now as the normal coordinate q, is to be considered as a 
parameter and we are only concerned with virtual electronic transitions, we have to 
delete the o)w term in the denominators of (2.13), (2.14) in(F). Moreover, in the 
absorption and emission terms of E, and E, we have incorporated the squared 
transition matrix elements of a, and a,*: 


h fy 
(om) Nous uae (11) 
In view of (9) and (10), for the static perturbation, both of these squared matrix 
elements should be replaced by 
rive ee eer ee (12) 


Consider first the static perturbation terms analogous to £;. Putting o,=0 
in the denominator of (2.13) of (F) is equivalent to putting o) =0 in the brackets in 
(4); and the replacement of (11) by (12) is equivalent to replacing both (v+1) and 
v in (4) by Mszwg,,”/h. ‘Thus we obtain from (4) the corresponding expression for 
the static perturbation as 


2M Flo sw 


2k, +w 
im 3n(2n)2h OP { (Hh? —w*)In 


| +4hye}. 


Expanding for long waves, we find 


16F Eko ne scree 
— Ful Qn) Msw*q,? = —2vFMs*wq,7. 
If we add this additional potential function due to the adiabatic deformation of the 
electronic state to the potential function of the mode q,, viz. $M(ws)? qy?, we 
obtain the modified harmonic potential : 


4M {w?s*(1 —2vF) qq, ae. (13) 


where we have completed the square by adding a term in F?, which is formally 
permissible in view of the perturbation procedure we have used. ‘The adiabatic 
consideration leads thus to exactly the same change in the sound velocity as (8). 
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- This shows that at least for the longer waves E, describes merely an adiabatic 
deformation of the electrons by the lattice vibrations. 

We have not yet considered the static perturbation terms analogous to £5. 
Such terms however vanish identically. For we have to replace v+ 1 and vw both 
by the same factor (Msw/h)q,,” in (1) and (2) and add the two expressions. 'Thesum 
vanishes, as we have shown that the two integrals are equal and opposite in sign. 
The £, term is thus not involved in the adiabatic deformation of the electronic 
state. 

The difference between the underlying physical mechanisms of the F, and EF, 
terms in Frohlich’s theory is now clear: FE is mainly the energy change due to the 
modifications in the lattice frequencies, caused by the adiabatic deformation of the 
electrons. £, is on the other hand a purely dynamical effect, representing the 
energy change of the electrons due to virtual collisions with the lattice vibrations. 

From a purely practical viewpoint, while the large lattice frequency changes 
following from F, are almost certain to be inhibited by some mechanism, £, is not 
to be similarly suppressed. For the latter is clearly related closely to the electron 
scattering responsible for electrical resistivity. In view of the essential difference 
between the mechanisms responsible for the two energy terms as well as the fact 
that £, depends on the lattice vibration quantum number, whereas F, does not, one 
is thus inclined to the belief that whatever mechanism ultimately inhibits £, will 
not seriously affect F,. In this connection, one may note the heuristic advantage 
of the use of the experimental resistivity values in the numerical estimates given in 
Frohlich’s paper. 

A discussion of the effect of the Coulomb interaction between the electrons is 
given at the end of the Appendix. Effects due to the correlation between electrons 
are still being investigated. 

It need hardly be emphasized that this note does not provide an unequivocal 
answer to the points raised by Wentzel. But the analysis given should make us 
better able to assess the situation and may serve as a pointer for further 
investigations. 

It is a pleasure for me to express my thanks to Professor H. Frohlich and 
Mr. M. J. Buckingham for stimulating discussions. | My thanks are also due to 
Mr. B. V. Paranjape who has undertaken some calculations in connection with 
the investigation. 


tig) 2p ea OO a Oe bd < 


THE ADIABATIC DEFORMATION OF ELECTRONS 
BY LATTICE VIBRATIONS 


For his considerations, Wentzel has employed Tomonaga’s boson treatment 
(Tomonaga 1950) of the electronic motion, which holds for a linear case if the 
electron density does not vary appreciably within a length of the order of the 
lattice constant (or average electron distance, to be more precise). We wish to 
point out briefly the simple underlying physical picture, as it further illuminates 
the above interpretation of £,. However, we shall omit all mathematical details 
and report only on the results of our investigations. 

It is found that ‘Tomonaga’s method is completely equivalent to the use of 
the Thomas—Fermi method in each small neighbourhood. That he can use it for 
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dynamical problems is the exact parallel of the usual practice in classical continuum 
physics. If one takes for instance his density waves for free electrons, one finds 
that the frequencies are exactly reproduced by the pressure waves corresponding to 
the compressibility of a free Fermi gas (in the completely degenerate state). 

When the interaction between the electrons and the lattice is introduced, two 
types of harmonic vibrations appear in the Tomonaga scheme corresponding to 
lattice vibrations perturbed by electrons, and to electron pressure waves perturbed 
by the lattice. When we examine the former type of solutions closely, the following 
simple picture emerges: When a lattice wave passes, we can think of the bottom of 
the potential well enclosing the electrons as being warped (solid line along the 
bottom of Figure). If the electron density is kept constant as before (top of the 


Adiabatic 
SS = Distribution 


\ <| Constant 
Density 


distribution indicated by dotted line in Figure), obviously no net energy change 
results from the interaction, for half of the electrons have their energies raised, and 
the other half have their energies lowered. ‘The reduction of the lattice frequencies 
in Tomonaga’s method arises from the circumstance that the electrons flow from 
the energy ridges to the energy troughs, so that the top of the distribution is kept 
level (solid line in Figure). In the sense of the Thomas—Fermi method, the 
electrons follow the lattice motions so that at every instant the electrons occupy the 
lowest states in the instantaneous lattice potential; thus the Tomonaga solution 
corresponds strictly to an adiabatic treatment of the electrons. As shown by 
Wentzel, it leads to the possibility of instability ; it is clearly equivalent to the FE, 
term which we have discussed. 

That such considerations should lead to a distinct type of stability condition is 
only to be expected. For it is associated with the flow of electrons from one 
volume element to another, a situation not to be realized in a static and homo- 
geneous elastic deformation. Clearly such instability owing to electron flow will 
in reality be prevented by space-charge effects. 

It is a straightforward matter to generalize the stability condition as following 
from. Tomonaga’s method to the three-dimensional case. We need only observe 
that the one-dimensional result can be realized in an extremely simple manner: 
we consider two equal pieces of metal, one is subject to a uniform dilation, the other 
to a compression, while free passage of electrons is allowed. Instability occurs in 
this ‘ two-piece model ’ if the lowering of energy due to electron flow exceeds the 
elastic strain energy in the lattice due to the compression and dilation. We have 
introduced an explicit lattice-electron interaction by assuming that the change in 
lattice potential for electrons is at every point proportional to the linear strain. 
With this explicit interaction, we have verified that the instability condition in the 
“ two-piece model ’ is the same as the condition for vanishing lattice frequencies in 
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the Tomonaga treatment. The three-dimensional analogue follows thus directly 
by replacing the linear metals in the ‘ two-piece model’ by two boxes. 

Using Tomonaga’s method, Wentzel has found that if the Coulomb interaction 
between the electrons is taken into account, the instability due to long waves is 
prevented. The three-dimensional case is similar: Mr. Paranjape has investigated 
this case by considering a density wave through the lattice and treating the 
electrons by the Thomas—Fermi method. He finds that the electron density 
bears to the ion density the constant ratio (— C + 47e?n,\/w?)/(26/3 + 47€2n,)/w?). 
It is immediately verified that if the wavelength 27/w is large compared with the 
lattice constant, the terms with e? predominate and the ratio is practically equal to 
unity. In other words, in the longer waves the Coulomb interaction prevents the 
electron flow as envisaged in the Figure and the adiabatic energy change is thus 
completely suppressed. 

Concerning the effect of the Coulomb interaction on the energy F,, we may cite 
the calculation for the resistivity of some metals made a long time ago by Bardeen 
(1937), who has used the self-consistent field of the electrons as the scattering field. 
His treatment thus incorporates the effect of the adiabatic deformation of the 
electrons whereby the Coulomb interaction is taken into account. He found 
resistivity values smaller, but not drastically so, as compared with results obtained 
by conventional methods; the latter, we note, are comparable with the method 
used in deriving E,. From this we may infer that the Coulomb interaction, while 
preventing the instability of the longer waves, will not drastically affect the value 
of E,. 

From (4) we readily find that waves with w comparable with hy are also subject 
to instability, for values of F somewhat larger than the critical values for the longer 
waves. We are investigating whether the exchange effect in a treatment such as 
that used in Bardeen’s resistivity calculation is capable of suppressing the instability 
of these shorter waves. 
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Scintillations from Organic Crystals: Specific Fluorescence 
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ABSTRACT. The scintillation response S of organic crystals depends on the nature and 
energy FE of the incident ionizing particle, of residual range 7. The specific fluorescence 
dS/dr is not in general proportional to the specific energy loss dE/dr. By considering the 
quenching effect of the molecules damaged by the particle on the ‘ excitons ’ produced by it, 
it is shown that dS/dr=(A dE/dr)/(1+kB dE/dr). Aand kB are constants, which have been 
evaluated for anthracene from observations of SS and L, and the range-energy data. Curves 
are computed for the relative response S of anthracene to electrons, protons, deuterons and 
a-particles of HE up to 15 mey., and these are shown to agree closely with the available 
experimental results. The method used for evaluating the relative response is applicable 
to ionizing particles of any nature or energy, and also to the other organic scintillation 
crystals. 


Se UN ROD Wi Cals Omn 

ONIZING radiations impinging on a fluorescent material produce short in- 
dividual light flashes, or scintillations. These scintillations can be detected 
with a photo-multiplier tube, and converted into electrical pulses, which can 
be counted and measured by standard electronic methods. ‘This technique of 
scintillation counting is being widely applied for the detection and measurement of 
nuclear radiations. The fluorescent organic crystals are of particular interest for 
scintillation counting, since they combine a reasonable fluorescent efficiency with a 
high transparency and a very short luminescent decay time, of the order of 10-® 
second. ‘These properties make them suitable for the detection of the more 
penetrating nuclear radiations, and for studies of fast nuclear and meson decay 
processes, and much of the previous work in this field has been primarily concerned 
with such applications. The nature of the scintillation process has received 
rather less attention, and it is, therefore, proposed in this and subsequent papers to 
consider various fundamental aspects of the fluorescence produced in organic 

crystals by ionizing radiations. 


§2. RESPONSE TO DIFFERENT RADIATIONS 


The intensity of the scintillations produced in anthracene and other organic 
fluorescent crystals depends both on the energy and on the nature of the incident 
ionizing particle. ‘The amplitude S (volts) of the voltage pulse from a photo- 
multiplier, operating under constant conditions and observing the crystal, is 
proportional to the number of fluorescent quanta produced, and hence S may be 
used as an arbitrary measure of the scintillation intensity. For electrons of energy 
greater than 125 kev., the scintillation intensity S from an anthracene crystal 
increases linearly with the energy E (Hopkins 1951), so that the fluorescent 
efficiency dS/dE (volts/Mev.) is constant. For electrons of lower energy, however, 


* Now at the Physics Department, Rhodes University, Grahamstown, South Africa. 
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the efficiency dS/dE is reduced. With heavier particles, which produce more 
intense ionization, a further reduction in dS/dE occurs, and the scintillation 
intensity is consequently less than that produced by an electron of the same 
energy. 

It has been noted previously (Birks 1950 a) that the response to different 
ionizing particles can be readily compared by considering the variation of the 
specific fluorescence dS/dr (v/cm. air equivalent), which is proportional to the 
number of fluorescent quanta emitted per unit path length, with the specific 
energy loss dE/dr (Mev/cm. air equivalent), where 7 (cm. air equivalent) is the 
residual range of the particle. At low values of dE/dr (electrons of H>125 kev.) 
the specific fluorescence is proportional to the specific energy loss, corresponding to 
S increasing linearly with E. At high values of dE/dr («-particles of <5 Mev.) 
‘dS/dr is practically constant and independent of the specific energy loss. The 
scintillation intensity S is thus proportional to the residual range 7 of the particle, 
rather than to its energy E. 

§3. THEORY 

The variation of dS/dr with dE/dr may be readily accounted for on the 
‘ exciton ’ theory introduced by Bowen et al. (1949) to account for the fluorescence 
of mixed crystals, excited by ultra-violet radiation. On this theory of the energy 
transfer processes in organic crystals, the electronic energy excited by the incident 
radiation is transferred from molecule to molecule within the crystal, until it is 
captured by a single molecule, which then either fluoresces or quenches the 
‘exciton’, depending on the nature of the molecule. This theory has been 
successfully applied to the energy transfer in mixed crystals of naphthalene 
and anthracene excited by «-radiation, where both molecular components 
fluoresce, but the anthracene molecules have the higher ‘exciton’ capture 
probability (Birks 1950b). It has also been used to explain the deterioration of 
the fluorescent efficiency of anthracene under prolonged «-particle irradiation 
(Birks and Black 1951). _‘In this case, the molecules damaged by the preceding 
«-particle irradiation act as quenching agents, having a higher ‘ exciton’ capture 
probability than the undamaged anthracene molecules. 

A similar effect will also occur with individual ionizing particles. ‘The passage 
of the particle through the crystal produces a local concentration of damaged or 
ionized molecules along its path. ‘These damaged molecules act as quenching 
agents for the ‘ excitons’ produced by the ionizing particle. ‘The number of 
‘excitons’ produced per unit path length is proportional to the specific energy loss, 
say AdE/dr. he local concentration of damaged molecules at any point on the 
particle track is also proportional to the specific energy loss, say B dE/dr molecules 
per undamaged molecule. Hence, if the ‘exciton’ capture probability of a 
damaged molecule, relative to an undamaged molecule, is k, the specific fluorescence 


aS AdE/dr 1 
OS HUD eel ie ee ii. We aa ie) 


At low values of dE/dr, (1) becomes dS/dr =A dE/dr and at high values of dE/dr, 
(1) becomes dS/dr = A/RB = constant, in agreement with the observed behaviour. 

The values of A and kB for anthracene have been computed from observations 
on the relative response to electrons and «-particles reported previously (Birks 
1950a). The residual ranger is taken from the range—energy data of Curie (1935) 
for electrons, and of Livingston and Bethe (1937) and Bethe (1950) for heavier 
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particles. 7 is expressed in cm. air equivalent at 15°c. and 760 mm. Hg pressure, 
but the small correction for the variation of the relative stopping power of anthracene 
with has been neglected.- This variation is approximately + 1% forr=0-5 cm., 


Mablewl 
Particles E(mev.) dE) dr dS/dr 
a-particles 2-0 1:8 10-7 
Electrons 0-5 0-00215 0-175 


and +2% for r=5-180cm. The specific energy loss dE/dr is expressed in 
Mev/cm. air equivalent. The experimental values used for substitution in (1) are 


listed in Table 1, giving for anthracene A = 82-5, RB=7-15. 


a-particles 
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The variation of dS/dr with dE/dr from (1) is plotted in Figure 1. The specific 
fluorescence dS/dr for different ionizing particles of residual range 7 has been 
calculated from (1) and the range-energy data. The results are plotted in Figure 
2. The variation of S with r has been obtained by integration of the area below 
the (dS/dr, r) curves. The scintillation intensity S produced by different 
ionizing particles of energy £ is plotted in Figure 3 and the response to low energy 
electrons is plotted on an enlarged scale in Figure 4. 


$4. COMPARISON WITH EXPERIMENTAL DATA 

The scintillation response of anthracene to protons of E<16-4 Mev. has been 
observed by Franzen et al. (1950), to electrons of E <3-0 Mev. by Hopkins (1951), 
and to «-particles of E <5:3 Mev. by Birks (1950a). In each case the observations 
_ agree closely with the predicted responses shown in Figures 3 and 4. 

In addition the relative responses of anthracene to different particles have also 
been reported. A comparison between the experimental response ratios and the 
predicted values is made in Table 2. The agreement is sufficiently good to estab- 
lish the general validity of the theory, which may be used to compute the response 
to any particle. 


Table 2 
Reference Particles Exptl. ratio Theor. ratio 
Harding et al. 8-78 Mev. a-particles iG oe 
(1949) 6-04 Mev. a-particles 
De . O- i 
[ 3 Mev. a-particles | 0-82 0-79 
Franzen et al. 0-624 mev. electrons 
(1950) 16:4 Mev. protons 16-4 1755 
{ 0°624 Mev. electrons 
. prot 
Hofstadter (1950) Se CN 2-0 2-16 


5 Mev. electrons 


5:3 Mev. electrons 
Hopkins ( ) 5:3 Mev. a-particles i 


Although the present discussion has been confined to anthracene, for which 
most experimental data are available, the response curves for other organic crystals 
to different ionizing particles will be very similar. It is known that the ratio of 
the response to 5-Mev. «-particles and electrons is approximately the same for 
anthracene, naphthalene (Birks 1950 a) and terphenyl (Hofstadter 1950), and 
similar results are to be expected for stilbene, phenanthrene, dibenzyl and the other 
fluorescent hydrocarbons. 
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One-Dimensional Lattice Model 
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ABSTRACT. The effect of changing the lattice parameter in Kronig and Penney’s linear 
lattice model is investigated theoretically. Méglich and Rompe have predicted a decrease 
in the width of the forbidden energy zone as the lattice parameter is increased. Itis shown ~ 
that this is true only if the width of the potential troughs is increased whilst the thickness of 
the barriers is maintained constant. However, the opposite effect, i.e. an increase in the 
width of the forbidden zone, is expected when the potential troughs are kept constant and 
the width of the barriers is increased instead. In the intermediate case of a uniform expan- 
sion of potential hills and troughs, a small net decrease in the width of the forbidden zone is 
expected. 


semiconductors have recently been investigated experimentally and efforts 

are being made to account for them on the current band theory of solids. 
Changes in temperature and pressure have been found to affect the fundamental 
lattice absorption (optical branch) and particularly the intrinsic conductivity of 
semiconductors. ‘The long-wave edge of the optical or ultra-violet absorption 
band of many insulators is shifted towards the red as the temperature is increased 
(by about 1 4/°c.) and conversely towards the blue on the application of a high 
pressure, (in CdS crystals for instance by about 0-02 a/atm. (Hohler 1949)), 

One cause of these effects is the change in inter-atomic distances produced by 
thermal expansion or by mechanical compression. ‘This change in lattice 
parameter has been considered by Méglich and Rompe (J 942) in an effort to calcu- 
late the magnitude of the temperature effect. Using the well-known linear model 
of Kronig and Penney (1931), they arrive at the conclusion that the width of the 
forbidden energy band which determines the absorption edge decreases as the 
lattice parameter is increased, the percentage decrease being equal to twice the 
percentage increase in lattice parameter. ‘This result seems rather surprising; 
it might be expected that at large inter-atomic distances such a lattice would behave 
as an insulator with discrete energy levels and that, as the atoms are brought more 
closely together, these levels would gradually broaden into bands, i.e. that the 
energy gap would gradually be reduced. It is the purpose of this note to re- 
examine the problem, as Moglich and Rompe’s conclusions continue to be used by 
subsequent investigators (Hohler 1949, Seiwert 1949). 

It will be recalled that Kronig’s model consisted of the periodic arrangement in 
one dimension of potential troughs (of width a and potential V =0), alternating 
with potential barriers (of width 5 and height Vj) corresponding to the barriers 
which impede the free passage of electrons from one atom. tothe next. Solution of 
Schrédinger’s equation with these boundary conditions yielded the following 
relation for the permitted energy levels (eigenvalues W) of the electron : 


y — fox 
2By 


T HE effects of various factors on the electrical and optical properties of intrinsic 


sinh yb sin Ba+coshybcosBa=cosa(a+b), ss see (1) 
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where the wave vector B=CW1?, and y=C(V,—W)!?, C? =87?m/h?, and « 
could assume any real value. Equation (1) yielded real values of W only for 
certain ranges of «: only if the absolute value of the expression on the left does not 
exceed unity does the solution correspond to a permitted (real) energy value and 
the electrons are in an allowed band; the alternative of imaginary energy values 
corresponds to electrons in a forbidden zone. Kronig andPenney discussed their 
solution under the following two simplifying assumptions: (i) B<y, i.e. the dis- 
cussion is restricted to those energy values of the electrons which are well below the 
top of the barrier (W<YV,); transition of electrons occurs by tunnelling through 
the potential barrier rather than by passing over it. (ii) by <1, i.e. the value of the 
parameter P=yab/2 = Vb aC?/2, which is proportional to the product of barrier 
height and width (the ‘ barrier resistance ’), remains finite. 
Under these assumptions, equation (1) reduced to 


(Psin 6a)jBo--cos pa =cosads =a = ree (2) 


Weare here interested in insulators and intrinsic semiconductors, particularly 
in the energy values corresponding to the upper edge of the first (full) band and the 
lower edge of the second (empty) band, their difference corresponding to the width 
of the forbidden zone. At these edges the right-hand sides of equations (1) and (2) 
assume the extreme values +1. Equation (2) is then fulfilled if either B,a=7, 
corresponding to the top of the first (full) band, or if 


PRacoL (Bea) 2) = 2 Paw Gt PU eee ese (3) 


corresponding to the bottom of the second (empty) band, the difference in the 
respective energy levels, W=?/C?, determining the width of the forbidden gap 
W,. 
The effect of altering the width of the potential trough a is found by differen- 
tiating W, with respect to a, whilst maintaining the barrier height Vy and width 6, 
1.e. the barrier resistance, constant. This always results in a negative value, 
dW ,,/da <0 (Figure 1(5)); as the space available to an electron between two adjacent 
barriers is increased, so the forbidden band becomes narrower, the model assumes 
a more and more metallic character. This is essentially Moglich and Rompe’s 
result, except that their actual formula considerably over-estimates the change in 
gap width. * 

Instead of maintaining b constant and varying a, we might keep a constant and 
vary b. ‘This can be done whilst still maintaining the useful approximations B <y 
and by<1.+ The gap is then found to increase with increasing b (Figure 1(c)); as 
the barrier becomes wider, the probability of an electron tunnelling across decreases 
and the situation in the limit of large separation approaches the case of tight 
binding, i.e. of discrete energy levels. 

The intermediate case of uniform strain along the lattice, da/a=65/b, leads to 
the width of the gap always decreasing with increasing lattice parameter, as in the 
case considered first, but the decrease is considerably less. 

In conclusion, expansion of the lattice may, according to this model, lead to 
either a narrowing or a widening of the forbidden zone, the former if the width of 
the troughs is increased and the latter if the width of the potential barriers is 


* Méglich and Rompe have not allowed for the variation of the cotangent term when in effect 
differentiating equation (3). 
+ The barrier resistance (Vb or P) does not, of course, remain constant in this case. 
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increased. It cannot easily be decided which alternative is the correct one; 
Kronig and Penney’s model is in fact too simple to be used for quantitative 
calculations as it does not take into account any finer details of the atoms making up 
the lattice. More quantitative results are obtained by using the appropriate 
atomic or ionic potential functions instead of square barriers. A good approxi- 
mation for the proper wave function is obtained by a linear combination of the 
atomic orbitals of lowest energy (I.c.a.o. method) in such a way that continuity of 
the wave function and of its first derivative is maintained when going from any one 
atom to its nearest neighbour (Wigner and Seitz 1933). This treatment is equally 
suitable for one-, two- or three-dimensional lattices and yields detailed results for 
the electronic band structure of metallic or covalent solids which can be used for 
comparison with experiment. The result of such a calculation for diamond is — 
shown in Figure 1(a), giving the energy bands as a function of inter-nuclear 


! 
| 
(a) 
| 


Figure 1. Electronic band structure of covalent solids: (a) l.c.a.o. method (Kimball 1935) ; 
(b) square barriers of constant thickness b ; (c) square barriers between potential troughs of 
constant width a. 


distance (Kimball 1935). It is interesting to note that according to this structure 
also—which is typical for many solids—either an increase or a decrease in gap width 
may arise on expansion of the lattice. For the covalent lattices calculated so far — 
(C, Si) the equilibrium distance is less than the spacing corresponding to the 
crossover point, Rj. Theory thus predicts a decrease of the energy gap with 
expansion of these lattices, in agreement with experiment. 
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ABSTRACT. For a black body inserted into a homogeneous scattering medium, the 
neutron density at large distances from the body is determined using the perturbation 
method for the case when the scattering in the scattering medium is assumed elastic and 
isotropic capture is absent and the following four types of the black body are considered : 
(i) a sphere of radius small compared with the mean free path in the scattering medium; 
(ii) a sphere of radius large compared with the mean free path; (iii) an infinitely long cylinder 
of radius small compared with the mean free path; (iv) an infinitely long cylinder of radius 
large compared with the mean free path. Combining the results for the above limiting 
cases with the results obtained by other methods the distribution is also estimated for a 
sphere and an infinitely long cylinder of intermediate radii. 


INTRODUCTION 


HE integral equation for the neutron density in a non-capturing isotropically 
scattering medium occupying a convex volume I, in the presence of black 
bodies inside this volume, is (Peierls 1939) 


eis n(r’ 
n(r) = ral | | ee = eR AV my eee (1) 
provided that we are dealing with the stationary case and that the mean free path 
of the neutrons is independent of energy. Here all distances are measured in 
units of mean free path, R is the distance between 7 and 7’, V(r) is the part of the 
volume V which is not shielded from 7 by the black bodies and n(r) is the neutron 
track length per unit volume per unit time in the vicinity of r. For brevity we 
shall refer to n(r) as the neutron density, though the use of this latter term is not 
truly justifiable, unless we assume that neutrons are scattered without change in 
energy. 1,(7) is the contribution to n(r) due to neutrons which have come directly 
from the sources or directly from outside the volume V, without having been 
scattered in V. 

It is known that ir the source free regions n(r) can be presented in the following 


form 
t Al =D ae ee Neg 5 a laa ee (2) 
in which mig (7) satisfies the diffusion differential equation 
V2 Natt (r) == ) see eee (3) 


while 1*(r) is of the order of e°, 5 being the distance from the point 7 to the 
nearest source, or the nearest point on the outer boundary of V, or on the surface of 
the nearest black body, whichever is the smaller. If the dimensions of V and the 
distances of the black bodies from the outer surface of V and between each other 
are large, it is usually sufficient to know ngjq(r) rather than n(r) itself. Equation (3) 
;S not sufficient in itself to determine 7,;_(7), one should know also the boundary 


t This could be seen by slightly extending the analysis given by Bothe (1942). 
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conditions satisfied by the latter at the boundaries of the regions containing 
sources, at the outer boundary of the volume V and at the surfaces of the black 
bodies. In the present paper we shall be concerned only with the boundary 
conditions at the surface of a black body. To specify these boundary conditions — 
it is sufficient to determine the ratio 


ney / “ae Par (4) 


(in which vis the direction of the outer normal) at the surface of the black body. _ If 
the black body is situated at a large distance from all sources, other black bodies 
and the outer boundary of the volume V, then the values of (4) should depend 
almost entirely upon the size and shape of our black body. 

We can determine an approximate value of (4) by solving the idealized problem 
in which the scattering medium. extends to infinity in all directions, containing 
no other black bodies apart from the one for which we wish to determine 
the boundary conditions, and the neutrons are assumed to be supplied from 


infinity uniformly from all directions. In this idealized problem ngjg(7) gives, 


at the same time, the asymptotic behaviour of m(r) at large distances from the black — 
body, and accordingly we shall denote it henceforth as ,,(7). ‘The value of the 
quantity (4) can be referred to as the linear extrapolation length, we shall denote it 
by A. 

The values of 4 for a sphere have been calculated by Marshak (1947 a, b) fora 
wide range of values of the radius using numerical approximation methods. It — 
has also been shown by Placzek (private communication) that for a convex black 
body of a fixed shape A should tend to the value 4/3 as the size of the black body 


tends to zero, while as the size of the black body becomes infinite the problem 


becomes the Milne one, for which A =0-7104 (see, for example, Hopf 1934, p. 85). 
It has been conjectured further that A varies monotonically with size of the black 
body. ‘This suggests that one might try to. determine expansions for X, valid for 
small and large sizes. ‘These expansions are of interest also for the spherical 
case, despite Marshak’s work. For a small sphere, in particular, the spherical 
harmonics method breaks down while the other method used by Marshak requires 
extremely laborious numerical work, but if use has to be made of expansions it is 
of course preferable to have an expansion of the correct function rather than that of 
an approximate one as obtained by Marshak. 

In the present paper we shall limit ourselves to the determination of the 
first few terms of these expansions for the spherical and the cylindrical cases. We 
shall consider the cases of large and small black spheres, and of an infinitely long 
black cylinder of large radius in Parts I, II and III of this paper respectively. 
The case of a black cylinder infinitely long and of small radius has been considered 
by Kushneriuk (1951) following essentially the same method as in Parts I to III of 
the present paper. For the sake of completeness we include, in the beginning of 
Part IV, a short summary of Kushneriuk’s article. 

In all cases we shall denote by 7 the radial coordinate and by a the radius of 
the black body, both expressed of course in units of the mean free path of the 
external medium. With these notations the expression (4) becomes 


A=,,(@))/ 1g Gd) 5 tx erein ee a ee (4’) 
For a large sphere and a large cylinder we shall determine A to the accuracy of 
O(a), i.e. we include terms of the order of a~*, but omit terms of the order of 
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a‘*log”a. For a small sphere we shall determine A to the accuracy of O(a") 
and for a small cylinder to the accuracy of O(a). The method of calculation used 
is essentially that of successive approximations, but the detailed arrangement of 
the successive approximations is different in each case. However in all cases, in 
order to avoid working out the terms which finally give no contribution within the 
accuracy intended, we start by examining the order of magnitude of various 
expressions arising in the application of the successive approximations method. 
Due to the limitations of space we shall give only a short account of some of the 
stages of the calculation while for the complete details the reader is referred to the 
original declassified reports. * 

In the second half of Part IV we give the complete curves of A(a) for the 
spherical and the cylindrical cases, obtained by combining the results of the 
calculations given in Parts I to IV and those for the intermediate values of a 
by means of the spherical harmonics method. 

Throughout this paper we use the notations 

E,(x) = | een dt: E(x)=E(x) 9 aseee (5) 
1 
Euler’s constant y=0-57722 and 6(x) is Dirac’s 6-function. ‘The symbols 
% 7, and f,, will have different meanings in each of the four parts of this paper. 


Levey ed I 


CASE OF A LARGE BLACK SPHERE 
§1.-METHOD OF THE SOLUTION 
In this case, on integration over angles the equation (1) reduces to 


rn(r) = ; ir r'n(r’) LE(|7 —7'|) — E{(r? — a?)1? + (72? —@) 88] dr’. woe CEE) 
If we put | TEAtX,T =A), THN) =G xX) n(@) os... es (1.2) 
it becomes 
1 ico) 1 r 00 
U*) = 5 |, q(y) E (|x — yl) dy — lk q(9) B {(a? + Zax)? + (y° + 2ay)"} dy, 
Sofas (1.3) 


with q(x) normalized so that q’( 00) =1, while the expression for m,,(7) can obviously 
be put into the form n,,(r) =r—"[r — a+ {q(x) — x}],—.0] 2(%) (notice that according 
to (2) and (3) q(«) behaves for large x as x), and hence from (4) 


{a(*) = X}|e=00 
a ay 1.4 
1a ql) — Jl BY. 
To solve the equation (1.3), put 9(x)=q)(«)+q,(x)+..... in which q(x) and 


4(x) (n > 1) are the solutions of the equations 
1 700 
gol2)= 5) gor) B(x—y) dy 
0 


i! oO 
anl)= 5) al) Bes) dy +fals) (n>) 
satisfying the conditions at 00, gy(00)=1 g,/(%)=0 (wm21). ....... (1.6) 
* Davison (1944, 1945); Seidel, Davison and Kushneriuk (1946). During preparation of 


this paper for publication a minor algebraic mistake has been discovered in Davison (1945). This 
is now corrected. 
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ise SS 5 | E 0.) FG Daye (9? a Day ae (1.7) 


The function q(x) is, of course, the solution of the Milne problem, and it has 
been already extensively studied. In particular it is known that 


x<qo(x)<w«+0-7104...., and that lim {qo(«) — x} =0-7104... * 


§2. LEMMA 
In determining other q,(x) it is convenient to introduce an auxiliary function | 
q(x, x) defined as the solution of the equation 


ols 0) = 5 [a0 %) Bey) d+ 8(09= 8) eee, (1.8) 


satisfying the condition at 00: E q(x, x) | ==) 4) a( Coptihite) quae ree (1.9) 


=) 


Then q,(«) for n >1 will be given by 
PAC | _ Fale) (0 8) te. (1.10) 


The function q(x,% 9) can be interpreted, except at x=x , as the neutron 
density in a half-space at the distance x from the boundary, due to a plane source of 
unit strength situated at a distance xy inside the medium, satisfying the condition 
that for finite x) the neutron current at infinity is zero. Then according to the 
optical reciprocity theorem q(x, x)) should be symmetrical in x and x). On the 
other hand if we make x, tend to infinity g(x, x9) will obviously in the limit go over 
into the solution of the Milne problem, i.e. q(x, 00) is proportional to q(x). But 
while g(x, 00) is normalized to unit current qo(x) is normalized to a value equal to 
4 at infinity. Collecting these results we see that g( 00, x) =q(x, 00) =3 qo(x) and 
hence the asymptotic behaviour of ¢,,(x) for large x, which alone will be needed in 
the final answer (1.4), according to (1.10) is given by 


In 0) =3 | i fi(8) qo &):dx= ner (1.11) 


§3. ORDEROF MAGNITUDE OF SUCCESSIVE CORRECTIONS 
Let us estimate now the order of magnitude of various q,,(x), f,(«) and q,,( 00). 
Let us assume that \g;,-4(x)| < Abe tee BS eer) See eee (1.12) 


(This is known to be correct for qo(x) with A) =0-7104...and By)=1.) Then 
one can easily derive from (1.7) 


1 
lfn(*)| S a fs Es (x? a 2ax)*? 
ees Beal. (ws Zan) ee z : 1 
+ 7p Bn 5 (Zanes pina t pent FEAR) OF othe CLokS) 
On the other hand, it is known that by increasing the free term in the inhomoge- 


neous equation (1.5) we shall necessarily increase the bounded solution (Hopf 
1934, p. 36). But if the free term is equal to [a A,_,+3a?B,_,]$£,(x) the 


* Hopf, (1934) pp. 28 ff. and pp. 47 ff. Notice that our q(x) corresponds to Hopf’s x+-q(x). 
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bounded solution is known to be equal to [a1A,_,+3a? B,_4] {qo(x) —%}, in 
which the second factor cannot exceed 0-:7104..... Thus 
lg,(x)| <0-7104... [a A, + 3a B,,_1], 
i.e. (1.12) is valid for any n provided we take for n > 1 
Boe (OglUe mde je O'/l0ttoa)s BO ssa. (1.14) 


This shows that our process of successive approximations converges for all 
a greater than 0-7104.... Using (1.13) one can show further that, if A’ and B’ are 


some positive numbers, 


an ee AGA 3B. A 4A. BO 9B eB 
i (A'+9B’)f,(x) de] < XS 4 St ty 


Combining this with (1.11) and (1.14) we see at once that 
Gal Coy Ol ae Rt). a frest Selmer 205.53 (1.15) 


In a similar manner one can also show that if in working out q,_,(x) we 
neglect a term of the order of C(a) x(log x)" the resulting error in evaluating 
(x) is of the order of C(a) a~”(log a)* while the error in evaluating g,,(00) will be 
of the order of C(a) a~”—*(log a)". 

Comparison of the order of magnitude of q,,(00) and of q,,(x) for unrestricted 
x, suggests that the part of q,(«) which is of O(a"), n>1, should come entirely 
from the free term of the equation (1.5) while the iterated term should be of a 
higher order of smallness. In fact one can easily show (see Appendix A) that 


| @n(*) —fn(x) | =O(a-""1 log a) +a-" O(w log x) (mn >1). ...... (1.16) 


Then, if in working out q,,,;(0) we replace, in the definition of f,,,,(x), ¢,,(y) by 
T,Ay), the resulting error, according to the previous remark, will be of the order of 


a" *loga. 
Estimating in this manner we obtain in particular 
{9(%) — ©} | aco = {go(*) — ¥} leo + G1( 0) + G2( 00) + O(a *loga) ...... (1.17) 


in which g,(00) is the expression obtained from the expression for q.(00) by 
replacing qi(y) by fi(y). 


§4. EVALUATION OF INTEGRALS INVOLVED IN SECOND 
AND THIRD APPROXIMATIONS 


Now, combining (1.7) and (1.11) we have 


an(c0)=— 5 [ - gals) gay) E{(a? + 2a)! + (9* + 2ay)!9} de dy 


Gu(co) =F [|| gala) gals) BA(o? + 2ax)!*-+ (32 + 2ay)!®}EA(92 + 2ay)# 
: + (22+ 2az)'*}dxdydz. ...... (1.18) 


To evaluate the integrals (1.18) for large a we need the expansion of g)(«) for small 
x. The derivation of this expansion from Hopf’s solution for go(x) is straight- 
forward though laborious. It has been carried out by Mark (1944), who obtained 
the following result : 


Q(x) = 5 {1—43xlogx+c,x+O(x?log?x)}, ...... (1.19) 
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in which c, = 1:278794. Substituting this into the expression for q,( 00) we get 
gi( 00) = — An | 2 E{(x? + Zax)? + (y? + Zay)1?} dx dy 
0 Jo 


- [- IL «(—4log x +¢,) E{(x? + Zax)? + (y? + 2ay)1?} dx dy + O(a~* log?a). 
OR-n0 
fo. ee 


Making the substitutions s =(x? + 2ax)12, t =(y? + 2ay)!? we can easily see that — 
the first term in (1.20) becomes 


sds A parted eu tee yi? i; 
-3)° ( E(s +t) —=—su5 (@+a)2 era al iF st E(s +t) ds dt + O(a) 
Al 
at Olas) fie See (1.21) 


Similarly, the second term in (1.20) gives —(1/4a*) (log 2a+2c,+ 2y—10/3)+ 
O(a log a), while for g,( ©) we obtain the expression 


Go( 00) = ial, bs (. stu E(s +u) E(u+t) ds du dt + O(a log a) 
JoJo 
= As (log 2 — 2) + O(a log a). 


§ 5-RESULTS 
Sean the above into (1.17) we obtain 


1 
{q(x) —#}| po =0°7104... — Gite EO DID (cy +y— 3) + O(e4 log a), 


while substituting the last expression into (1.4), expanding in powers of a~} and 
using the numerical values of c, and y quoted above, we finally obtain 


A=0-7104... +0-5047 a4. 0:2336 a-2— da-3 log a—0-1704 a8 — O(a-4 log?a). 


Seas (1.22) 
It should be pointed out that while the series g(x) = qo(«) + q,(*) + ... 1s, as we 
have seen, absolutely convergent for all a greater than 0-7104..., the series 


obtained by expanding (1-22) will, presumably, be only asymptotically convergent 
since, for instance, the expansion of 


| | E{(x? + Zax)? + (y2 + 2ay)"?} dx dy 
Jo Jo 


in powers of a~1 (see (1.21)) is only asymptotically convergent. 


PA Rela 
CASE OF SMALL BLACK SPHERE 


§1. GENERAL OUTLINE OF THE METHOD AND 
EXAMINATION OF CONVERGENCE 


We start, as in the preceding case, with the equation (1.1). However, instead 
Ese Lp n(p)=1-In(co) 4G) a (2.1) 
so that the equation (1.1) gives 


a(r)=3{ a(t) {Er—t) EG?) 4 (8 a2] dtte(r,a) (ra) 
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a ahieh e,(7, a) = 3{E,(r — a) — E,(r? — a®)¥2-—aE,(r—a)} (ra) ...... (233) 


and g(r) is that solution of equation (2.2) which is bounded at infinity. Notice 
that our present (7) and «,(7, a) are identical, except for the sign, with Marshak’s 
QO(r) and ¢,(r) (Marshak 1947 a). 
n,.(7) will obviously be given by ”,,(7) =r 1n(00) [7 + q(<0)] and hence, from (4), 
Nei (Goo) = a eattipe Wit Ni) Fe Weta (2.4) 
Turn now to the solution of equation (2.2) for small a. This can be done by 
means of the perturbation method. In order to apply this method, it will be 


convenient to extend the definition of g(r) and «,(7, a), for the region 0 <r<a. 
This we shall do by putting there 


g(r) = 3|. dt q(t) {E(|r—t]) -E(r+2)} (O74), 


ae ae (25) 
e,(7, a) =0 (0 <r<a). 
For brevity we shall write 
AB) = 5] dep(e){B(lr—tl) Ee +9), 
Ap)= 5) deg{Blr—tl) Ee +), 
Avg, 2 2)1/2 2 2\1/2 
See all dt $(t)[E{(? —a2)12 4 (2 —a2)12}—E(r+2)] for r>a, | 
0 for” 0 <7r<a. | 
fee (2.6) 
Then the equation (2.2) can be re-written in the form 
g(r) =Aq(r) — (Ay +As) g(r) +es(7,@). sn ee (2.7) 
Put now Qe Ged 1) seta 44 VE aniecess (2.8) 
in which q,(7) is the solution, bounded for r— 00, of the equation 
(1—A)q,(7) =f,(7) (Hee elocate cee) | saint. ccos (2.9) 
in which lt ean Lp biol) et Nat Aa) ot), he L). ane ee (2.10) 
One can easily establish that, for 7 >a, 
€,(7, 4) =O(a*e™”); Ayd(7) = Max |4(7)| O(a e~”), | (2.11) 
A2$(r) =Max|4(7)| Oe’), 


while for r comparable with a, or r <a, 


€,(7, 4) = O(a); Ay(7) =Max|4()| O(a), Asd(r) =Max | (7) O(a), ... (2.11) 


in which ¢(r) is any function for which the operators A, and A, exist. 

Now, since the kernel of the equations (2.9) is independent of a, Max]q,,(7)| 
should be obviously of the same order of magnitude as Max|f,(r)|. ‘Then 
according to (2.10) and (2.11’) q)(r) =O(a), 9¢,(r) =O(a?) and so on. Hence it 
follows that the series (2.8) will, for sufficiently small a, converge and its sum will 
be the solution of (2.7). 
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The solution of each of the equations (2.9) can be obtained, for instance, in the 
following way: [q,,(7) —f,(r)] can be regarded as 7 times the neutron density in an 
infinite (isotropically) scattering medium (occupying the entire space, i.e. in the 
absence of any black bodies) due to a spherically symmetrical system of distributed 
sources, the source density at the distance 7 being f,(7)/r. Then if p(r) is the 
neutron density due to a point source of unit strength in an infinite (isotropically) 
scattering medium in the absence of black bodies, ¢,,(7) will obviously be given by 


alr) = ||] Z2r—rol) fale) rot halt) noe (2.12) 
and p(r) entering into (2.12) is already known (Bothe 1942). 


§2. SIMPLIFICATION OF CALCULATIONS 
The problem of determining A to any prescribed accuracy for small a is thus 
reduced to the evaluation of certain multiple integrals. We shall consider now 
how to minimize the complexity of these integrals. If we wish to determine A to 
the accuracy of O(a?) then, in view of (2.4), and since g(00), as we shall see later, is 
O(a’), q(©) should be determined to the accuracy of O(a’*). 
It has been shown by Marshak* that 


3a? ee 
q(o)=- = +3] g(r)ex(r,a)dr 
a 
combining this with (2.7) we have 


(oo) =— =F +3 | “ Lesra)Pdr—3 | es(rsa)(Ar+As)a(r) dr 


43 | a(t, a) Aad (2.13) 


It follows from the estimates (2.11) and (2.11’) that the second term in (2.13) 
is of O(a*) and the third term is of O(a*). Let us now show that, provided we 
refer to the quantities of O(a” log* a) simply as O(a”), the fourth term in (2.13) 
is also of the order of a. 

We estimate first Af,,(7). Using (2.10), (2.11) and (2.11’) and understanding 
O(a’) in the sense explained above we can easily derive 


NF, (2) =O Ga) (for’any, value Of 7 )> ae (2.14) 
On the other hand applying the operator A to (2.9) we obtain (1 — A) Aq,(r) =Af,(7), 


and hence by the argument already used before Max | Aq,(r) | should be of the same 
order of magnitude as Max|Af,(r)|. Using (2.14) then 
NGG) —O(e 7?) (for any valuéof 7) (2.15) 

while, combining this with the estimate of the order of magnitude of ¢,(7, a), we 
see, in fact, that that last term in (2.13) is of the order a*. If we are working out 
g(©) to the accuracy of O(a*) this means that we shall need only the leading 
terms in the third and the fourth integrals in (2.13), and consequently we can 
replace g(r) there by q)(7). Further it follows from (2.15) that (A, + A,) Ago(7) is of 


a higher order of smallness than (A, + A,)e,(7, a) and because g)(7) =Aqo(r) + €3(7, @) 
we can replace the third term by 


= i é,(7,-4) (A, Ag) e,(z, a) deat Ola) eee (2.16) 


* See Marshak (1947 a) formula (23) bearing in mind that our g(r) and «,(r, a) differ from his 
QO(r) and «;(r) by the sign. 
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Let us now consider the last term. We will first approximate the expression 
for qo(7), so as to avoid the evaluation of the rather unwieldy integral (2.12). Let 
us start by obtaining an approximate expression for f,(7), i.e. for €3(7, a). 

For r>a, expanding it into the powers of a and using the recurrence 
relationships between the functions F’,(x) we obtain 


«,(7,a) = — oe {1+0(@+0(5 5) brea, tee (2.17) 


For r comparable with a the approximation (2.17) is no longer valid, the error 
becoming comparable with «,(7, a). However, the approximation (2.17) can be 
used to evaluate the leading term in AAe,(7,a) for any r. This can be seen as 
follows: According to the definitions (2.6) 


Al o,+ | =5[° [at a+ | (r-21)- E(r+t)}dt. ...(2.18) 


If r>a, then, for ¢ comparable with a, {E(|r—t|)— E(r+1t)} is of the order f. 
Now {te,(t, a) +4a*e‘} = O(a") and the contribution from this range of ¢ to the 
integral (2.18) is of the order a*, while for t> a the integrand (according to (2.17)) 
is of the order a?. Thus for r>a the integral (2.18) is of the order a’. 

Similarly one can show that for r comparable with a, the integral (2.18) is of 
the order a?. ‘Then, using a derivation similar to (2.14), one can show that for 
any value of r 


AA Ee a)+ | =2¢)9 eae (2.19) 


Remembering now that p(r), which occurs in (2.12), satisfies the equation 
(1 —A) p(r) =e""/4nr, we see that (1—A)AA[q)(7) +7 a@*p(r)] = O(a?) and hence, 
similarly to (2.15) AA [q)(r) + 7a? p(r)] = O(a’), so that the last term in (2.13), on 
omitting the quantities of O(a’), can be re-written as 


3 | ie ee a) de ora! | “PACA AVi we ries Cao e (2.20) 


The last term in (2.20) can be transformed further. Substituting for A from 
(2.6) and changing the order of integration we have 


—3na* AS) p(r) AAe,(r, a) dr 
0 


which in view of (2.19), and omitting terms of O(a°) becomes 


Tr a! ite pir) AA — dr. Re... (2.21) 


Collecting the above results, and abandoning our convention of writing 
O(a”) instead of O(a™ log* a) we obtain 


goo)=— +3] fes(rsa)PPdr+3 |” s( 4(r, a) (A—Ay—Ag) e5(1, a) dr 


+ 2 a | pir) AA — dr +O(a@log*a). ...... (2.22) 
0 
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§3. EVALUATION OF THE INTEGRAL (2.21) 

Now f(r) is known only in terms of its Fourier transform: it has to be determined 
by the inversion formula, which would involve numerical integration. However, 
in order to evaluate the integral (2.21) we do not have to know p(r) explicitly ; it is 
sufficient to know its Fourier transform, if we make use of Plancherel’s theorem 
(Titchmarsh 1937, p. 72). According to this theorem, if 


AQ) = J? | F a(r) sinur dr and B,(u) = J2) b(r) sin ur dr 
710 WAY 
are the Fourier sine transforms of a(r) and 6(r), then 
| a(r)b(r) dr= | A,(u) B,(u) du. 
0 J0 


Now the Fourier sine transform of e~"/r is (2/77)! tan~! u, while the Fourier sine 
transform of Ad(r) is equal to u-1tan-1u times the transform of ¢(7). Then from 
the equation (1—A) p(r) =e-"/477, it can be seen that the Fourier sine transform 
of p(r) is utan-4u/(27)??(u—tan-1u), while the Fourier sine transform of 
AAe*/r is u(2/7)1? (u-+tan-1u)®. Therefore integral (2.21) is equal to 
3a4 ip (tan-1u)* du 
Sr Jo u(u—tan— lu)” 


On making the substitution tan-1u =@, it becomes 


Sat ae 64 dé 

‘Sir if sin@[sin@—@cos@]) Cz) 
This last integral can be easily evaluated numerically, and the value of the 
expression (2.23) was found to be (7233704 ee ee eee (2.23’) 


§4. EVALUATION OF THE OTHER INTEGRALS ENTERING (2.22) 

To simplify the evaluation of the other integrals entering (2.22) it will be con- 
venient to replace «,(7, a) by some approximation in terms of more elementary 
functions giving the desired accuracy for small values of a. To obtain such an 
approximation re-write the expression for ¢,(7, a) in the form 


—r 


(7? —a?)12 
€,(7, 2) = 3 | exp a—exp {r—(7r? —a?)"?} —(7? — a?) | = exp (r— | 


and expand exp a, exp{r—(r? — a”)""} and exp(r—f) into powers of a, {ry — (7? — a?) 1} 
and (r—t) respectively, retaining the appropriate number of terms. Notice 
further that if we wish to evaluate a 
e4 brayar, 

a 

say to an accuracy of O(a*tlog* a), and if F(r, a), for r>a, is of the order of a® we can, 
in this integral expand e~?” into the powers of r and retain only the few first terms. 
In this manner one can show that 


Res r@ fate = a®[3r + (7? — a?)1?](1 — 27) 
pe 
ie [e3(7, a)| d= x { 16r2 1672 [r me (7? 7. @) i] 
a r—a@a  rta at 
~ Afr + (P— a1} Ee log ag A. fr+(r— aoe dr + O(a? log*a). 
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The integrals entering into the last expression, apart from the one leading 


to E-functions are either elementary or expressible in terms of the series 
ss ; 
2X ,=1 1/n? =7n?/6 , and the ultimate result is* 


(ea) 3 
3 | [ex(r,a)Pdr= + Sa! | ogat+y + 5 | + O(a log*a). ...(2.25’) 


ones 
16 
Turn now to the integral — | fs €,(7, a) (A, +A.) €3(7,a)dr. Notice that, 
according to the definitions (2.5) and (2.6), A,e,(7,a)=0. Inevaluating 
| ex(rs@) A es(r,a) dr 


we can replace e,(7,a) by — a?/{r + (7? — a”)"?} + O(a?), on the grounds of (2.24), and 
similarly replace the function [E{r? — a?)"? + (¢? — a?)1?} — E(r + t)] occurring in the 
definition of the operator A, by the leading term of its expansion for small 7 and 
t, i.e. by log (7 +2)/{(7? — a?)1? + (#2 — a?)"?}. Thus, we have 

ay { €4(7, 4) (A, + Ae) «g(r, a) dr 


dr dt r+t 


a Se ee ee or 
~~ 8 Ja Ja HP eR i+ (P—aye 8 aes (Pays FO) 
Ie ees (2.26) 
on evaluating the integral by means of series we get 
Sills 5 37? 2 é 
Bara = (2n+1)8 cy a} GiOlG?)> wpeecieesiest (2.26) 


Consider now the last remaining integral, ie. 3 | €3(7, @) A €,(7, @) dr. 
On the grounds of the considerations developed above it can be re-written as 
Lay San fae Cee ard : 
4), |, eRe aE aay A Ble +} + OFe8 
or alternatively as 


3 ; 0 o-T © pt 
=a I. a —[E(t—r) — E(t-+r)]at dr 


Se ae 1 1 ) t+r 3 
50 |e ae eae em — np ei, +O 


The second term in (2.27) can be evaluated essentially by the same means as 


(2.26) and the result is 
3 52? I 1 Cy al 
fy) | See, RE wr =i,-8 eee 3 2 
a {3 128 Se gm log 2 76 18 2+ ea ra SM Celeeue ete (2.28) 


To evaluate the first term in (2.27) notice that for small r 
oO pt 2 
| —[E(t-7)- E(t +r) dt= = + O(rlogn), 
0 


* Tn the original report the arrangement of calculations differed in one respect: the term involving 
[4(r?—a®) log {(r+-a)/(r—a)}—ra] was excluded from the formula corresponding to (2.25) and 


was considered in connection with fesAeg dr. 
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and hence this first term becomes 
Soe rem cea oc . = . 
35 4 E(a) + %a* | “1 =a [E(t—r)— E(t+17)] dt— = hdr O(a®). 2 .( 22a) 
0 a 

The integral appearing in (2.29) after some rearrangement can be reduced to a 
form which can be evaluated by means of a series, and the expression (2.29) — 
becomes 32 sam 1 
32 a’ E(a) —54a = (2n+1)° i . 

Ss aRhS Ul nS 

After substituting into (2.22) the expressions (2.23), (2.25), (2.26’). (2.28) 

and (2.29’), we ultimately obtain 


3 a 38/a he 
g(%o) = — Z q* + ee Ver — 1) atlog a—0-87095 + O(a? log*a), ...... (2.30) 
while substituting this into (2.4) we have 
A =$—2a—0-97827 a? log a—1-4002 a? + O(a? log* a). ...... (2.31) 


PA Kel aetsl a) 
CASE OF A LARGE BLACK CYLINDER 
§1. TRANSFORMATION OF THE EQUATION 


To reduce the equation (1) to a form suitable for the present case we should | 


first of all perform the integration over the surfaces r’=const. occurring in the — 
right-hand side of that equation. It will be convenient to take for the three co- 
ordinates describing the position of the point 7’, 7’, z’ and p where 2’ is measured 
along the axis of the cylinder, 7’ is the distance of the point 7’ from the axis of the 
cylinder, and p is the distance between the projections of the points 7 and 7’ on 
plane z=0. Then the equation (1) becomes 


rae rata) | 0 ee 
7 = P= 0—r ear p pe, MOM 
si the ae (3.1) 
where K,(¢) is the modified Bessel function of the second kind. If we expand 
{(7 +1’)? — p?}-1? in powers of {p? —(r—7')?}/4r7’, and put 
r=4405. 7 =ays. 72n(r) = Clg), a aeeeee (3.2) 
(in which C’ is some constant) it can be re-written as 
Ii 2 Po(O 
q(x) = 9} | F q(y) dy an nat ie Aa 
x{F,(|x—y])—F,[|x—y]5 (2ax+02)2+(2ayty?2]}, 2... (3.3) 


in which P,,,(s) are the Legendre polynomials and F,, and F,, are defined by 


~ De (ee ee) 
Fy(n;k)=2 |” (pty |” Kygdedp; (A>h>0) | 


TT J is } 
F,(h) =F,,(h, h). 


The functions F’,(h) can be shown to be connected with the E-functions (defined 
in (5)) by the relationships 


Fy(h) =H); F(A)=Qn-1)[- tF, Od. .....8.5) 
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§2. NORMALIZATION OF THE SOLUTION AND THE 
ARRANGEMENT OF THE SUCCESSIVE APPROXIMATION 
Before proceeding to solve equation (3.3) let us examine the behaviour of 
q(x) for x> a, remembering that in our present case a in itself is very large (a> 1). 
Using the formulae (2) and (3) of the introduction, noticing that the general 
solution of the equation ( Fee ae) 


aaa e 5) n,(7) =9, 


is 7,,(7)=A+Blogr, and expressing here A and B in terms of the current at 
infinity and the linear extrapolation length A we get at once 
nr) =((7/a) 2'(r)\|,n0 {a@logrjat+A+Ofe*™®)} i... se (3.6) 
Then if we write the constant C’ in (3.2) in the form C’ =(1/C v/a) [rn'(r)]|,-. 
in which C is another constant to be determined later, the expression for q(x) 
becomes 
Qa)=Clal aja) log (1-Fx/a) + NI 4/2) +O *)). 222. (3:7) 
In this expression the only quantities that are not already known apart from the 
common factor C are the constant A and the terms of the order of e~*. Both could 
be determined, for a>1, by examining the variation of g(x) in the region x <a. 
Now the kernel of the equation (3.3) for large |y — x| is essentially of the order of 
exp(— |x —y|) so that if x is small compared with a, the contribution of the region of 
y either comparable with a or large compared with a, to the value of the right-hand 
side of (3.3) should be essentially of the order of e-. While the contribution of 
the region of y <a to the values of g(x) for x <a will consist of terms of the order 
of unity, of the order of a! and so on. 
This shows that, in so far as the determination of q(x) for «<a is concerned, we 
can replace (3.3) by the equation 


q(x) = ; ia q(y) E(|x — 9) dy 


Gre P, (0 - 
ap a)! a Fy( (xy) — Fully] Qax+ 2)! + Cay-+y°)!91} dy 


1/” | P20) 
a5 ay) - (x+y) =e Ti a=y))o dy Bae 2 ea (3.8) 
We solve this equation by the perturbation method, i.e. by putting 
UX) = o(%) +. ga(%) + 99(4) + 2.52. ars. (339) 
in which q(x), g2(x), q(x) etc. satisfy the equations 
pee 
gol®= 5) WO) Ey ey ee (3.10.1) 


qo(X) = an go(y) E(|x —y]) dy + an qo(¥) | a | Fey) Alley: 
(Zax + x?) + (Zay +") "i} dy POSHGCRY (3.10.2) 
q3(*) = An qa(y) E(|x— yl) dy — AN qo(¥) (x+y) ead F,(|x—y|) dy 


Tree Me 
Z 31, qx) Fol lx —y]; (Zax + 2°) + (Zayt+y*)?]dy, os... (3.10.3) 


and so on. 
In the equations (3.8) and (3.10) we have extended the range of integration up 
to y = 00, rather than integrating only over the range y <a, since this change in the 
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range of integration can affect the result for x <a only by a quantity of the order of 
e~*, as in the original equation (3.3). 
o(x) is the solution of the Milne problem, and is known to be asymptotically | 
linear in x. We have a normalizing condition, achieved by making q9'(«) =1. 
In so far as the other q,,(x) is concerned, obviously if ¢,,*(x) is a possible solution of 
the equation for q,,(x), so is {q,,*(x) + Ago(x)}, where A is anarbitrary constant. We 
can, however, eliminate this arbitrariness by stipulating that the asymptotic 
expression for q,,(x) for «> 1 should not contain a term proportional to x. Since 
the normalization of g9(«) was chosen in an arbitrary manner, the introduction of 
the above condition on q,,(x) for n>0 obviously will have no effect on the shape of 
q(%) = U0 In(*)- 
The reason for combining the terms 
(|P.(0)|/4a") Fy(|v—y]) and = — Fo[|x—y], (Zax + x?) P+ (2ay +9") ] 
in the expansion of the kernel of the equation (3.8) is as follows: With our 
restriction on the solutions of inhomogeneous equations of the type (3.10), the 
solution of the equation 
qoi(X) = ; | qe(y) E(|x— yl) dy + : [ qo) ee F,(|x—y)) dy, 
J0 J0 


for «<a is obviously of the order of a. On the other hand, as is shown in 
Appendix B, the solution satisfying our restriction of the equation ; 


1/2 
Go(*)= =|  ao(y) E(|x—y|) dy 
2J0 


is ~ 
os ae go(¥) Follx —y|; (2ax+«2)12 + (2ay+y?)2] dy ...... (3.11) 


is, for x >a“, also of the order of a*. Thus q2(x) = q21(*) + qoo(x) is of the order of 
a for a1<x <a, and this is just the range of values of x for which we need to 
know g(x) in order to determine A from (3.7). 

It will be also shown in Appendix B that the appropriate solution of the equation i 


1 - 00 
qo(8) =| aul) B(x —yl) dy 
17% : 
25 (20) Polley]; 2ax+ x)? + (2ay+y)P]dy oo... (3.12) 


is, for x>a~, of the order of a~*. So that q(x) is for x in the range a-!<x <a of 
the order of a? and so on. 


§3. EXPRESSION FOR A IN TERMS OF 
ASYMPTOTIC PARTS OF q(x) 

Considering the equations (3.10) for nA0, one can easily show by induction 
that the free term of the equation for q,,(x) is a polynomial in x of the order (n —1) 
plus an exponentially decaying term O(e~*); and hence q,(x) in itself should be 
equal to a polynomial in x of the order ( + 1) plus an exponentially decaying term 
O(e=*). 

Remembering also that q¢,(x) for nA0 is assumed not to contain a term 
proportional to x, we can thus write 


an(*)={2? Onat+#Onet +210, .34+G,(8) (e2) .-.- (3.13) 
in which q,,(x) is a bounded function and Q,, ,,, are some constants ; while g)(x) can 
be written as qo() =x + qo(«) in which gp is also bounded. Substituting these into 
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(3.9) and noticing that for x tending to infinity 7,,(x) tend to constant values we get 


Qx)= 2% g,(co)+x+x* 2D Ono+..3 + O(e-*) 
n=0 n=2 
while (3.7) can be re-written as 
r Cr MS 
q(x) =CA+xC 1+5, tat gate 4 e”); 


Comparing the above two formulae we see at once that we should have 
C(1+A/2a)=1;CA= &%* 9 q,,(«) and hence eliminating C and solving with respect 


to A we get os 
Ee) (oa) 1 (oe) m 3 14 
A= aS = » { >) n (ee) } © 6:6 16-0 (0 
ae ee m=0 (2a)™ sry ant ) oY 
2a nap” 


$4. SECOND AND THIRD APPROXIMATIONS: DERIVATION OF THE 
EXPRESSIONS FOR 3,(©) AND 3,(0) 

It is evident from (3.14) that in order to determine A to the accuracy of O(a~*) 
we should determine &7_)q,,(00) to the same accuracy. While since q,(x) is of 
the order of a~” for a!<x <a, it follows from the definition of ¢,,(«) that it is of 
O(a") for all x>a~1; and in particular ¢,,(00) is of O(a-"). Also qo( «) is already 
known and is equal to 0-7104...(Hopf 1934, p. 85). Thus, all that we need to . 
determine is ¢,(00) and q3(00) to the accuracy of O(a-°). 

The polynomial parts of the free terms of the equations determining qo() 
and q(x) can be easily worked out and are {x+q)(o)}/12a? and 
— {6+5x[x + Go(20)]}/30a* respectively. Then one can easily verify that the 
formulae (3. is forn =2 andn =3 become 


1 - 3 
q2(x) = Aa 2 ee a 3x? Go() | a qo(X), q3() i 2443 [x4 <7 2a Jo()] a 93(X). 


Substituting these into (3.10), and expressing all the integrals of the form 


| EREDAR GE 


in terms of E,, functions we get for g,(x) and J,(x) the equations 
q2(x) = An g(y) E(|* — yl) dy — ; : gol) Fo dy 
+ apga le? 3B) — dolce) fe-* 3B) 
aS . [do(¥) —¥ — Go(%0)] [Ea(|x —y]) +|*—y|Ea(|x—yl)] dy} 
and Gy(s)= 3 [ H)EUe—yddy— 5 ar) Foay 
+ spcp 3 —m)e*— 15 Bylx)] Jol 20)l(2— ae*— 9E4(3)] 


. ie [qo(¥) —¥ — Gol 20) (~ + y)[Ea(|x—y|) + |x —y |E2(|x —y])] dy}, 
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while using the formula (1.11) in Part I of this paper, i.e. 
7.0) =3 | o(x){the free term in the equation for ¢,,(x)} dx, 
0 


we finally get (00) ={Ly+Lo(a)}/a2; Gq( 00) ={L3+ La(a)}/a?, ...... (3.17) 
in which 


0) 


L,(a) = — 5 a? ee “| qolx)aoly) Folly: (Zax + x?)¥° + (2ay +y*)¥"] dx dy, 
3 
2° | 


Efa)=— 50 | | gls)aaly)Follx—31s Cae +ae)!*-+ Cay +9%)**] de dy, 
0 


while the expressions for L, and L, are self-evident from the comparison of 
formulae (3.16), (3.17) and (3.18). 

The expression for L,(a) depends upon q,(y). However, as it will be shown in 
Appendix B, of the two parts composing qo, viz. of go, and qo, the contribution of 
4o1(y) to L,(a) will be of O(a~4). Also if we substitute for q.9(y) its expression from 
(3.11), only the free term of (3.11) will give a contribution of the order of a~* while 
the contribution of the iterated term will be of O(a~*). Thus in our present 
approximation we can replace the second formula (3.18) by 


~ 3 pO rao ow ~ ie 
Lfa)= 5a) | | aola)aole)Falle—yls 2av+se)!" + (2ay +9*)!* Folly — 5 


2ay + y?)12 4 (2az+2°)2]dxdydz+O(a*). — ...... 3.18% 
yay 


§5. EVALUATION OF THE INTEGRALS APPEARING IN THE 
EXPRESSIONS FOR @,(«) AND 4,(«) 


Consider first L, and L;. One can easily verify the identities* 


an qo(9)Ex(|x —y]) sign (y— x) dy =}, 
5), aor) Ealle—yl) dy = 4 +G9(0)}, 
I ra eg he ei - (3.19) 
5 an qo(¥) E.(y) dy = = sel, [o(*) — * — gol c0)] dx, | | 
‘1G Lee ic ms 
31, a0)EsO) y= 3o(ce)— 5] slau(s)—2— Fue) dy. | 


* These and other similar identities can be derived as follows, put : 


S,0)= 5 [alo B, (ley) sign Vay" ay, 
Then obviously Si, 4 (x)=S_, (x) —(2n)~ "qo(x)i1 —(—1)"1, 
while for x — «, from the asymptotic expression for g)(x) we get 
S,(x) ~{x+G(~)}/n if n is odd, Si,(x%) ~1/(n+1) if n is even. 


Combining these with (3.10) one can easily determine S;,,(x) by induction, and thus verify (3.19). 
Some of the identities (3.19) have been given by Hopf (1934, p. 30). 
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Using these identities and the recurrence formula for the F,, functions, one 
can eliminate all the double integrals in the expressions for L, and Ls, as well as all 
the integrals involving F,, functions, obtaining 


one iz ee . ge ie 
Ly = 40 qo( 20) + 3 dol 00) | ‘ [4o(*) — * — go( 20) ] dx — all - aL o(%) — # — Gol 20)] 
L. = a lg 2 lL. Hes [ = ld ' 
2 = Fp lGo(©)] ey x ga(*) —¥ — Go( )] a | 
The integrals appearing in (3.20) have been already considered by Mark (1944) 
who has found that 
iii = Noegs bx 
J ele) aco) v= 4 = Flo), | 
(0) 


and ie x dol) — 2 —G,(co)ldx = —0-025993. | 
0 


Turn now to the integrals L,(a) and L,(a). These are quite similar to the 
integrals (1.18) occurring in Part I of this paper. Following the same method as 
there, and noticing that 


F, (\(@+s°)!®—(@ +22]; s+) = = | Ko(z)log — dz +0), 
Y st+t i 


we obtain after a slight rearrangement 


2 i Rd a ie se S*logs log?a 
NOs = le st E $5 (log 24526) — +0/( = )| G(s +1) dsdt 
ee (3.22) 
and 
Lp) = =| | le ds, ds, ds G(s, +5») G(5_ +53) +O (“2 ‘) Hate (3.22') 
in which G(s) = lg K,(2)log(z/s)ds. saa (3.22") 


An exact expression can be obtained for the integrals occurring in i (a) and 
L,(a) can - expressed in terms of rapidly converging series as follows: 
5 loga £37 
L,(a) = — a aie eae ate § 2er+2y +1og8)} +0 (8 ) mations) 


in which c, and y have the same meaning as in Part I of this paper and 


log?a 


eetl>: <= 1 Sagem cta es | 5log?2 5log2 
ts 167? 22 (2n+1)® 167, 23 as 967?  —- 192 
347 19 log a ry 
- seam — 79g +0 (“E*). 0:.0).0\0. 8.0 G23") 


6s eRES WLS 
Carrying out the summation in the series (3.23’), using the known values of 
Vy €1, Fo( ©) and (3.21), substituting (3.23) and (3.20) into (3.17) and then the latter 
into (3.14) we ultimately obtain 


A=0-7104.. +0°2524..a724+.0:0949 a2 — aa log a—0-01634 a-8 +0 fie ‘) 
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PPAR ey, 


SMALL BLACK CYLINDER AND BLACK SPHERE AND BLEACH 
CYLINDER OF INTERMEDIATE SIZE 


$12 CASE OF A SMALL BEACHK-CYLPINDER 
The solution of the problem for the case of a small black cylinder has been given 
by Kushneriuk (1951) and we shall limit ourselves here to a short abstract of that 
work, to which the reader is referred for all the details. 
The equation governing n(r) is the same as the equation (3.1) which, on putting 


Z (2 —a?)LP 4 (472 — a?) dp 

BG ae lar Poorer yaar |, Kloet 
Eye (4.1) 
can be re-writtenas n(r)= | TMT ) KG, 7 Gy di 0 eed) ae (4.2) 


Following the same method as is adopted in Part II of the present paper, the 
definition of n(r) is extended to 0 <r <a by putting 


Sut) = (3 Cnr Kr, 4 5 ONGE® SF a) a eee (4.2) 


and the solution is sought in the form 
nr) =, (7 )s-14(1) + a i ee (4.3) 


in which the successive ,,(7) satisfy the equations 
Ny (7) = | r'no(r’) K(r, 7’; 0) dr’ | 
0 


(oa | ri Ar) Kr 7: 0) drei 7) gn) 
~ 0 
in which 


ar eS) a a WKGir a) Raya [rq alr VK, r', 0) dr’ 


and S(r— a) =1 forr >a and zero for r<a; n(r) is normalized so that 
‘(jan @)\,..= 1. (Se ae (4.6) 


Consistent with this normalization one can take m)(r)=1 and normalize various 
N,»(7) for m>1 so that, for ro, n,,(r) behaves as 


n,,(71) =B,,(@) log ria O(€7)., | nee (4.7) 
Then, as one can easily see, A will be given by 


To solve the equations (4.4) for m > 1, we can use the following formula: 


inl?) = ||| PCR) fn(r0) do dyy day thy?) seve (4.9) 


(this is analogous to (2.12) of the second section) in which R is the distance between 
the points (x, y, %) and (x9, Vo, %) and p(R) is the neutron density due to a point 
source of unit strength in an infinite isotropically scattering medium, at the distance 


vo 
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R from the point source. Substituting the explicit expression for p(R), choosing 
the additive arbitrary constant entering that expression so as to satisfy (4.7), one 
can re-write (4.9) as 


tn(?) = By (a) log r/a-+ fy(a)—3 | {log (|r) f(r?) dr 


fea) ice) “270 
a ral ay) dy | 4’ flr") dr" | Kofy(o? +12 —2rr! cos 6)!}d0 (4.10) 
1 0 0 
in which © AREy== 6 ee (CRUE gis eee (4.11) 
0 
fs fil Se hfs, Ir i)} 
and a(y)= {1 Mag Beet + re (= + log a ge (4.12) 


In applying these formulae, it is necessary to determine beforehand the order of 
magnitude of various terms, so as to avoid carrying out lengthy calculations of 
quantities which would be subsequently found to be negligible. The main result 
in this direction is as follows. If, for some m, f,,(r) satisfies the inequalities 
lfin(71)|<Ame/r for r>a, |f,,(r)|<A,, for r<a, then f,,4:(7) satisfies the same 
inequalities, with A,,.,<C A,, alog a—}, in which C is some constant independent 
of m. Using this result one can show, for instance, that if we wish to determine A 
to an accuracy of O(a) i.e. neglecting terms of O(a? log? a~4), it is sufficient to evalu- 
ate B,(a) and B,(a). Working out of these two coefficients gave B,(a) = #a + O(a’), 
B,(a) =(3/7?) — (9/16) a log a+0-1217 a? + O(a? log a) and hence 


A=4/3 + (1 — 16/37?) alog a—0-2164a+ O(a? log?a). 


§2: BLACK SPHERE AND BLACK CYLINDER OF 
INTERMEDIATE SIZE 

The preceding Parts of this paper dealt only with the determination of A for 

the cases of small and large values of the radius, rather than with the case of the 
radius comparable with unity. In this latter case the exact determination of A is 
impossible, and one has to resort to approximate methods, to which reference has 
been made in the introduction. However, supplementing the results of the 
application of the approximate methods with the values of A for small a and for 
large a from Parts I to IV and from Kushneriuk’s paper, one can judge fairly 
accurately what the exact values of A are for intermediate values of a. ‘This is 
illustrated in Figures 1 and2. Figure 1 refers to the case of a black sphere. The 
curves I and II correspond to the expressions for A for small and large a derived in 
Parts I and II of the present paper. ‘The curves III(a) to III(d) represent the 
results of Mark’s calculation of X by the spherical harmonics method using P; to 
P, approximations, the application of this latter method to the case of a black sphere 
has already been described by Marshak*. From the comparison of these curves, 
the curve IV was drawn, joining smoothly curves I and II for small and large a, 
and having the same shape as the curves III for intermediate values of a; there is 
* See Marshak (1947b). Actually the method that had been adopted by Mark differed slightly 
from that described by Marshak. Namely, instead of the boundary conditions used by Marshak, 
Mark used the alternative boundary conditions, outlined in the footnote on p. 444 of Marshak’s 
paper. In determining even order spherical harmonics approximation, Mark was allowing 
discontinuities in the even order moment at the boundary, but was stipulating that the discontinuity 


should be proportional to the infinite part of the derivative of the true value of the corresponding 
moment. 
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very little doubt that the curve IV represents the dependence of A upon a quite 
accurately. 

Figure 2 refers to the case of a black cylinder. Here the curves I and II 
correspond to the expressions for A derived in Part III, and in Kushneriuk’s 
paper, while the curve III represents the result of calculating A by the spherical 
harmonics method in P; approximation. Though higher order spherical 
harmonic approximations have not been calculated for the case of a black cylinder, 


0 1 2 3 7 5 
Radius C 
Figure 1. Extrapolation length A, in terms of radius of the sphere: I, small black sphere; I], 
large black sphere; III (a) (6) (c) (d), spherical harmonics (Ps, P;, Ps and P,; approximations 
respectively); IV, probable shape of dependence for a sphere of finite radius. Units in terms 
of mean free path. 


0:8 


0-6 
0 i 2 3 4 ei 
Radius 


Figure 2. Extrapolation length A, in terms of the radius of the cylinder: I, small black cylinder; 
II, large black cylinder; III, spherical harmonics (P; approximation); IV, probable shape of 
dependence for a cylinder of finite radius. Units in terms of mean free path. 


the comparison with the spherical case indicates that shape of the A(a) curve for a 
comparable with unity should be given quite accurately in the P, approximation. 
The comparison of the curves I, II and III then leads one to expect that the correct 


dependence of A on a for the case of a black cylinder should be as indicated by the 
curve LV of Figure 2. 
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Proof of the Inequality (1.16) 


Consider simultaneously with (1.5) the equation 
1 mc!) 
(x)= 5 | Wy) Bey!) dy + E(x) 


imposing upon h(x) the condition h’(oo)=0, similar to (1.6). Both h(x) and 
Gn(*) are, in Hopf’s terminology, N-solutions, and by alemma proved by Hopf 


(0) =f) =| Qala) BO) de = 5 FC) Ma) de. eee, (A1) 


On the other hand A(x) is obviously of the form h(x) = E(x) + a bounded function. 
Then proceeding as in the derivation of (1.13) and (1.15) we can easily show that 
the right-hand side of (A 1) is of the order of (a~-”"!log a). This establishes (1.16) 
forx=0. ‘To prove this for x40 from (1.5) we obtain 


1 ce) 
dn(*) fal) =[40(0)-Fa(O)] + 5] dal) B(x yl) -EO)] dy. 
On the grounds of the estimates obtained previously, this means 


an(e)—falt) =O(*™ Hoga) +a" Ol |” {LE(#—y)- 20) dy | 


and (1.16) then follows at once. 


AL PP EN-DE XB 
The Order of Magnitude of qo2(x) and q5o(x) 

Notice that 

Fr[|x—y]; (2ax+ x2)!” + (2ay +y?)¥?] < Fo{(2ax + x)? + (Zay + y2)¥? ; 
(2ax-+ x*)12 + (2ay + y?)1?} = E{(2ax + x*)1? + (2ay +y?)¥?}. 

Then the free term of the equation (3.11) is of the same order of magnitude, 
for each particular x, as f,(x) in Part I. Also, with the free term decreasing 
exponentially as x tends to 00, the stipulation that q..(”) should not contain a term 
proportional to x for large x is equivalent to the stipulation that qo9(x) is a 


bounded solution of (3.11). Then the order of magnitude of q,.(x) and of the 
iterated term in (3.11) should be also identical with the order of magnitudes of 


gx(x) and of al qi(y) E(|x—y|)dy of Part I respectively. From the analysis 
0 


given in Part I we see then that q..(x)=O(a*loga) for x>a' and that the 
iterated term in (3.11) is of the order of a? log a for any value of x. This proves 
the estimate of g..(«) quoted in connection with (3.11). 

To prove the estimate quoted for qs9(x) put q3o(*) = Gga1(*) + Gg20(*%) 1n which 
9321(*) 1s defined by the equation (3.12) with q,(y) in the integrand of the free term 
replaced by {q2;(y) +the iterated term of (3.11)} and q35.(x) is defined by the same 
equation with q,(y) replaced by the free term of (3.11). Now it is obvious that the 
free term of the equation for g32;(x) is of the order of a~* log a as compared with the 
free term of (3.11) and consequently q3.,(”) = O(a log @) qoo(*) = O(a~* log?a) 
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except for x comparable with a~', while the free term of the equation for q390(x) 
is of the same order as 


[| B(Qax +32)! + (2ay+ 9")? Bf(ay+ 9")! + 2aa+ a) dy da 


which for any particular x is of the order of a~! as compared with the free term of i 
(3.11). Then q392(x) = O(a~1 qo9(x)) and this proves the estimate of q39(x) quoted in 
connection with (3.12). . 
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The Excitation Function for the Production of “Be by the 
Bombardment of ’C by Protons 


By J. M. DICKSON anp T. C. RANDLE 
Atomic Energy Research Establishment, Harwell, Didcot, Berks. 


Communicated by D. W. Fry ; MS. received 9th February 1951, and in amended 
form 15th May 1951 


ABSTRACT. ‘The excitation function for the production of 7Be from '2C has been 
determined by an activation method, using the internal proton beam of the 110-inch 
Harwell cyclotron. The threshold energy has been found to be 32 Mev. and the cross 
section at 156 Mev. is 11:0+1-3 x 10-2? cm?. 


Sot LN ERO DU ChION 

N a previous paper (Dickson and Randle 1950) evidence was presented for the 
| production of 7Be by proton bombardment of #C, and a preliminary value 

of the cross section was given. Further experiments have been made by 
bombarding carbon targets at various radii in the internal beam of the Harwell 
cyclotron, in order to determine the excitation function up to 156 Mev. 

Although it is convenient to refer to the process as the reaction 2C(p, 3p3n)’Be, 
the protons and neutrons are not necessarily ejected as six separate particles. 
It will in fact be shown that they are certainly not ejected in this way at the lower 
proton energies. 
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§2. EXPERIMENTAL METHOD 
The carbon targets were semi-circular discs } in. thick by ?in. radius, and 
were known to contain very little impurity. For each irradiation one of these 
discs was mounted with its diametral edge vertical and facing towards the centre 
of the cyclotron. In addition to 7Be, 'C was produced in the target by the 
reaction *C(p, pn)#4C. Since the cross section for this reaction has been 
published, from the threshold (18-5 Mev.) to 350Mev., by Aamodt, Peterson 
and Phillips (1949), the “C activity could be used to monitor each bombardment. 
After irradiation the target was enclosed in an aluminium sheath, whose 
thickness was twice the range of positrons from UC, and the activity was measured 
by an end-window Geiger counter. The radiation detected consisted of 
annihilation radiation from the positrons of !C (half-life 20 minutes) and nuclear 
y-radiation from “Be (half-life 53 days). Activities with half-lives of a few hours 
were also detected, but the counting rates due to these were low. From the 
_decay-plot of y-activity against time the counting rates N(1'C) and N(‘’Be) due 
to UC and “Be at the end of the irradiation were calculated. These are related 
to the disintegration rates D(#!C) and D(’Be) as follows: N(?2C) =2D(UC)EHC), 
N(‘Be) = RD(‘Be)E(’Be) where E(4C) and £(*Be) are the efficiencies of the 
Geiger counter for the two y-rays and R is the branching ratio of the 7Be decay. 
Since the periods of irradiation were comparable with the half-life of 'C 
a correction factor {1—exp A(1'C)¢}/A(7Be)t was applied to the observed ratio of 
the two activities, where A(MC) and (Be) are the disintegration constants of 
UC and “Be and ¢ is the duration of the irradiation. The cross section o(*Be) 
for the production of 7Be was evaluated relative to the cross section o(14C) for 
the production of "C by the formula 
1—exp A(UC)t N(*Be) ECC) 2 


ee) Chey NCEA RB 


i The value assumed for the branching ratio R was 0-123 and the counter 
efficiency ratio, E(#1C)/E("Be) was 1-04. 


Effective Proton Energy 


The thickness in Mev. of the } in. thick carbon targets varied from 2-3 Mev. 
for 175 Mev. protons to 9 Mev. for 35 Mev. protons. ‘The beam, however, traversed 
the target several times (Cassels, Dickson and Howlett 1951), causing an even 
wider spread in bombarding energy. he effective thickness of the targets was 
18-5 Mev. at 175 Mev. and 13-5 Mev. at 35 Mev. ‘The cross sections were evaluated for 
the effective proton energy, i.e. the bombarding proton energy minus the effective 
target thickness. A few % in. thick targets were irradiated in such a way that the 
beam was prevented by a shutter arrangement from traversing the target more 
than once; all those protons which had been scattered by more than a very small 
angle during the first traversal of the target were absorbed by the shutter. These 
irradiations, at energies from 35 to 50 Mev., were used to determine the threshold 
energy. 

Half-life of UC 

The “C could not be counted until 4 to 5 hours after irradiation on account 
of its initial high intensity, but the counting rate after five hours could be 
estimated to an accuracy closer than 1%. Several short irradiations of carbon 
and polythene, producing relatively weak sources of “'C with negligible quantities 
of any other activities, were used to measure the half-life of "tC. The average 
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value of six runs, each extending over at least six half-lives, was 20-0 + 0-1 minutes. 
Values of 20-0, 20-35+0-08 20:5+0-6 minutes are given in the literature by 
Siegbahn and Bohr (1944), Smith and Cowie (1941) and Solomon (1941), but 
it was decided to use the value 20-0 minutes to calculate the initial ''C activity 
from the counting rate after 5 hours. |The error introduced in each measurement 
of the 7Be cross section, by an error in the half-life of “C, will not affect the shape | 
of the excitation function but only the absolute values of the cross section. 


Half-life of “Be 
Measurement of the half-life of 7Be gave the value 53 +2 days; the generally 
accepted value of 52:93+0-22 days, measured by Segre and Wiegand (1949), 
was used. 
Branching Ratio 


7Be decays by K-capture either to 7Li in the ground state or to “Li in the 
480 kev. excited state with the subsequent emission of a y-ray. The branching 
ratio R for this process has been measured by Turner (1949) and by Williamson 
and Richards (1949). Their results were 0-:118+10% and 0-107+0-02, 
respectively. Turner’s work was less liable to systematic error since it involved 
comparing counting rates and computing a counter efficiency ratio for two y-rays 
(0-51 and 0-48 mev.), while the measurement by Williamson and Richards involved 
two absolute calibrations of counters. Turner, however, used the value of 
20-5 + 0-6 minutes for the half-life of "C. He also had to find the ratio of the 
efficiencies of his counter for annihilation radiation and 0-48 Mev. y-rays. When — 
the latest values for the mass absorption coefficients (Snyder and Powell 1949) 
are inserted in Turner’s formula for the efficiency ratio, together with the value 
of 20-0 minutes for UC half-life, the value of the branching ratio becomes 
0-123+5%%. ‘This value has been used for evaluating our results. 


Counter Calibration 


The counter efficiency for 0-51 Mev. y-rays was measured by counting the 
annihilation radiation from a thin !C source and later counting the positrons 
in a defined solid angle f-counter. The value found was 6:0x10-4*. The 
efficiency for 0-411 Mev. y-rays was measured by counting the y-rays from a 
standardized 1%Au source, and found to be 5:2x 10-4. By interpolating these 
results the counter efficiency ratio, H(11C)/E(7Be) was calculated to be 1:04+5%. 


$3. RESULTS 

Twenty-five carbon targets were bombarded at effective energies between 
35 Mey. and 156 Mev. to determine the excitation function. Six carbon targets 
were bombarded for very short periods to determine the half-life of MC. 

The yields of "C and “Be in the 25 carbon targets per min. of irradiation 
were about 8 x 10’ and 500 counts per min. respectively at 156 Mev. effective 
proton energy, falling to 2-5 x 107 and 100 counts per min. at 50 Mev. 

The excitation function is plotted in the Figure. The yield of 7Be rises 
sharply above 32 Mev. The cross section is very small (less than 0-1 x 10-?? cm?) 
at 30 Mev., reaches a maximum of 13-8 x 10-27cm? at about 55 Mev. and falls 
slowly to 11 x 10-°7cm? at 156 Mev. 

The threshold energies for the reaction, calculated for all combinations of 
(3p3n), range from 24-4 to 59-0 Mey. (see Table). It is evident that the (3p3n) are 
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not produced as six separate particles and that the first three groupings in the 
Table are the only possible products of the reaction near the threshold energy. 


o (cm?) 


0 20 40 60 80 100 120 140 160 
Proton Energy (MevV.) 


Excitation function for the production of *Be by the bombardment of !2C by protons. 
The deviations shown in the figure are the compounded sum of errors due to counting 
statistics and errors quoted for the #C(p, pn)!!C cross section. In addition there is an 
uncertainty of 8% which can only affect the curve as a whole, without changing its shape. 
This arises from uncertainties in (1) 7Be branching ratio, (2) #4C half life and (3) the y-ray 
efficiency ratio of the Geiger counter. When combining this error with those shown in the 
diagram, the normal square law should be used. 


Products of reaction O (Mev.) Threshold (Mev.) 
"Be +S —22-6 24:4 
*Be+ *He+?H — 24-1 26:1 
"Be +*He+p-+n . — 26:3 28°5 
*Be+ ?He+ "2H — 38:5 41-7 
*Be+°H+2p+n —46-2 50-0 
7Be+ 9H+p+2n —46°9 50:8 
*Be+-"H--2H-+p —44-0 47-6 
"Be + 32H — 48-0 52-0 
*"Be+3p+3n — 54:6 59-0 
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ABSTRACT. The gamma-ray spectra from the nuclear reactions °B-+d, “B+d, and — 
27Al+-p have been measured with a pair spectrometer. Lines are found from B+d at 
8-88+ 0:06 mev., 6°7+0:15 Mev., 6:-4+0-15 Mev., and 4:38+0-05 mev.; from “B+d at 
4-44+0-05 mev.; and from ?’Al+p at 12:12+0-1 mev., 10-46+0-07 Mev., and 
7-62+0-1 Mev. The correlation of these gamma-rays with the energy levels of the nuclei 
involved is discussed. 


Silo UN AOD NUKE AN WON 
HE study of the energy and intensity of the high energy gamma-rays emitted 
in nuclear reactions has received comparatively little attention because of 
the lack of suitable techniques, particularly in cases where the gamma-ray 
spectrum is complex. However, the development of a pair spectrometer by 
Walker and McDaniel (1948) has made it possible to measure the spectra of 
high-energy gamma-rays with greatly improved resolution and efficiency. 

We have, therefore, applied this technique to the measurement of the 
gamma-ray spectra from a number of nuclear reactions for which our particle 
accelerator, which supplies a comparatively large beam current at potentials up 
to 800 kv., is particularly suitable. The results of an investigation of the proton 
capture radiation from !°F have already been published (Rae et al. 1950). 
Measurements on the radiations produced when 1°B and ™B are bombarded by 
deuterons, and ?’Al by protons, are presented in this paper, together with a 
description of the experimental technique. 


§2. APPARATUS 
(1) Generator 
A cascade generator was used to accelerate protons or deuterons to a maximum 
energy of 800kev. A magnetically resolved beam current of up to 150 ya. was 
available at the target. A complete description of this generator and its 
radio-frequency 1on source will be published shortly. 


(ii) Target Arrangements 


The general arrangement of the target assembly and pair spectrometer is 
shown in Figure 1. The resolved beam from the accelerator tube emerged from 
the resolving chamber down the axis of a glass tube inclined at 30° to the vertical. 
A shallow water-cooled brass cup was sealed by a demountable rubber gasket 
joint to the end of this tube. ‘T’argets were normally prepared on thin brass or 
copper discs, which were clamped to the inside of the cup. The cup was insulated 
from ground and was connected to the input of a current integrator of conventional 
design. 


Gamma-Ray Energy Measurements QO7 


A Geiger counter with thick aluminium walls was mounted at a fixed distance 
from the target and was used to monitor the gamma-ray intensity. All runs were 
made for a fixed number of counts in this counter. The current integrator 
readings were used mainly to check target deterioration and to estimate the 
relative yield of the various reactions. 


(iii) Pair Spectrometer 


The vacuum chamber of the spectrometer consisted of a circular copper 
cylinder 12 in. in diameter and 21 in. long with circular steel end plates 14 in. 
in diameter and 4 in. thick sealed to the cylinder on one side with a soldered joint 
and on the other by a rubber gasket. The gamma-rays from the target were 
collimated by the lead shielding blocks and allowed to fall on a thin lead foil 
placed at the centre of the vacuum chamber. This foil was mounted on a light 
frame which could be rotated about an axis through one edge by means of a rod 
passing through a rotary vacuum seal. Thus, it was possible to take measurements 
of the background coincidence counting rates by rotating the foil through 90°, 
so that it lay against the vertical end plate of the vacuum chamber. These 
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Figure 1. General Figure 2. Block schematic diagram of spectrometer and 
arrangement. recording circuits. 


background rates were negligible when measurements were made on proton 
reactions, but were often serious in deuteron reactions because of the difficulty 
of providing adequate screening of the counters from the gamma-rays produced 
by the accompanying neutrons. 

The chamber was placed between the 12 in. diameter poles of an electromagnet 
of conventional design which provided a field of up to 7,000 gauss. This field 
was sufficient to measure gamma-ray energies up to 30Mev. ‘The field was 
measured by means of a flip coil and fluxmeter. The relative accuracy of this 
measurement was estimated to be within approximately 0-3°%. An absolute 
calibration of the instrument was made by measuring the spectra of the gamma-rays 
from fluorine and lithium under proton bombardment and will be described in 
detail later. A magnet current of 50 amperes at 50 volts was required for the 
maximum field and this was supplied by a D.c. generator. The current was 
stabilized and controlled to an accuracy of within 0:1% by a feedback circuit 
operating on the field current of the generator. 

The electron—positron pairs which were ejected predominantly normal to 
the foil by the gamma-rays were deflected through 180° by the magnetic field 
and detected by two sets of five Geiger counters placed in the plane of the foil. 
The counters were screened from direct gamma-radiation by the lead blocks. 
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These counters, which are described elsewhere (Rae 1950), were constructed 
from brass tube, 4 in. in diameter and of 4 in. wall thickness. The tube walls 
were cut away over 180° along their operating length and replaced by 0-001 in. 
copper foil to allow electrons or positrons to enter the counters with a minimum 
of scattering. The counter wires were supported on Kovar seals soldered to 
one end of the counter tubes. The other ends of the tubes were sealed up and 
provided with short helical springs which were located in blind holes recessed 
in the side plate of the spectrometer. ‘The ends provided with Kovar seals 
projected through clearance holes in the other side plate of the spectrometer and 
were sealed to it with ‘Q’ compound. A gap of Lin. was left between this plate 
and the magnet pole to allow connection to be made to the counters. 

Baffles were placed in the vacuum chamber, as shown in Figure 2, to reduce 
the background due to electrons produced by gamma-rays from (n, y) reactions 
in the counter walls and the lead shielding behind them, which could otherwise 
traverse the spectrometers and cause spurious coincidences. 


(iv) Recording Equipment 

The output pulses from the counters were fed via cathode followers and 
shaping circuits to twenty-five coincidence circuits each with a resolving time 
of about 2 .sec. The arrangement of these circuits is shown in the schematic 
diagram of Figure 2. The output pulses from the coincidence circuits were 
grouped into nine different energy channels, and the output of each channel 
was recorded on a mechanical counter. 

The pulses from the monitor Geiger counter were recorded on a scale of 100. 
The output of this scaler was fed into a master control unit. This unit contained 
relays which switched on all the recording units and scalers simultaneously at 
the turn of a switch. It also contained a counter which could be preset to switch 
off all the recorders after any desired number of counts from the monitor counter 
scaler. Thus, all runs could be made for a fixed number of monitor counts. 
The time of each run was also recorded on a counter operating from a clock which 
gave pulses at half-second intervals. 


§3. EXPERIMENTAL PROCEDURE 

In measuring a particular gamma-ray spectrum a series of runs was made at a 
sequence of values of the magnetic field to cover the desired energy range. Each 
magnetic field setting gave counts in nine different energy channels, the spacing 
between channels corresponding to 5-8°% of the Hp value of the centre channel. 
The interval between the successive values of the field was usually of the order 
of 10% or less, so that there was a considerable overlap between the runs. The 
results were then normalized for channel number and flux and plotted as a 
function of the mean Hp of a single particle. The precise energy of a particular 
gamma-ray was determined from the Hp value of the half-amplitude point on 
the high-energy side of the peak, as this could be measured more accurately 
than the Hp value of the peak itself and also should be less dependent 
upon the effects of scattering and energy loss in the foil. An analysis similar 
to that of Walker and McDaniel (1948) showed that a relation of the form 
k=600 Hp(1+2-5 ?/k*) should be accurate to within 0:5% over the 
energy range covered by this work, where k is the energy of the gamma-rays 
and yz the rest energy of an electron. No absolute calibration of the fluxmeter 
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coil was made. The spectra of the 6-13 Mev. gamma-rays from the 340 kev. 
resonance of the reaction !°F(p, «, y)!6O and of the 17-6 Mev. gamma-rays from 
the 440 kev. resonance of the reaction ’Li(p, y)&’Be were measured and used. to 
calculate the conversion factor between fluxmeter readings and Hp values. The 
values of the conversion factor so obtained from the two measurements differed 
by only 0-5°% and the mean value was used in the calculation of subsequent 
results. 

The spectrum of the 6:13 Mev. gamma-rays from F(p, «, y)!®O is shown in 
Figure 3 to illustrate the resolution of the spectrometer. 
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§4. EXPERIMENTAL RESULTS 
(i) The Reaction B+d 


A target of natural boron was prepared on a brass plate and bombarded with 
600 kev. deuterons. The measured gamma-ray spectrum is shown in Figure 4. 
At the time of this measurement the spectrometer had not been fitted with the 
baffles described in § 2 (iii) and the background coincidence rate due to neutrons 
was of the same order of magnitude as that due to pairs produced by gamma-rays 
in the foil. This background was subtracted to obtain the true spectrum shown 
in Figure 4, but, as a result, the statistical errors were rather large. 

A target of 1°B, prepared in the Harwell electromagnetic separator, was then 
bombarded with 600kev. deuterons and the consequent gamma-ray spectrum 
is shown in Figure 5. The background was considerably smaller on this run, 
partly because of the very much smaller neutron flux from the 1°B(d, n)4C 
reaction as compared with that from “B(d, n)#C, and partly because of the 
baffles which had been fitted to the spectrometer. The gamma-ray yield was 
also greater because of the use of an isotopic target. The statistical accuracy was, 
therefore, very much better on this run than on the previous one. 

A third run was then made with a target of separated "B under similar 
conditions and the gamma-ray spectrum is shown in Figure 6. The background 
on this run was of the order of magnitude of one third of the true coincidence 
rate at the peak of the spectrum and has been subtracted to give the spectrum 
shown. 
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The results of these three runs are summarized in Table 1, together with 
the measurements of previous workers on the same reaction. It will be seen 
that the majority of high energy gamma-rays from B+ d reactions are due to the 
10B isotope and that 4B+d does not produce any gamma-rays of energy greater 
than 4-44mev. The spectrum from 1°B+d is clearly complex in the region 
between 6 and 7 Mev. and is not completely resolved. However, both the B+d 
and 1B+d runs indicate two gamma-rays of the same order of intensity at | 
6-4 and 6:7mev. There is also a pronounced ‘tail’ on the high-energy side of | 
these peaks which suggests the presence of a gamma-ray of smaller intensity — 
at about 7 Mey. 7 
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Figure 5, Gamma-rays from °B-+d at Eg=600 kev. Figure 6. Gamma-rays from 14B+d at 
E4=600 kev. 
Table 1. Summary of Gamma-ray Measurements 


Halpern and Crane 


Ree Present work Gaerttner et al. (1939) (1939) 
E, (Mev.) Rel. int. £,(Mev.) Rel.int. E,(Mev.) Rel. int. 
8-88 + 0:06 1 S104. il 8-6 1 
6:7, +0°15 0- 
Bad pee ear eOL 04s 6-0 2 
438+ 0-05 2°9 4-4+0-3 10 4-2 6 
8-88 + 0-06 1 
Ord sees 0:8 
10 = 
Ea 6-4 40-15 1-0 
4-50+ 0-05 0-9 
uB+d 4-44+ 0-05 
122122 0-4 1 
27Al+p 10-46+ 0-07 HS 
F-O2eta ei $305 


We can now consider the interpretation of these results. In the case of 
1B +d, three reactions are energetically possible as shown in the energy level 
diagram, Figure 7 (a) (Hornyak et al. 1950). These are °B(d, p)"B, O = 9-24 Mev. ; 
10B(d, «)®Be, O=17-81 Mev. and !B(d, n)"C, O=6:53 ev. 
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The reaction !B(d, p)"B leads to a number of excited states in™B. The 
Most precise measurements of the energies of these proton groups have been 
made by Buechner and his co-workers using 1-5 Mev. deuterons and his results 
| are shown in Table 2, taken from the recent review by Hornyak et al. (1950). 


Table 2. Energy Levels of “B from the B(d, p)“B Reaction 


Proton measurements + Gamma-rays 


ae sone ~ Rutherglen and Smith Present work 


E (level) E (level) Rel. int. E,, (Mev.) Rel. int. 


(Mev.) (Mev.) 
0-000 
2-141 
4-457 4-50 1-9 4-50 0-9 
5-033 5-01 0-3 
6-4 1-0 
6-752 6-72 Of 6-7 0-8 
6-802 
7-297 7-26 0-5 
8-565 8-51 0-1 
8-921 8-93 1-0 8-88 1-0 
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We also include in Table 2 the results of measurements on the relative intensity 
of some of these proton groups obtained with a magnetic spectrometer in this 
laboratory (Rutherglen and Smith, to be published) using 600 kev. deuterons. 
It will be seen that the energies of the measured gamima-rays at 4:5, 6-7 and 
8-9 Mev. are in good agreement with the energy levels predicted by the three most 
intense proton groups. The intensities of the remaining proton groups are too 
small to give observable gamma-rays in this work, but the proton group 
corresponding to the 7-3Mev. level may be responsible for the shape of the 
gamma-ray spectrum on the high energy side of the 6-7 Mev. peak. 

The !°B(d, «)®Be reaction is known to lead to levels in *Be at 7-5 Mev. and 
2-9 mev. (Smith and Murrell 1939) but it is almost certain that these levels decay 
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to the observed gamma-radiation. 
The neutrons from 1B(d, n)#4!C have been studied by Gibson (1949) at 
E, =0-93 mev. and by Swann and Hudspeth (1949) at Ej3=1-4Mev. The former 


find neutron groups corresponding to the ground state and to a 2:0 Mev. excited — | 


state in "C. The latter find evidence for only two levels in #1C at 4-5 Mev. and 
6:7 Mev., but there is a possibility that the corresponding neutron: groups were 


due to contamination. If a level exists in the region of 6-7 Mev. it may be 


responsible for the complexity of the gamma-ray spectrum observed in this region. 
We may summarize these reactions by saying that the majority of the high 


energy radiation from !°B + d arises from transitions between excited states of "B | 


and the ground state, following the 1°B(d, p)"B reaction, but that there may be 
some gamma-radiation of energy about 6-4 Mev. from the 1°B(d, n)""C reaction. 
In the case of “B+d, we must consider the reactions “B(d, n)?C, 


OQ =13-78 mev.; “B(d, p)2B, Q=1-1Mmev. and ™B(d, «)®Be, Q=8-03 Mev., as 


shown in the energy level diagram of Figure 7 (6) (Hornyak et al. 1950). The first 
reaction has recently been studied by Gibson (1949) and is found to yield neutron 
groups leading to excited states of #C at 9-72+0-15 Mev. and 4-47-+0-10 Mev. 
with relative intensities of 0-67: 1. ‘The second group corresponds well with the 
gamma-ray energy of 4-44+0-05 Mev. measured in this work. However, there 
is no evidence of a gamma-ray at 9-7Mev. This is not unexpected, since this 
level can decay by the emission of an alpha-particle to the ground state of *Be, 
with an energy release of 2-3 Mev., and one would expect the lifetime for this 
process to be very much shorter than that for the emission of a gamma-ray. 
Malm and Buechner (1951) have also observed this level in the reaction ® N(d, «)¥C 
and find that it has a greater width than the lower levels in ?C, which is consistent 
with the assumption that it decays by alpha-emission. The possibility of a cascade 
transition through the 4-5 Mey. level seems to be ruled out by the fact that a 
gamma-ray of energy 5-2 Mev. is not observed. 

The reaction /B(d, p)”B need not be considered because of its low Q-value. 
The reaction “B(d, «)®Be has recently been studied by Van Patter and Buechner 
(Hornyak et al. 1950) and has been found to lead only to the 2-4 Mev. level in °Be. 
It would, therefore, not be expected to contribute to the gamma-radiation 
observed in this work. 

We may summarize these conclusions by saying that the high energy 
gamma-radiation from “B+d arises from the transition from the 4:47 Mev. 
level in #C to the ground state and that the 9-7 Mev. level in ®C decays by 
alpha-particle emission to ®Be. 


(ii) The Reaction ?’Al+p 


The tneasured gamma-ray spectrum from a thick aluminium. target 
bombarded with 750key. protons is shown in Figure 8. The energies and 
relative intensities of the main peaks are shown in Table 1. The thin target 
excitation curve of the gamma-rays from Al+p has been measured by a number 
of workers and shows a number of sharp resonances. In order to obtain 
sufficient intensity in the present work it was necessary to use a thick target, 
so that all the resonances below 750 kev. were excited. Even so, the total 
gamma-ray yield was only about one-fifteenth of that produced by the 340 kev. 
resonance of °F(p, «, y)"®O. It will be seen that three gamma-rays of energies 


| 
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: 12-12 +0-1 Mev., 10-46 +0-07 Mev., and 7:62 +0-1 Mev. are clearly resolved, but 
the background both above and below the 7-6 Mev. peak suggests the presence of 
other weak gamma-rays. ‘The low intensity of the reaction prevented a more 
detailed investigation of this region. 

_ The interpretation of these gamma-rays can be discussed in terms of the 


energy level diagram of ?8Si shown in Figure 9. The only reactions which are 
80 
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Figure 8. Gamma-rays from ?7Al+p at Ep=750 kev. 


si” 


Figure 9. Energy level diagram for ?’Al-+p. 


energetically possible in this case are ?’Al(p, y)?8Si and ?7Al(p, «)?4Mg and only 
the former has a sufficiently large Q-value to produce the observed gamma-rays. 
The most energetic gamma-ray is presumed to correspond to the direct transition 
to the ground state of *Si. From the results of Brostrom ez al. (1947) on the 
relative intensities of the resonances below our bombarding energy of 750 kev. 
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we assume the mean effective bombarding energy to be 630 kev., which corresponds 
to an additional excitation energy of 610kev. in centre-of-mass coordinates. 
Although it is known that the spectrum is not identical for all resonances (‘T'angen 
1946), the above assumption cannot be in error by more than +100 kev. We 
thus obtain a value of 11:51+0:2Mev. for the Q-value of the ground state 
transition. This Q-value may be calculated by two other methods. First, the 
reactions 27Al(d, p)?8Al, O =5-45 + 0-05 Mev. (Pollard et al. 1949) and.*®Al(8")?*Si, 
O=4-80 + 0-05 Mev. (Benes et al. 1948) predict a Q-value of 11-70 + 0-1 Mev. for 
27Al(p, y)?8Si. Second, the reaction ?7Al(d, n)?&Al, Q=9-08+0-2mev. (Peck 
1949) predicts a Q-value of 11:31+0-2mev. It will be seen that there is no 
serious discrepancy between these three determinations and they yield a weighted 
mean value of 11:55 + 0-07 Mev. . 

The gamma-ray of energy 10-46+0-07 Mev. is almost certainly due to a 
cascade transition which proceeds through a level of 78Si in the region of 1-8 Mev. 
On this assumption our measurements give a value of 1:66+0-2 Mev. for this 
level. The reaction ?8Al(8~)?8Si proceeds via a level at 1-86 + 0-05 Mev. (Benes 
et al. 1948) and the reaction ?’Al(d, n)?8Si proceeds via a level which has been 
measured as 1-:78+0-13 Mev. from the neutron groups (Peck 1949) and as 
1-72+0-08 Mev. from the gamma-radiation (Alburger 1949). Thus, all the 
determinations are consistent with a weighted mean value of 1-75 + 0-04 mev. for 
the energy of this level in 28Si. 

The interpretation of the 7-6 Mev. gamma-ray is less certain, but it seems_ 
probable that it is due to a cascade transition which proceeds through a level 
in 78Si at 4:50+0-15 Mev. ‘This interpretation is consistent with the level at 
4-47 + 0-13 Mev. found by Peck (1949) in the ?”Al(d, n)?8Si reaction. 
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ABSTRACT. ‘The ratio of penetrating particles to electrons in extensive cosmic-ray 
showers has been investigated experimentally at sea level. For the great majority of showers 
this ratio, in agreement with previous results, is found to be 2°5%. A small number of 
showers was detected in which the proportion of penetrating particles was much greater. 


§1. INTRODUCTION 

T is well known that the particles of the extensive air showers consist 
| predominantly of electrons but that there is also present a small fraction of 

more penetrating particles. The ratio k of the density of penetrating particles 
to electrons has been found to be of the order of 2°% (see §6 for references), but 
evidence for a variation in k has also been found within a limited range of k values 
(Cocconi, Tongiorgi and Greisen 1949 a,b). The present experiment was carried 
out in order to detect such variations if they exist over a wider possible range of 
values. 


§2. THE EXPERIMENTAL ARRANGEMENT 

Three trays of ten Geiger-Miiller counters, of area 1,750 cm?, were situated at 
the vertices of a triangle of side 5 metres, housed in a wooden hut 4 metres below 
a roof of mean thickness 5gm/cm?. For the greater part of the experiment these 
trays were shielded by 20 cm. Pb above, with 15cm. Pb at the ends and sides and 
2-5cm. Pb below. When a run was made with the trays unshielded, they were 
placed above the lead piles with a layer of wood 8:5 cm. thick between the counters 
and the lead to avoid a possible spurious increase of the density of the particles at 
the counters by back scattering from the lead. When unshielded, the counters 
were covered by an aluminium box 1 mm. thick. 

In addition to these three trays which provided a threefold coincidence as 
master pulse, an array of six counter trays of areas 1,400 cm?, 700 cm?, 350 cm?, 
175cm?, 80cm? and 37cm?, enclosed in an aluminium box 1mm. thick, was 
situated within the triangle formed by the master trays. The four largest of these 
trays consisted of 8, 4, 2 and 1 counters respectively, identical with the counters 
comprising the master trays. ‘The remaining two consisted of single counters of 
smaller size. All the counters were separated from each other by 0-5cm. Each 
of the six counter trays was connected to a neon lamp in a hodoscope so that by 
examination of the photographic records it could be determined which of the 
hodoscope trays had been struck by each of the extensive showers detected by the 
three master trays. 

The resolving time of the threefold coincidence circuit providing the master 
pulse was 0-6 microsecond so that in the total time of recording less than one 
casual coincidence would be expected. 
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§3. ANALYSIS OF THE HODOSCOPE DATA 
(i) Notation 

It is convenient to consider one fortieth of the area of the master trays as the 
unit of area in analysing the hodoscope results. The areas of the hodoscoped 
trays are then 32, 16, 8, 4, 1-8 and 0-85 units, while the master trays are 40 units of 
area. For ease in calculation the areas of the two smallest hodoscoped trays were 
taken as 2 and 1 units. Errors involved by this approximation will be negligible. _ 
It will be observed that any combination of the hodoscoped trays set off in 
coincidence with the three master trays may then be denoted uniquely by one of 
the integers from 0 to 63. For example, a coincidence in which none of the 
hodoscoped trays is struck is denoted by ‘0’, one in which all the trays are discharged 
by ‘63’; a‘ 48’ coincidence is one in which trays of areas 32 and 16 onlyare set off 
and no others, and so on. 

The numbers of coincidences recorded involving the various possible combina- 
tions of the unshielded hodoscoped trays will be given in the following two sections, 
with the shielded and the unshielded master trays respectively providing the 
shower selection. The observed distribution of the numbers of different types of 
coincidences will be compared with the predicted distributions. 


(11) The Shielded Distribution 


The expected counting rate of any arrangement of counters may be calculated — 
using the well-known integrals (Broadbent and Janossy 1948). This calculation is 
based on the following assumptions: (i) that the average particle density in a given 
shower is constant over the experimental arrangement, (ii) that the spatial distri- 
bution of particles is purely random, (iii) that the rate of incidence of showers with 
particle density equal to or greater than x per unit area is N(x) =Bx~’, y and B 
being constants, and (iv) that k, the ratio of penetrating particles to electrons, is 
a constant. 

In Table 1, column 2, are given the observed numbers of coincidences involving 
the particular combination of the hodoscoped trays indicated incolumn1. In the 
succeeding columns the calculated distributions of the various events, normalized 
to the observed total number, are given using the density exponent y=1-425 as 
found by Broadbent, Kellermann and Hakeem (1950) and for the values of k 
indicated. 

The combinations given in Table 1 are those of which two or more instances 
were recorded. It will be observed that there were seven other combinations 
which occurred to the extent of once each. ‘These have been grouped together in 
the category ‘others’ and are compared with the predicted numbers of events 
other than those listed. 

In order to find the calculated distribution which best fits the observed distri- 
bution, we evaluate the sum of the squares of the deviations of the observed numbers 
of events 0; of the different types tabulated from the corresponding calculated 
numbers ¢;, ‘Thus wu? =X(0,—c¢;)?.__ We may consider the best agreement to be 
afforded by that calculated distribution for which u? is a minimum. It will be 
seen from the values of wu? given in the last line of the Table that a minimum is 
reached corresponding to the calculated distribution for k=2°5/40=6:25%. The 
expected value of u? is given by u2 = N— Q/N where Nis the total number of events 
recorded and O=%c,?._ For a single observed distribution differing from the 
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predicted distribution only because of statistical fluctuation of the individual 


numbers of the events recorded, the deviation o of u” from the expected value u? has 
the mean square value o? given by* 


where K=2c,?._ Thus for the observed distribution to be in agreement with the 


Table 1. Shielded Distribution 


Values of k 
Counter NO." ot 0-5 foe dss 0 28S 30 S40 
combination events — — —— ees —— —— eee = 
40 40 40 40 40 40 40-40 
£63” 122 286 221 175 143 120 +102 ~# 62 4 
“62 53 44 CTs ose 7 
“61° 16 11 fo! 2.893 gg) gg i95 9g 3 
«60. 48 14 26M bis 6 DI 477 e500 al sae 
«58° 5 5 Tiitlie maniounw 1S roel 6 
57° 7 3 ee ee a a 3 
(56? 29 2 Oe Te 20 95 90 a 40 6 
54? 2 1 ee Se a ne 2 
«50° 4 1 ee ee ek ae es ed 5 
‘48? 20 3 Bee ch 6S Seay OR) Ld) 44s eS 
‘44? 2 oe 0 1 ae rene eee: 4 
«40° 7 1 Genesio atags) Teg iN) ate (ny eats 
36° 3 0 ari 1 7 Avcone, 33 8 
39° 13 0 ee Ve TT Peeples. Lets 
16° 7 0 0 uaed i Tae aie vary 
42? 2 0 Oe me LORS OR DeROra Os PANO 1 
<0 26 0 1 iene mee aN enn Cones ee oy 
Fathers’ 7 4 100 thes ote ag 90, al S60) 
Total 373 


ee BOs PON. Ah teeh PsAey arias) APS) Sysisks) Aah 7/s 


Table 2. Unshielded Distribution 


(1) (2) (3) (1) (2) (3) 
§63’ 9 12 “40” 46 51 
SOD Bil 20 aCON 14 23 
Gib: 16 8 34? 8 11 
£60’ 60 41 SB 182 180 
£58’ 25 iy *24’ 12 17 
9572 10 8 e207 7 8 
e50n 88 78 NG? 71 70 
£54’ 3 6 De + 3 
ASye 30 3) S28¢ 33 31 
50k 14 15 a2 9 15 
£48? 131 128 ‘others’ 48 62 
°44? 7 11 oan Oe 246 D5 


Total 1,104 ; uw?=1,266 ; expected value u2—979 +460 


(1) Counter combination ; (2) observed number of events ; (3) calculated number 
of events. 


* We are indebted to Professor Schrodinger for his investigation of the statistical problem here 
involved and for making available to us this expression. 
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calculated distribution with k=6-25°% one expects u? to be within the range 
u? + (o?)"* =307+182. In fact we find w2=1,778 for this case hence, while 
k=6-25 % may give the best agreement with the observed distribution, this agree- 
ment is unsatisfactory and, therefore, at least one of the initial assumptions (1) and 
(iv) used in the derivation of the calculated distributions cannot be strictly valid. 


(iii) The Unshielded Distribution 


Table 2 gives the observed and calculated values for the unshielded rates; 
it includes all combinations of which five or more instances occurred. As before 
the sum of the squares of the deviations was evaluated. The expected value is 
u® + (o?)"? =979 + 460, and observed value u2=1,266. Agreement between the 
observed and calculated distributions in this case therefore is satisfactory. 


§4. THE THREEFOLD COINCIDENCE RATES 
It is possible to obtain the fraction k by a number of different methods. One 
method has already been discussed. Other methods are to compare the rates of 
threefold coincidences of various combinations of shielded and unshielded trays. 
The three combinations used were (a) threefold unshielded, (b) two unshielded 
and one shielded, (c) threefold shielded, and the number of coincidences per hour 
observed were 
(a) Ray =33-9 + 0-97, (6) Roysig=2°36+0-10, (c) Rgg=0-1159 + 0-0057. 
The rates (a) and (b) have been corrected for barometric variations to the mean 
atmospheric pressure over the run (c), using a barometer coefficient of —11% 
per cm. Hg which is known to apply to observations with counter arrangements (a) 
and (b) (Millar 1950, 1951). From any two of the counting rates given above a 
value of kcan be deduced. ‘The values so obtained are as follows: from (5) and (a), 
k=2-55 + 0-10%, from (c) and (a), k = 1-86 +0-07 %, from (c) and (6),k =1:140-3%. 
It will be seen that these three values of k differ markedly, and even considering the 
first two values as two independent determinations of k, the difference between 
them is greater than can be explained by statistical fluctuation of the counting 
rates alone. Moreover, none of these values is in agreement with the value, 6:26%, 
giving the best agreement with prediction for the hodoscope data with the three 
shielded trays. 
The reason for these discrepancies must be sought in the original assumptions 
(1) to (iv) which provide the basis of the calculations, and these assumptions will be 
examined in the light of the observed data. 


§5. DISCUSSION OF THE EXPERIMENTAL RESULTS 
(1) Unshielded Data 


A calculation of the average densities triggering the different combinations 
shows that the range of densities (3 to 400 particles per square metre) is much the 
same as that investigated by Broadbent e¢ al. (5 to 500 particles per square metre). 
It will be seen that the predicted and observed rates are in satisfactory agreement 
when Broadbent’s value for y is used. 

There is some evidence that the observed number of high density events is 
larger than the calculated number. _ If this effect is genuine (none of the individual 
fluctuations is greater than 2-5 times the expected deviation) it is probably due to 
the variation of the shower density over the apparatus as suggested by the authors 
mentioned. 
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" (ii) Shielded Data 


The great majority of the showers falling on the apparatus have a mean value of 
k=2:55+0-10%. During the total running time with the three shielded trays 
about 100,000 showers capable of triggering the three unshielded trays fell on the 
apparatus. Of these 373 triggered the shielded arrangement and it is of these 
showers that a significant proportion show variations in R. 

The most significant difference between the predicted and calculated rates is 
that for the ‘0’ events. ‘Twenty-six showers were observed in which three 

shielded trays were discharged without an accompanying discharge of any of the 
_ unshielded trays. In a further thirteen events only the largest of the unshielded 
trays was struck. It would appear that these showers can only be explained by 
assuming a value of k of $ to 1. 

Considering now the high density end of Table 1 it will be seen that 122 events 
were recorded in which all trays were discharged. If k is equal to 2:5°% one 
expects 330 such events to be recorded. (The threefold shielded coincidence 
rate must be increased by 50%, from 0-116 to 0-174 per hour, to bring it into 
agreement with the threefold unshielded rate for R=2:5°%.) The average 
density here involved is greater than 1,000 particles per square metre and it is 
therefore possible that in this density region k is less than 2:5 %% or that y is greater 
than 1-425, both of which variations would lead to the recording of fewer high 
density showers than are expected for k constant and y having these values. 
However, it is also possible that for such high densities the variation of the density 
of incident showers over the apparatus becomes more important. 


§6. COMPARISON WITH PREVIOUS EXPERIMENTS 

Previous determinations of the value of k have in general involved one of two 
methods: (a) the comparison of the coincidence rate obtained from unshielded 
counters with the rate obtained from unshielded counters plus one or more 
shielded counters in coincidence, or (b) the comparison of the coincidence rates 
obtained with all counters unshielded and with all counters shielded. Using the 
first method Broadbent and Janossy (1948), Chowdhuri (1948) and Cocconi, 
Tongiorgi and Greisen (1949 a) obtain at sea level a value of k of approximately 
2-5 % which is in good agreement with our value obtained by the same method, i.e. 
comparing the twofold unshielded plus one shielded rate with the threefold 
unshielded rate. At about 3,000 metres altitude both methods have been used, 
the first by Cocconi, Tongiorgi and Greisen (1949 a), and the second by Treat and 
Greisen (1948) and Ise and Fretter (1949). As pointed out by Ise and Fretter 
the second method appears to give a higher value of k than the first. ‘This is in 
direct contradiction to our results, since we obtain a lower value of k from the 
second method (comparing our threefold shielded and threefold unshielded rates) 
than from the first. The results at high altitude, however, were obtained from 
experimental arrangements involving different areas, separations and dispositions 
of the counters and comparison of these results is complicated by this fact. 

Evidence that & is not constant but varies with the mean electron density to 
which an arrangement of counters responds has been deduced by Cocconi, 
Tongiorgi and Greisen (1949 a) who find that k is proportional to x. The 
problem of the variation of k in relation to the lateral spread of the showers has 
also been investigated by these authors. They find that at each of two distances of 
48m. and 96m. from the shower core detected by a ‘ core selector’ there is no 
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evidence that k is not constant for electron densities over a range varying by a 
factor of more than 100; there is evidence, however, that k increases with distance 
from the shower core and they deduce that the value may be twice as great at 
100 m. from the core as it is at 5m., a result which is in agreement with the weak 
dependence of k upon the mean particle density previously detected by these 
authors. . 

Such a slow variation of k with distance from the shower core could explain 
the discrepancies in Table i for the high density regions corresponding to showers 
detected near the core. The discrepancies in the low density regions, however, 
cannot be satisfactorily explained in the same way in terms of the lateral structure of 
the showers since, in particular, our observed ‘0’ events must arise from shower 
regions in which the density of penetrating particles is at least of the same order of 
magnitude as the density of the accompanying electrons. We are led therefore to 
interpret our low density events in terms of the longitudinal structure of the 
extensive air showers and conclude that with an arrangement of three shielded 
counters in coincidence an appreciable fraction of the recorded coincidences are 
due to extensive air showers in which the soft component has been largely absorbed 
in the atmosphere before reaching the apparatus, leaving only an extended tail 
consisting of non-interacting penetrating particles and, presumably, their 
accompanying knock-on electrons. 


§8. CONCLUSION 
It has been shown that the assumption of constant k is invalid from analysis of the 
extensive penetrating shower data here obtained. It is possible to interpret these 
results on the basis that the majority of showers have a mean value of R=2:5%, 
but it is also necessary to postulate that there exists a small number of showers with 
k approaching unity. If these are in fact the tail ends of extensive air showers 
whose soft component has died out above the apparatus as seems likely, it would 
lend support to the hypothesis that all cosmic-ray particles observed in the lower 
atmosphere are indeed parts of extensive air showers as recently suggested by 
Cocconi (1950). 
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ABSTRACT. The molecular orbital method is employed in calculations of the 
diamagnetic anisotropies of a series of highly condensed aromatic systems having hexagonal 
symmetry and containing up to fifty-four carbon atoms. A new method of calculation 
is introduced in which the z-electron susceptibility of the whole molecule is determined 
by an integration, no knowledge of the unperturbed energy levels being necessary : the 
inclusion of overlap between adjacent atomic orbitals is shown generally to have little 
effect on the calculations. 

The results are compared with those calculated using the method of Pauling: this 
semi-classical method is found to be quite misleading when applied to large molecules 
or to molecules possessing stray C—C bonds. Qualitative valence bond considerations are 
advanced in support of these conclusions, and certain general results are tentatively 
suggested. 


§i. REDUCTION OF THE SECULAR EQUATION: AN EXAMPLE 

N this paper the method of molecular orbitals is employed, in the manner 
] outlined in Part I and illustrated in Part II of an earlier paper (McWeeny 

1951, to be referred to as I and II), in a calculation of the diamagnetic 
anisotropy of various carbon structures containing as many as 54 carbon atoms. 
The secular equations involved in calculations of this type cannot often be 
dealt with by analytical methods (cf. II) but they can be made manageable in 
certain cases of high molecular symmetry by an application of group theory: 
for this reason our attention will here be confined to structures possessing 
hexagonal symmetry, in which case a great simplification can be achieved. An 
application to the benzene molecule affords an excellent example of the method 
to be used; we shall therefore consider this case in some detail, although the 
results (obtained otherwise) are well known (e.g. London 1937). ° 

The point group of the benzene molecule with a magnetic field along the 
6-fold axis is the cyclic group C,, the reduction of symmetry from the normal Dg, 
resulting from the exclusion of reflections by the field (cf. I, Appendix I). The 
character table for C, is given below. 


i E= 62) C G2 (ON Gs Ge 
ea 1 1 1 1 1 1 
Rey = 9% 1 —1 1 —1 1 —1 
Ts Xs 1 w uw w* w4 w? 
ES il w?> wt w? w? w 
PR Ge 1 aw? w4 1 w w4 
1B 1 w4 ow? 1 w* w 

cee (1) 


Since the group is Abelian, an irreducible representation by unitary matrices 
is necessarily unidimensional: consequently the representative ‘matrices’ D,(s) 
(s being a group operation and r denoting the rth representation, I',) are 
synonymous with the characters y,(s) given above. 
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In constructing molecular orbitals from the atomic 2p orbitals 4, ... . ¢¢ of 


the six carbon atoms it is possible to write down immediately six distinct linear — 


combinations providing a basis adapted to the complete decomposition of the 
six-dimensional (d-space) representation into the direct sum of six one-dimensional 
representations. For, from above, these combinations, #, .. . . , Say, (%, being 
the basis vector for the invariant (uni-dimensional) sub-space providing the 
representation I) must transform under the operations of the group according 
to o,>w,’ = x,(s)¢,. Here s is regarded as a coordinate transformation leaving 
the Hamiltonian invariant and 7,’ is the function constructed in the rotated 
coordinate system by the same process as ¢, is constructed in the original system. 


Now if we take b,=Xx,(C™) bn 
wehave ——x,(Ciby =Bx,(C)xe( Cm = D(C" 


=x XAC")h n—1 — x Xol Oi) ra a by 


where 7,’ is constructed in the frame where ¢,, is replaced by ¢,,_), i.e. rotated 


by C’; x, therefore transforms under all operations C’ in the way required of a 


basis function for the 7th representation. 
The basis functions for the completely reduced representation are thus 


5 


1 
Yr, Li Fel tt Paget ea) 
1 
Yr, Ya 76 Pin te tee Past on Ps) 
1 
Ps bs = /6 Co ac ws ts why of why st wds a w° be) 


1 
Dy x= 6 (f1 + why + wids + w8d, + wh; + whe) 


ee ltt ota + 0}, +h,+ 0%; + of, 


1 
Vs pe= 1/6 (fp) + w4hy + ws + db, + wid; + wdg) | 


In the absence of the field the representations [,, [, and T's, [', are equivalent 
to two pairs of degenerate representations but this degeneracy is removed by the 
field (I, Appendix I). 

Using the fact that the only non-vanishing matrix elements of the Hamiltonian 
are those between vectors whose direct products provide a basis for a 
representation containing the totally symmetric I',, the usual secular equation 


of a m.o. treatment based on yw, .. . . %s becomes 
a+22(B')—E 
a+22(—p')-E 


“2+22(wB')—E 
a+22(wp'*)—E 
4+ 22 (w23’)—E 
a+ 2FR(w*s'*)—E 


i 


; 
! 
| 


Diamagnetic Anisotropy of Large Aromatic Systems: III 923 


all off-diagonal elements vanishing. Here (see I) B’ =e” where 6 is related to 
the field strength H and to S, one sixth of the area of a benzene ring}, by 
6 =2n(el/hc)HS. 

We therefore have six distinct energy levels: putting x=(H-—«)/B and 

expanding as far as terms in 6? we find 

Py 4 =2—77 

TP, x=-24+6 

TP, x=1—4/30-—16° 

Dr xal)4/30— i? 

DT; x=—-1-1/30+40 

De x=—14+/30+40 
the equations being labelled by the representations to which the associated 
wave functions belong. The occupied states belong to T,, I's, ['y, and since the 
total z-electron susceptibility is proportional to the coefficient of 6? summed 
over the occupied states it follows that, in appropriate units, yp = —2. 

For comparison purposes it is then sufficient to evaluate Xx, for the various 
molecules from the occupied energy levels, x =x)+%,0+x,6. It will be noticed 
that those levels which are degenerate in the absence of the field give equal 
contributions to the susceptibility—in accordance with I, Appendix I. 


§2. APPLICATION TO-CORONENE AND OTHER STRUCTURES 

The largest molecule to be considered (Molecule 7, Figure 1) contains 
54 carbon atoms whose 2p orbitals ¢,,....¢5, provide the basis for a representation 
of dimensionality 54: it is, however, not difficult to pick out a set of functions 
adapted to a complete decomposition into 54 uni-dimensional representations. 
For every atom belongs to a group of six equivalent atoms in the molecule (in 
Figure 2 the ‘first’ atom in each of the groups, A, B, C,... . is labelled) and 
if the 54 required m.o’s are divided into corresponding sets, based on A-atoms, 
B-atoms, etc. then the representation will be the direct sum of nine six-dimensional 
representations: if now the six functions within each set are constructed from 
the available ¢’s exactly as in (2), then the representation based on the resulting 
54 functions will be completely reduced. Grouping together functions belonging 
to the same irreducible representation, the secular determinant will be reduced 
into a product of six determinants of ninth order involving functions A,, B,... I, 
where C, is conveniently used to denote that function belonging to set C and to 
irreducible representation I’,. The typical ninth order determinant is found to be 


—x+22(yk) 1 
1 nL pe Rr 
Ri —x yk k* 
Re y*h x k 
Nx, v) = k —*x | acl =(3) 

Re =e ke R*8 

ae i. oe 
Rr aye yk? 
ks sel saa —Xx 


+ cf. II. Sis here re-defined so as to give a numerically simpler secular equation. 
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where k=e**., From (3) six determinants are derived by putting, in turn, y=1, 
y=—-1, y=, y=w*, y=w*, y=w**, Here the functions are arranged in order 
Ay, By,.... 1, As, Ba, ....1,,.... (see Figure 2) and the corresponding, secatam 
determinants for a whole series of molecules may be obtained by deleting 
various rows and columns. 

The structures we shall consider are shown in Figure 1: of these only benzene 
and coronene are known molecules but, in principle, there appears to be no reason 
for the non-existence of several of the others. On the other hand structures such 


(8) (9) 
Figure 1. Molecules 1 to 9. 


Figure 2. One-sixth of Molecule 7, showing labelling of typical atoms. 


as Molecule 2 are admittedly hypothetical. This is, however, of little consequence 
since it is our purpose at this stage to examine the dependence of magnetic 
anisotropy on molecular characteristics in general, rather than to effect a detailed 
comparison of experiment and theory. 

Algebraic expansion of the various determinants, up to the ninth order, is not 
difficult but the subsequent calculation of Xx, becomes tedious for the larger 
molecules and a mew technique (see §3) is desirable. Here, however, we shall 
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illustrate the direct treatment by an application to coronene. This molecule has 
been considered by Squire (1938)* and has been the subject of some controversy 
(Squire 1938, Rogers 1947); Squire’s calculations, however, contain a number 
of errors. 

(i) Coronene 


The secular determinant, derived from (3) by retaining elements in the 
first four rows and columns only, is easily expanded giving 


A, (x,y) =(—*x+22(yk))[ — (x? — 3) + 2B(yk?)] —(x®—-1). «0.00. (4) 
Putting in turn, y=1, —1, w, w*, w?, w**, and neglecting terms in 03, 04...., 
this gives the six secular equations 
TD, A,(x,1)  =(«*=2x8 —4x?+ 4x45) +(x? + 46x —100)0? =0 
LT, A,(x, — 1) =(x*+ 2x3 — 4x? —4x% +5) +(— x3 —46x— 100)6? =0 
TD; A,(%,w) =(xt—x° —4x%+-2x+4 2) 
+ 4/3(x8 + 4x — 8) + 4(x? + 46x — 8)02 =0 
T, A,(%,@*) =(x*—x*®—4x7+2x+42) 
+ 4/3(— x8 —4x — 8)0 + 4( — x8 — 46x — 8)62 =0 
T; A,(x,@?) =(x4+ x8 —4x? —2x +2) 
+ 4/3(x3 + 4x + 8)6 + 4( — x? — 46x — 8)6? =0 
Dg Ag(%, w**) = (x4 + 98 —4x? — 2x +2) 
+ 4/3(— «3 —4x%4+ 8)0 + 3(x3 + 46x — 8)6? =0. 
Since the solutions of T';, Ty, and likewise of I';, Ig, yield equal contributions 


to Xx, (1, Appendix I), it is sufficient to consider the four equations Tj, T', Ts, T;. 
Each equation is of the form 


A(x, y) =Ao(a, y) +fi(%¥)0 + fal%, VP =0 wae, (5) 
the suffix zero indicating the zero-field secular equation whose roots (x9) are 
supposed known. Neglecting terms in 6°, 64, .... approximate solutions of the 


form *=%)+%,0+4%,6? may be found by making Taylor expansions of A)(x, y), 
Six, ), fo(%,y) (taking x=x,+h with h=x,0+«x,0") and equating to zero the 
coefficients of @ and 0? in (5). 

Calculating the x,’s in this way for all the occupied energy levels and summing 
we find, for coronene, Xx,= —19-592. Consequently Ay/Ayp =9-:795;. 


(11) Other Structures 
This method may be applied to each molecule: the results of our calculations 
are summarized in the Table. So far the non-orthogonality of adjacent atomic 
orbitals has been neglected but the results are easily modified (Brooks 1941) to 
admit an overlap integral 5=0-25 (the value usually assumed): the modified 
results are also given in the Table. 


Diamagnetic Anisotropies Ay/Ay,, 


Molecule 1 2 3 ot 5 6 dl 8 9 

No overlap 1:0 0:04055 0:589, 9-795, 0-771, 14:94, 35-03, 5:266) 17°85, 
Overlap=0°:25 ) 1-00-02935 0-604, 10°13, 0:827, 15°36, 36:65, 5-307, 17-81, 
Pauling theory Oy ta) 1:0 13:87 14:57 29:20 68:98 34:86 65-29 


* The method of calculation is not fully described. 
PROC. PHYS. SOC. LXIV, IO—A 60 
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§3. AN ALTERNATIVE TREATMENT OF THE REDUCED SECULAR 
EQUATIONS 


Whilst the above method of calculating Xx, by summation of one-electron 
contributions is quite simple in the case of small molecules, particularly when 
the unperturbed energy levels are known, it becomes extremely tedious when 
applied to very large molecules. We shall now show how Xx, may be obtained 
directly by contour integration: in this way the diamagnetic anisotropy may be 
calculated, without any knowledge of the energy levels or wave functions of the 
unperturbed molecule. 

Consider the nth degree secular equation (in general with 60) f(x) =0, 
with n real and distinct roots (x=w’, x’’,.... x’). In the m.o. approximation 
the total energy of m electrons doubly occupying the negative energy levels is 
(in units 8 above the zero level nx) 2X. where the summation is over the occupied 
levels and may be obtained from the identity (Coulson 1940) 


= 5 =|. {oe - a} de. Seale (6) 


Here the contour consists of the y axis and an infinite semicircle to the right* 
of this axis. The inclusion of —n in the integrand conveniently ensures (at least 
in the case 0=0) the vanishing of the integral over the circular contour. The 
substitution z=zy then gives 


ix=— xl. ane =n} dys i) Se ieee (7) 


If f(x) is an even function the integrand is real and numerical integration is 
straightforward; otherwise rationalization is necessary. Now in our case 
Lx = Ux) + GLx,+ FLUx, and we require Ux, only: clearly this quantity may be 
obtained directly as the integrated coefficient of 6? in the expansion of the 
integrand in powers of @. 

Since our equations are not, individually, even functions af x it would at first 
seem necessary both to rationalize and to expand in powers of 0, thereby obtaining 
a very complicated integrand. It appears, however, that our six equations may 
be paired and the pairs, 


(ae z 1) Sie Ge 6 ay Ge = Tae oe (8) 


may each be written in the form 


y= {Aale)+ Ar(2)} + {B,(0) + Ble)}0-+ (CoG) + CY)}M=0 1g 
A_ = {A,(«) — Ap(x)} + {B,(x) — Ba(x)}0 + {Cy(x) — Gone Olea 


where A,, B,, Cy and A,, By, C, are respectively odd and even functions of x. 
If now we take f(x) =A_,A_a relatively simple even function of x is obtained and 
the integration outlined above gives the summed coefficients Uw,, arising from 
the roots of both A, and A_. Also, as we have already seen, the last two pairs in (8) 
must give equal eopuribunaee ux. It remains therefore, to evaluate two infinite 
integrals for each molecule. 


* We require Xw over the positive roots since our is a negative quantity. 
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We find, on putting 
— A=A2-AZ, B=2(A,B,—-A,B,), C=2(A,C,—A,C,)+(B2—-B,2), 
(so that f(x) =A + BO + C6?) 


Lei AGAC AAB* — BBA? 
Lk — = {(=S"5) af (“3 — }} dy ee eeee (10) 

where A=(xdA/0x), B=(xdB/Ox), C=(xdC/dx) and x=zy throughout. 

In the cases with which we deal the coefficient of 0? in the integrand of (6) vanishes 
at least as strongly as -s-?—-hence the absence of any contribution from the 
circular contour. The simple polynomials A, B, C, A, B, C, contain up to nine 
terms in the case of Molecule 7, so the integrand of (10) is best evaluated 
numerically without further algebraic expansion. It should be noted that the 
integral associated with the first pair of equations in (8) is relatively simple, the 
second term in (10) vanishing. Also each of the two integrals required for a 
given molecule is determined mainly by the region 0<y<4: an asymptotic 
expression is easily derived for the integral from 4 to oo. If the roots xy are in 
the first place unknown the determination of Xx, by this method may involve 


considerably less labour than proceeding by way of the numerical solution of 
the secular equations. 


The Inclusion of Overlap 
The contour integral method of obtaining the total 7-electron energy can be 
modified to take overlap into account. The treatment given by de Heer (1950) 
may be applied directly to magnetic problems, for, with overlap 5 and 
replacement of 8 by y=6 —6a, the modified total energy is found using 


ae al. ; — ae =n} AG TOS 7’) 


instead of (7). Now in obtaining Xx, with neglect of overlap, we have simply 
expanded the quantity in brackets in powers of @ and selected the 6? term as a 
measure of the 7-electron susceptibility; since the factor (1 +6*y?)-1, introducing 
the overlap correction, is common to every term in the expansion, it is clear 
that the modified Xx, may be calculated directly from (10) provided the integrand 
is weighted with the factor (1+6°y")-1. The smallness of 5% suggests immediately 
that overlap should have little effect on our calculations; it appears moreover, 
as de Heer has shown for the total energies, that the modification is roughly 
equivalent to multiplication by a factor which differs little from molecule to 
molecule. The ratio Ay/Ayp, might therefore be expected to show very little 
dependence on overlap. 

This result is a quite general one provided we use simple m.o. theory and 
deal only with pure hydrocarbons; it therefore explains immediately previous 
observations that the calculated anisotropy is not strongly dependent on overlap 
(II, Brooks 1941). Whenever the secular determinant, involving a perturbation 
parameter of any kind, can be made formally similar to that obtaining without 
overlap by the change of variables x—x/(1+6x), B+ 8—d« then the total 
perturbation energy, in any order, is only weakly dependent on overlap. ‘The 
conditions under which this reduction is possible have been examined by 
Chirgwin and Coulson (1950); it is important to notice that a similar result 
cannot be established for molecules containing hetero-atoms. 


60-2 
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§4. DISCUSSION 
In the Table our results are accompanied by the values of Ay/Ayp calculated 
using the semi-classical method of Pauling (1936). Both types of calculation 
predict a steady rise of Ay/Ayp for Molecules 1, 4, 6, 7, although quantitatively 
there is considerable divergence; but startling contradictions appear elsewhere. 
It is evident that the Pauling method is, in general, likely to be misleading: in 
seeking an explanation we shall examine separately the different kinds of molecule. 


Molecules 1, 4, 6, 7 

The only structures considered hitherto have been of this type* (II, London 
1937); they are completely condensed cyclic hydrocarbons in which every carbon 
atom is included in at least one six-membered ring and there are no ‘holes’ in 
the structure (cf. Molecules 8, 9). Although for the smaller molecules considered 
in these papers and elsewhere (up to, say, coronene) the Pauling theory has shown 
itself of value as a crude substitute for the more complete theory, our results 
show clearly that, even in the most favourable circumstances, it must be 
abandoned in dealing with very large molecules. On the other hand the results 
of m.o, theory appear to be entirely plausible and our contention that the 
computed values of Ay/Ay, should be more or less independent of overlap is 
well substantiated. Indeed, this near-independence is observed with every 
molecule considered, with one exception: for Molecule 2 the integral (7’) has 
an abnormally low value, arising from cancellation of negative and positive areas 
—in this case the whole integral is not insensitive to changes in 6. 


Molecules 2, 3, 5 

The Pauling treatment is evidently quite inapplicable even to small molecules 
when these possess ‘stray’ C—C bonds (i.e. bonds not contained ina ring). It is, 
however, still desirable to seek some physical interpretation of our results. 

The sudden drop in susceptibility in going from Molecule 1 to Molecule 2 
would be explicable immediately if there were considerable ‘bond fixation’ 
but calculations of resonance energy and bond order preclude this possibility: 
in the language of valence-bond theory this means that, even in the ground state, 
the single ‘fixed-bond’ structure (all stray bonds double) must be supplemented 
by important excited structures. ‘The extremely small z-electron susceptibility 
cannot therefore be attributed to bond fixation, but the very fact that the only 
possible unexcited structures show the stray bonds double also places a restriction 
upon the types of excited structures which can be drawn; and the existence of an 
appreciable susceptibility is linked with the occurrence of certain particular 
types. ‘Thus Brooks (1940) has shown that an appreciable susceptibility can 
occur only when there exists a group of excited structures by means of which 
an electron can be ‘ passed round the ring’. 

Now in the case in point there are clearly no structures of this type, for no 
ring-atom can possess an unbonded electron in any valence bond structure; 
if the atom is not shown bonded to either of its neighbours in the ring it must be 
bonded to the adjacent stray carbon atom which can always provide an electron 
for such a bond. It should be noticed that excited states showing unbonded 
those atoms not contained in the ring do not contribute very much to the 
susceptibility because (see Brooks 1940) the relevant resonance integrals then 
involve non-adjacent atoms. The extremely small value of Ay/Ay, predicted 
for Molecule 2 by m.o. theory thus finds support in qualitative valence bond 


* With one exception : Evans, de Heer and Gergely (1949). 
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considerations : this is a very pleasing result in view of the impossibility of dealing 
with the vast numbers of valence bond structures involved in the larger molecules. 

This interpretation also explains a similar effect noted by other writers (Evans,, 
de Heer and Gergely 1949) in an application of m.o. theory to p-benzoquino- 


_dimethane : in this case there are two carbon atoms within the benzene ring 
_ which cannot show free electrons in any valence-bond structure and therefore 


serve as ‘barriers’ tending to prevent the passage of electrons round the ring. 
Passing now to Molecule 3, it is evident that structures in which the peripheral! 


_ atoms are bonded are no longer incompatible with the presence of free electrons. 
on ring-atoms; indeed, if structures of this kind only were involved we should 
_ expect the system to show little difference, magnetically, from the simple benzene 


ring. ‘The considerable rise in Ay/Ay, is thus to be anticipated. 
The addition of stray carbon atoms to coronene might be expected to annul 


the effectiveness of the outermost complete ring of carbon atoms, exactly as in 
the case of Molecule 2, so that Ay for Molecule 5 should have a value only of the 


order Ayp. Our results support this conclusion and it seems likely that the effect 
is quite a general one: bordering a condensed structure with stray carbon atoms 


reduces the magnetically effective size of the system. 


The above results may be significant in discussing the effect of different 
kinds of ‘edge’ upon the magnetic properties of carbon particles. ‘Taking as a 
“standard edge’ that in which the free carbon valencies are satisfied by hydrogen 


atoms, we suggest that (i) replacement of hydrogens, —H, by — OK leads to a 


marked reduction in the magnetically effective size of a carbon particle, roughly 
equivalent to striking out all those rings with one or more such groups attached 
(cf. Molecules 2, 5); (ii) replacement of hydrogens, —H, by ‘actylenic tails’, 
—C=CH, has little effect on the susceptibility of the original particle. This is 
suggested by the fact that the secular determinant for such a system, 


e.g. benzene with six ‘tails’, is formally identical with that for a system of a 


type already considered, Molecule 2. 


Molecules 8, 9 

The removal of atoms from the centre of a condensed system evidently has a 
very marked effect on Ay/Ay,; moreover, the effect may be entirely the opposite 
of that predicted on semi-classical grounds*. Thus, on Pauling’s theory, the 
magnetic anisotropy of ‘ring molecules’ such as 8 and 9 should depend jointly 
on the area of the whole ring (roughly that of the molecule) and on the 
‘conductivity’ of the ring, the latter presumably being larger for Molecule 9 
than for Molecule 8 because of the increased number of ‘conducting’ bonds. 
The reason for the failure of this theory is primarily the assumption of independent 
ring currents ; the necessary interdependence of these currents is easily illustrated 
by valence bond considerations. 

Unfortunately a thorough valence bond treatment of Molecules 8 and 9 is 
precluded by the vast number of available structures; we shall merely make the 
customary assumption that in general the higher the order of an excited structure 
the less its contribution to the ground state wave function. Consider now 
Molecule 8: by far the most abundant of the less highly excited structures are 
those in which two opposing atoms in each ring possess unpaired electrons. This 


* On Pauling’s theory the susceptibility/carbon atom for Molecule 9 is not less than that for 
Molecule 7 but greater. 
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group of structures makes it possible for electrons to ‘hop’ round the individual 
benzene rings, but it cannot give a net flow between the rings. In classical language 
the two conducting paths (inner and outer) which encircle the whole molecule 
must carry equal and opposite currents: on the other hand the Pauling theory 
predicts a large current along each path in the same direction. If only the excited 
structures above were involved in the ground state wave function we might 
therefore expect a value of Ay/Ay, of the order 6-0, corresponding to that for six 
unconnected benzene rings. Actually other structures, in which the connecting 
bonds are shown double, must be involved if there is appreciable conjugation 
between the rings, but it is not difficult to see that these also cannot be responsible 
for electron-hopping around the entire periphery of the molecule, and their 
presence therefore serves only to reduce the weight of those structures giving 
benzene-type ring currents. The value Ay/Ay, =5:°3, given in the Table conforms 
quite well with these conclusions. 

In Molecule 9 the restrictions on the types of excited structure possible are 
considerably less severe: there is now a much greater variety of possible 
electron-hopping paths encircling more than one benzene ring. We may 
therefore expect Ay to be considerably greater than for twelve independent 
rings although less than that for the whole molecule (7) in which the restrictions 
are even further relaxed. The value Ay/Ay, =17-8 thus appears quite plausible, 
whereas the semi-classical result (65-3) even fails to suggest the marked difference 
between Molecules 7 and 9. . 
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APPENDIX 


Parts I and II. Errata 


Several errors had arisen in I and II (McWeeny 1951) during preparation 
of the MSS. They are: 

The exponent should be reversed in sign in equation (4), p. 265 and in 
line 9, p. 267; it should be halved in equation (6), p. 265: it should be halved 
and reversed in sign in line 8, p. 267. 

Last line, p. 274, should read: U=—4,HYx'+}3HXy'=4R.1r'xH=—-R.A’, 
and in line 4, p. 275, exp{—27i«xR.A} should read exp{27aR.A’}. 
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Paramagnetism of Cerium Ethylsulphate at Low 
Temperatures 


Single crystals of cerium ethylsulphate, Ce(C,H;SO,);.9H,O, have been investigated 
at liquid helium temperatures by measurements of susceptibility and by paramagnetic 
resonance experiments. This substance shows no paramagnetic resonance at temperatures 
above 20° k., the spin-lattice relaxation time being so short that the absorption lines are 
broadened to such an extent as to be undetectable (Bleaney and Scovil 1951). According 
to the theory of Elliott and Stevens (unpublished) the ?F;,. ground state of cerium ion is 
split by the crystalline field in this salt into three doublets with J,= +1/2, + 3/2, + 5/2, the 
2 axis being the hexagonal axis of symmetry of the crystal. The +3/2 level is some 
130 cm above the others, so that its population is negligible at temperatures below 
20°. Our experiments have enabled us to find the magnetic properties and energy 
separation of the other two levels. 

The resonance experiments were made in the usual way, the absorption in the specimen 
being measured by its effect on the damping of a cavity resonator and displayed on a 
cathode-ray tube as a function of the applied magnetic field. The resonant wavelength 
was 3:25 cm. Experiments were made on a dilute salt (a mixed crystal of one part cerium 
ethylsulphate and 200 parts of lanthanum ethylsulphate) and on a concentrated salt. The 
liquid helium used was obtained from a large Simon expansion liquefier with transfer 
siphon, which has just been brought into use in this laboratory. 

At 4-2° kK. two resonance lines were found in the diluted salt, one due to the + 5/2 level 
and the other to the +1/2 level. The results may be expressed in terms of the following 
spectroscopic splitting factors, or g-values (defined by hv=gfH, where h= Planck’s constant, 
8=Bohr magneton, v=resonant frequency in a magnetic field H); g, and g, are the g-values 
in a magnetic field respectively parallel and perpendicular to the crystal axis : 


gy(5/2)=3-704 0-02, g,(5/2)=0°3540-15, g4(1/2)=1-02+0-02, g,(1/2)=2-17+0-02. 


At an angle @ to the axis the g-value is given by g?=g, cos? 6+ g,? sin? 6. 

At 4:2°x. the J,= +1/2 line was about ten times as intense as the J,= + 5/2 line, both 
measured near the angle at which their g-values are equal. From the change in their 
relative intensities as the temperature was lowered to 2:5°K. it is concluded that the 
+ 1/2 level lies lower by 3+1cm-. An interesting feature of the results is that in this salt 
we are able to observe absorption due to an excited state. Moreover, this absorption arises 
from a transition in which AJ,=5. 

In the concentrated salt no resonance was detected at 4:2°k. except for a line from 
neodymium present as an impurity. The width of this line was markedly dependent on the 
orientation of the crystal, being narrowest with the magnetic field perpendicular to the 
crystal axis, and three times broader at 15° to this direction. In accordance with recent 
work by Bleaney, Elliott and Scovil (1951), the broadening of the neodymium line is 
proportional to the g-value of the broadening agent, which is therefore cerium in the 
+5/2 state. On cooling to 2:5° x. a very broad line (line width approximately 700 gauss) 
appeared with g-values depending on angle in a similar way to that of the +5/2 state in 
the dilute salt. It is concluded that in the concentrated salt the +5/2 state lies lowest, 
contrary to the case of the dilute salt. 

This conclusion is confirmed by the susceptibility measurements, which were made 
on four single crystals mounted together with their hexagonal axes parallel. Measure- 
ments were made with this axis (a) parallel and (b) perpendicular to the measuring field, 
using the apparatus and method described elsewhere (Benzie and Cooke 1950). Only 
the temperature-dependent paramagnetism is measured. ‘The results indicate that in 
the temperature range 1° K. to 20° k. there are two levels contributing to the susceptibility, 
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and that the +5/2 level is the lower. In terms of g-values the susceptibility y of such 
a system should be, to a first approximation, 


ies NB? 85/27 +81/27e AIT 
ART 1+e-A/f 


where RA is the energy separation of the two states. Good agreemént with experiment is. 
found by taking g, (5/2) =3-80+ 0-04, 2, (5/2) =0+ 0-4, 9,,(1/2)=1-0+ 0-2, 9,(1/2) =2°25+ 0:2, 
ang N= Dest Danke 

Previous measurements of the susceptibility have been made at temperatures down 
to 14°x. by Fereday and Wiersma (1935) and to temperatures down to 1:37°K., in the 
parallel direction only, by Becquerel et al. (1938), who used the Faraday effect and assumed. 
that the rotation of polarized light was proportional to the susceptibility, as verified by 
measurements above 14°K. The results of Becquerel etal. in the parallel direction are in 
good agreement with ours. In order to fit the results of Fereday and Wiersma in the 
perpendicular direction a temperature-independent term must be added. This would 
be measured in their experiments but not in ours. The presence of such a term may be 
accounted for theoretically (Van Vleck 1928) ; it is due to off-diagonal elements in the 
magnetic moment matrix connecting the J,= + 3/2 state with the ground states. 

The final expressions for the susceptibility in E.M.U. per gm.mole, valid up to 
20° K., are’: 

0:094 14:44+¢-75/2 0:094 5:0e-7 5/7 
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At higher temperatures the effect of the J,=+3/2 level must be taken into account, and 
results at 60° k. already lie off the curve. 

The results show that this substance is of particular interest for experiments below 
1° k. by the adiabatic demagnetization method. At such temperatures only the + 5/2 state 
is occupied, and our results show that the susceptibility is at least a hundred times greater 
along the crystal axis than at right angles to it. We propose, therefore, to make experiments 
by magnetizing a single crystal in a large magnetic field at 1° k., and then to cool it by 
adiabatic demagnetization. The magnetic axis of the crystal, the magnetic field used for 
cooling, and the coils used for measuring the susceptibility will all be parallel (vertical). 
After such a demagnetization a horizontal magnetic field may be applied without much 
effect on the temperature of the specimen. In effect the specimen has been cooled by 
rotation of the magnetic field. Cerium ethylsulphate is therefore of value for cooling 
substances on which it is desired to make magnetic experiments at temperatures below 
1° x., for example the investigation of the phase diagram of a superconductor, or nuclear 
resonance experiments. 

We wish to thank Mr. H. E. D. Scovil for supplying the crystals and much helpful 
advice; also Mr. A. J. Croft and Mr. G. R. Hercus, who provided the liquid helium. 


The Clarendon Laboratory, G. S. BoG Le. 
Oxford. A. H. Cooke. 
25th June 1951. S. WHITLEY. 
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The Paramagnetism of Cerium Ethylsulphate: Theory 


The theory of the resonance phenomena in rare earth ethylsulphates as reported in @ 
previous letter (Elliott and Stevens 1951) has now been taken to second order in 
perturbation theory and applied to cerium ethylsulphate. The experimental results on 
this salt have been reported by Bogle, Cooke and Whitley (1951) who have obtained 
accurate g-values, particularly for the dilute salt. We find that their results can be explained 


Letters to the Editor 933 


t 
i 


| by using the following values for the crystal field parameters, as defined in our earlier letter : 
=|elA,%”?=0, —|e|Ay°r4=—70 cm-}, —|e|A,'r4=8 cm7}, —|e|A,°r®=—130 cm-, 
—|elA’ry*=1300 cm-!. The Cayeesitine potential used is slightly different from that ee 
previously in that we have included a term of the form A,°2(x*—3xy*) to account for 
departures from D3, symmetry. As can be seen from the value given above, this extra 
term is small, but it proves to have an important effect in this salt because it so happens 
that the levels corresponding to J,=+4 and J,=+5/2 are near together and a term of 
this type has matrix elements epupbee the two fevele! The second order calculation has 

been made by assuming that only the important contributions come from the other 
multiplet (J =7/2) which we have taken to be 2500 cm~! above the ground multiplet. 
In the following Table we compare the experimental and theoretical results. 


J,(approx.) Theoretical g Experimental g Theor.energy Exp. energy 
Sil 8&1 Sil 81 (cm~) (cm~) 
+4 1-02 217, 1:02+0:02 2-173 0-02 0 0 
+§ 3:70 0-42 3:70+0:02 0:35+0-15 3 Sap i 
+8 DS 0 — — 160 —_ 


The values of the parameters are in reasonable agreement with those obtained from a 
model of the ion and its surroundings but they cannot be regarded as known with any great 
certainty as it is possible to use appreciably different values and still remain within the 
experimental errors. If the value of g, (+ 5/2) or the properties of the + 3/2 level were 
accurately known the crystal field values could also be more closely determined. 

In the concentrated salt the uncertainties in the g-values are much greater and there 
is an even wider range of possible values for the field parameters. Thus at present there 
is little to be gained by applying the theory in detail to the concentrated salt. 

A full account of this work is being prepared and will be published shortly. 


Clarendon Laboratory, RoJ. Buurorr. 


Oxford. K. W. H. STEVENS. 
4th July 1951. 
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Dipole-Dipole Interaction in the Rare-Earth Ethylsulphates 


Paramagnetic resonance spectra of several rare earth ethylsulphates show interesting 
features in line shape which can be related to spin-spin and spin-lattice effects in the 
crystal. The former gives usually a broadening of the resonance lines which is due to 
the magnetic fields in neighbouring ions, randomly superimposed. ‘This broadening 
consists really of a closely spaced fine structure, some detail of which can be resolved if 
the crystal is diluted with an isomorphous diamagnetic salt (Bleaney and Ingram 1950). 
The ethylsulphates are, however, sufficiently dilute themselves to give some resolution. 

The spectrum of neodymium ethylsulphate at 20° kK. is shown in Figure 1. With the 
magnetic field parallel to the hexagonal axis a triplet structure is visible ; this is not due to 
hyperfine structure, which gives lines of only 14% relative intensity. It is well accounted 
for by dipole-dipole interaction with the next neodymium ions, which lie at a distance 
c=7a. away along the hexagonal axis. A simple classical picture shows that the resultant 
field of these two neighbours is + 2 g,6/c*, or 0, the latter being twice as probable as either 
of the former. Thus we expect a triplet line, with a spacing of about 50 g,=180 gauss, 
the satellites having half the intensity of the centre line. The rotating fields of the neighbours 
give a resonance interaction which is small owing to the low value of g, ; this should slightly 
increase the splitting. With a field perpendicular to the axis, the classical splitting should 
be 25 g,=50 gauss and the resonance displacement is more important; the triplet is not 
resolved though the shape indicates the presence of the satellite lines. Broadening by the 


more distant neighbours gives a root mean square width (Van Vleck 1948, Pryce and 
Stevens 1950) of approximately 45 gauss for either direction, which is consistent with 
experiment. 

The resonance in a concentrated gadolinium cthylsulphate crystal was observed as a 
wide line some thousands of gauss broad. The elementary picture here is complicated 
by the fact that the neighbouring gadolinium ions are now dipoles with eight possible 
orientations instead of two. Thus the nearest neighbours alone will give a large number of 
lines stretching out to +14g¢/c?=+700 gauss from the centre of each line of the 
gadolinium fine structure. The seven lines of the latter are therefore spread out into one 
broad line. 

The effect of spin-spin broadening on the spectrum of an impurity was observed in a 
crystal of neodymium ethylsulphate containing about 0:2°% of the gadolinium salt. At 
90° x. the neodymium spectrum is invisible owing to the very short spin-lattice relaxation 
time, and for the same reason the Gd spectrum (Figure 2 top) shows only the width normal 
for high dilution ; this width is due to the nuclear ‘moments of the protons in the water of 
crystallization. At 20° x. the relaxation time of Nd is so long (>10~® seconds) that the 
full effect of its dipole field is observed in the Gd spectrum (Figure 2 centre), each line 
being split into three. With dissimilar ions there are no resonance effects and the classical 
picture is a better approximation than in the concentrated Nd salt. Thus the intervals 
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Figure 1. Triplet structure of ggg TE 


the neodymium line in 
neodymium ethylsulphate : 
Top, magnetic field parallel 
to crystal axis ; bottom, 
magnetic field perpendi- 
cular to axis. 


Figure 2. Gadolinium impurity spectrum in neodymium ethylsulphate, 
Top, at 90° k.; centre, at 20°k.; bottom, theoretical spectrum at 
ZR 


of the triplet should be 180 gauss when the spectrum is observed parallel to the axis. 
As the fine structure intervals in Gd are about 300 gauss, some of the satellite lines overlap 
so that 13 lines are observed in all. The root mean square width should be about 40 gauss 
which is of the order observed experimentally. 

Similar experiments have been made with gadolinium impurity in cerium and samarium. 
No transitions due to either of these two ions have been observed down to 14° k. in the 
ethylsulphates. In the samarium salt the Gd impurity spectrum appeared to have the 
normal lire width, but in the cerium salt a small broadening of the lines parallel to the 
axis was visible at 20° k., which increased considerably on pumping down to 14° x. At the 
same time saturation effects, which are usually present in gadolinium at these temperatures, 
could not be obtained. 

‘The absence of magnetic resonance lines due to cerium and samarium can be accounted 
for in two ways : (a) there may be no allowed transitions in the Kramers doublets occupied 
at these temperatures (this occurs when one of the principal values of the g-tensor is zero) ; 
(6) the relaxation time may be so short that the lines are too broad and too weak to detect. 
In the former case the gadolinium impurity lines should be broadened in directions where 
the g-values of the Kramers doublets of the predominant ion are not zero. In the latter 
case no broadening will be observed until the temperature is lowered so far that the 
relaxation time of the predominant ion approaches the Larmor precession period of the Gd. 
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‘The results indicate that (b) is the correct explanation for both samarium and cerium. 
The absence of saturation effects is due to the shortening of the relaxation time of the 
Gd impurity through transitions induced by the fluctuating magnetic fields of the rapidly 
relaxing cerium and samarium ions. In the case of cerium the correctness of this 
explanation has been confirmed by Bogle, Cooke and Whitley (1951), who found 
paramagnetic resonance lines at helium temperatures. 


Clarendon Laboratory, __B. BLEANEY. 
Oxford. Reet ig IDE oni ane 
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A Geiger Counter Spectrometer for Soft X-Rays 


The spectroscopy of soft x-rays, of wavelengths which may be as great as several 
hundred Angstrém units, has been made possible by the use of ruled gratings in grazing 
incidence settings. ‘The x-rays are detected photographically. Wavelengths can be 
determined with precision, but intensity measurements are of more doubtful accuracy. 
‘They usually involve the assumption that the reciprocity law is valid and that the output 
of the x-ray tube is constant. -Chalklin (1948) has used a method which avoided these 
assumptions, but the uncertainty of the photographic photometry was considerable. Yet 
intensity measurements are needed to give the precise shapes of the soft x-ray bands that 
involve the valence electrons and which have afforded striking evidence on the electron 
theory of solids; and intensity measurements are needed in the study of the absorption 
of soft x-rays by gases and solids. We have therefore attempted to construct an instrument 
which, with the photographic plate replaced by a Geiger counter, would make direct 
measurements of intensity without sacrificing resolution. ‘The purpose of the present 
Letter is to indicate that this attempt has eventually proved successful. In undertaking 
the task we were encouraged by the success of the early unfocused counter spectrometers 
of Neufeldt (1931) and Messner (1933). 

Owing to the high absorbability of the soft x-rays the spectrometer is enclosed in a 
vacuum box. The optical system consists ofa fine slit as effective source, a one-metre 
Siegbahn concave grating and a ‘detector’ slit, all mounted on the circle of focus of the 
grating. The grazing angle of incidence is 2:5°. The detector slit, with the Geiger counter 
mounted behind it, is moved along the circle of focus, and its position measured, by 
operations conducted outside the vacuum box. To permit the entry of the soft x-rays the 
counter is provided with a very thin window consisting of a double celluloid film of total 
thickness 1,000. 

Tests on the Geiger counter with gamma-rays show that, although the counting rate 
is proportional to the intensity for rates below 3,000 per minute, there is an increasing 
deviation for higher speeds of counting. Allowance for this has been made in the curve 
of the Figure, which represents the first order of the molybdenum line Myy, y—-Nnu, m- 
The three sets of points in the Figure correspond to different radiation intensities, and for 
comparison the B set have been multiplied by 1:12 and the C set by 0°565. The 
wavelength scale has been obtained from the dispersion curve, and by matching the 
wavelength of symmetry of each set to the value of 64:35 a. given by Siegbahn and 
Magnusson (1934). The agreement of the three sets (if we except the single point at the 
highest counting rate of 10,000 counts per minute) indicates that the counting rate, when 
corrected, does indeed measure intensity. The smoothness of the curve indicates the 
reproducibility and steadiness of the apparatus. 

For the instrument to be fully successful it is necessary that the counting rate be 
great enough to allow accurate readings to be taken at reasonable speed when the slits 
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are sufficiently narrow to yield satisfactory resolution. The slit widths at present in use 
yield a theoretical resolving power of roughly 4A (where A is the wavelength in A.), and, 
with this resolution, it is seen from the Figure that the counting rate is more than adequate 
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Intensity curve for the first order of the Mry, v—Nu, m line of molybdenum 
(uncorrected for instrumental line width). 


for emission line measurements. Since the highest x-ray tube current was 1 ma., it is 
clear that it could be increased sufficiently to permit satisfactory measurements to be made 
of the continuous radiation. Hence it appears that the apparatus should be of value for 
absorption as well as for emission spectroscopy. 


Physics Department, Ie Ib; 
Canterbury University College, F. C. CHALKLIN.. © 
Christchurch, New Zealand. 
29th June 1951. 
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On the Oscillator Strength of the x7II—a?A Transition 
of the CH Radicle 


Stephenson (1951) has recently employed the |.c.a.o. method in calculations on the 
CH, x*JI—a?Atransition, He finds that the oscillator strength f is about 0-008 and states 
that this seems to indicate that the value 0:06 assumed by Dunham (1941) is too high and 
that (considering the nature of the approximation) it is in good agreement with the value 
0:002, given both by the experimental work of Lyddane, Rogers and Roach (1941) (as 
corrected by Dwyer and Oldenberg (1944)) and by the independent theoretical work of 
Herzberg (quoted by Stromgren (1948)). Astrophysicists eagerly await a reliable determin- 
ation of the quantity in question. "They may, unfortunately, be led into thinking that new 
evidence has been provided, and that several different approaches now favour the view that f 
is low as opposed to only a suggestion that it is moderate. Quite apart from the fact that 
0-008 is not significantly less than the geometric mean of 0:06 and 0-002, and from the fact 
that simple calculations on perpendicular-type transitions should be treated with special 
reserve, it must be pointed out that the summary of the literature whichis given is misleading. 

(i) The value 0:06 is not based on a mere assumption, it is based on measurements 
performed by Dunham (1941) in collaboration with King and, moreover, is in concord with an 
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admittedly indirect deduction made by Gaydon and Wolfhard (1949) from a spectroscopic 
study of the low-pressure flames. 

(ii) Lyddane, Rogers and Roach (1941) did no experimental work; like Stephenson, 
they used quantum mechanics. They, however, studied a number of transitions and, being 
doubtful (as well they might be) of the accuracy of the usual approximation to molecular 
wave functions, they attempted to improve their results by arbitrarily multiplying them by a 
correction factor which was chosen so as to bring the computed oscillator strength of OH, 
x?IIT—a?X+ down to the value obtained in the laboratory by Oldenberg and Rieke (1938). 
When the latter value was later revised by Dwyer and Oldenberg (1944) it was Herzberg 
(1950) who made the necessary alteration to the original correction factor and hence estimated 
f for CH, x?IT—a?A to be 0-002; he did not carry out a theoretical investigation. There is 
thus only a single determination, and it is but semi-empirical. 

The present letter seeks merely to prevent the growth of any misconception amongst 
astrophysicists as to the reliability of our knowledge of the oscillator strength of a transition 
of such importance. It must not be inferred that the writer favours a moderate value; on 
the contrary he favours a low value. 


Department of Applied Mathematics, D. R. Bates. 
Queen’s University, Belfast. 
8th August 1951. 


DunuaM, T., 1941, Publ. Amer. Astr. Soc. 10, 123. 

Dwyer, R. J., and OLDENBERG, O., 1944, 7. Chem. Phys., 12, 351. 

Gaypon, A. G., and WoLFuarD, H. G., 1949, Proc. Roy. Soc. A, 199, 89. 

HERZBERG, G., 1950, Spectra of Diatomic Molecules, 2nd edition (New York: Van Nostrand), p. 386. 
Lyppang, R. H., Rocrrs, F. T., and Roacu, F. E., 1941, Phys. Rev., 60, 281. 

OLDENBERG, O., and RIEKE, F. F., 1938, 7. Chem. Phys., 6, 439, 779. 

STEPHENSON, G., 1951, Proc. Phys. Soc. A, 64, 666. 

Srromcren, B., 1948, Astrophys. 7., 108, 242. 


On Certain Applications of the Fermi-Thomas Model 
to the Nucleus 


Recent publications (Yang 1951, Ivanenko and Rodichev 1951) on the so-called ‘magic 
numbers’ of nuclei (Haxel et al. 1949) prompt the following remarks : 

(i) ‘The Fermi~Thomas model by its very nature cannot explain shell formation in its 
energetic aspect. 

(ii) The other main aspect of shell theory is the prediction of nuclear spins. We may 
take as basic assumption that the total angular momentum J of an even—odd nucleus is the 
total angular momentum 7 of the odd nucleon. This assumption does not nccessitate 
j-j coupling, and the Fermi~Thomas model, which, of course, deals with orbital angular 
momenta only, is at any rate not inadmissible on these grounds. 

(iii) If we assume a certain nucleonic density distribution (as, for instance, does Yang 
(1951)) the Fermi model yields the number of nucleons of a given orbital angular 
momentum, in an approximation which is the W.K.B. approximation to the individual 
particle solution for the corresponding potential. Further, it shares with the latter all the 
faults of any individual particle model and, in addition, requires a translation from the 
continuous to quantized angular momenta which is arbitrary within certain limits. 

(iv) The number 7; of nucleons of a given orbital angular momentum / amounts to one 
in the statistical model for a total number of nucleons N, which is not, in general, at all the 
same as N,+1, where N,j is the total number of nucleons corresponding to n being exactly 
zero (e.g. in the case of uniform spatial density the values are N,=21-1 as compared with 
Ng 16:6) for/=3))! 

(v) The first appearance of certain total spins docs not coincide with the higher magic 
numbers; on the shell theory (Mayer 1949) we expect /=1, 2, 3, 4 first to occur for a proton 


938 Letters to the Editor 


or neutron number of 3, 9, 21, 41. Using the shell interpretation, we find these / values 
empirically first in those regions. 
(vi) For uniform density the Fermi model determines the absolute values n, by the 


equation = [8{U+ 1} /a] iy? 1+ 1) 2 U4 1) PH? cos [+ P29, 


where y=(97N/8)!"°, N being the total number of nucleons, with no further adjustable or 
empirical parameter. 

The atomic numbers corresponding to first appearance of successive orbital angular 
momenta calculated from this equation (i.e. for a ‘square distribution’ of nucleonic density) 
are too small by a factor of nearly two. 


Department of Mathematical Physics, H. R. PANETH. 


University of Birmingham. 
2nd August 1951. 
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Angular Correlation in the Reaction °Li(d, p)’Li*(y)’Li 


The branching ratio in the slow neutron reaction 1°B(n, ~)*Li suggests a spin of 5/2 
and odd parity for the 478 key. excited state in 7Li (Hanna and Inglis 1949) ; however, 
recent measurements have given isotropy within the experimental errors for the angular 
correlation between the alpha-particles and gamma-rays in this reaction (Rose and Wilson: 
1950) and for the angular distribution of the gamma-rays in the reaction *Li(p, p’) *Li*(y)*Li 
(Littauer 1950). This indicates a value 1/2 for the spin but is not quite conclusive (Inglis. 
1951), since chance values of the nuclear matrix elements, which cannot be calculated, . 
might give very slight or zero correlation for higher spin. It is unlikely that this will 
happen with all the reactions leading to the excited state, and, therefore, it was thought worth 
while to measure the angular correlation between the protons and gamma-rays in the 
reaction ®Li(d, p)*Li*(y)7Li. 

A thin target of ®Li was bombarded with 630 kv. deuterons from the Cavendish 1 Mev.. 
high tension set, and the protons, detected by an anthracene screen and E.M.I. 
photomultiplier, type 5032, were observed at 90° to the incident beam. Another 
scintillation counter, using a solution of terphenyl in toluene, and capable of rotation in the 
plane perpendicular to the direction of the deuteron beam, was used to detect the 
gamma-rays. Proton—gamma coincidences were obtained with a resolving time of 
15 millimicroseconds ; the circuit was similar to that previously described (Curling and 
Newton 1950), but used secondary emission valves (E.M.1I. type E.2133) as pulse limiters,. 
the undelayed and delayed channels being fed from the anodes and secondary cathodes. 
respectively. The average true-to-chance ratio during the Rexpecument was 2:1 and the 
true counting rate three per second. 

The number of proton—gamma coincidences per proton was determined at eight 
angles @ in the range 90° to 160°. No individual reading took more than two minutes, and 
the various angles were taken in as random an order as possible, so that any long-term 
drift in the apparatus should cancel out. In fact no such drift was observed. The results. 
were corrected for absorption in the target holder, scattering and mechanical imperfections, 
by measuring the angular distribution of the gamma-rays from the _ reaction 
*Li(p, p’)*Li*(y)’Li. ‘This must be isotropic in the plane perpendicular to the proton beam. 

From the known lifetime of the excited state (Elliot and Bell 1949), the average velocity 
of the recoiling nucleus at the instant of emission, and hence the centre-of-mass and 
Doppler corrections were calculated. After all corrections had been applied the result 
was fitted to a curve of the form A+ Bcos?@ by the method of least squares, and the 
correlation coefficient B/A calculated to be —2:7+ 2°5°%. (If no centre-of-mass and Doppler 
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correction is applied the result is —1°8+2:5%.) Thus within the limits of accuracy of 
this experiment there is no evidence for the angular correlation being anything other than 
isotropic, or the spin of *Li* having any other value than 1/2. 

If the spin of ‘Li* is other than 1/2 some correlation is expected since the protons 
are known to have angular momenta up to 3 (Dunbar and Hirst 1951). In particular, if 
it is 5/2 (odd) there should be a correlation of the form 1+C cos?6+D cos*6, most likely 
with C>D, because the magnetic dipole strength is expected to be about 200 times larger 
than that of the electric quadrupole (Inglis 1951). 

Since three experiments of this type all suggest isotropic correlations to quite a good 
order of accuracy, the evidence for the spin of the 478 kev. excited state of ’Li being 1/2 now 
seems quite strong. This value is also in agreement with that predicted by the shell model 
(Mayer 1950). 

I wish to thank Sir John Cockcroft and the Atomic Energy Research Establishment 
for providing the separated °Li target. 
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Cambridge. 
1st August 1951. 
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A Study of Metal Transfer during Sliding, using 
Radioactivation Analysis 


A number of methods have been used for detecting the fragments transferred when 
one metal slides on another, and a comprehensive review has been given by Burwell (1947). 
One of the most common methods involves the sliding of a. radioactive metal over a 
non-radioactive surface, the transfer of activity to this surface being measured. I have 
recently had occasion to use a method based on a different principle, that of radioactivation 
analysis. 

A copper surface was slid over a steel surface, and as a result some copper fragments 
were transferred to the steel. The steel specimen was then activated in a nuclear pile, the 
B.E.P.O. at Harwell, and within a few hours of its removal from the pile an autoradiograph 
was made of its surface, an exposure of ten seconds being given. The autoradiograph showed 
a light fogging corresponding to radiation from the iron in the steel, and a few intensely 
black spots produced by the copper fragments adhering to the steel surface (Figure 1). 
It will be seen that, if the method is to be successful, there must be a large difference in the 
specific activities of the two metals. Here the difference is by a factor of more than 4,000 
(copper : 50 mc/gm., iron: 11 wc/gm.) (Isotope Division, Harwell, 1950). 

In this particular case it also proved possible to detect the steel fragments transferred 
to the copper surface during sliding. The copper surface was also activated in the pile, 
and an autoradiograph was obtained some two weeks after its removal from the pile. By 
this time the activity due to the copper had practically disappeared (half-life of 
84Cu=12°8 hours) whilst that of the steel had been reduced but little (half-life of 
55Fe = 3 years, of **Fe=47 days). An exposure of three days showed up the steel fragments 
quite well, there being a background due to impurities in the copper (Figure 2). It is 
important in this experiment that the metal present in bulk does not contain appreciable 
amounts of impurities that are strongly activated in the pile, or else the background will 
be too high. 
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The autoradiographs, of clean surfaces of area 14 cm?, slid under a load of 2kg. ata 
speed of 0-1 cm/sec. for a distance of 0-5 cm., show that copper transfer to the steel has 
been in the form of some ten fragments ranging in size from 3X10 to 10-8 gm., as 
estimated from the amount of blackening by the method of Rabinowicz and Tabor (1951). 
The steel fragments in the copper have been far more numerous and widespread, there 
being perhaps 100 with mass in the range 3 x 10-8 to 3 x 10-® gm. The two sets of fragments 
do not correspond spatially. 

It is usually considered that the softer metal of a combination, in this case the copper, 
will wear away much more than the harder one, the steel. In these experiménts the copper 
fragments have the greater mass whilst the steel fragments are more numerous, and it is 
not possible to draw definite conclusions, except perhaps that hardness is not the only 
factor that determines the wear. 

The procedure used in this expcriment has several advantages over the standard method 
of measuring metal transfer during sliding: (a) It facilitates the use of large sliding 
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Figure 2. Autoradiograph of coppert 
surface showing a large number of}! 
steel fragments as well as background 
fogging. 


Figure 1. Autoradiograph of mild steel surface showing 
a number of copper fragments; some background 
fogging is visible in the corners. 


specimens, since strongly radioactive specimens need at no stage be handled in the friction 
laboratory. (6) It eliminates the problems which arise when, as so often happens, it is 
required to clean a specimen by removing metal to produce a fresh surface. If the specimen 
is radioactive, the procedure may be very hazardous. (c) Since very little handling of the 
surfaces after irradiation is required—they have merely to be placed on the photographic 
film and later removed—it is possible to carry out sliding experiments at leisure and later 
activate the various specimens at one time. This is of importance when the sliding 
experiments are time-consuming and the required isotope has a short half-life. (d) If, as in 
these experiments, fragments in both surfaces can be detected, it becomes possible to 
measure metal transfer in both directions during one experiment. 

I wish to thank Dr. F. P. Bowden and Dr. D. Tabor for their encouragement of this 
work, and the Department of Scientific and Industrial Research for a maintenance grant. 
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The Theory of Oriented Overgrowth of Crystals 


Smollett and Blackman (1951) have investigated the stability of a monolayer of an alkali 

~ halide, on the surface of another alkali halide, with its ions in the positions of the continuing 

lattice of the substrate, except for a uniform normal displacement. The result of the calcul- 

ation is that these configurations are unstable in a number of cases which are known to form 

oriented overgrowths. This result is interpreted as showing that the oy of Frank and 
van der Merwe cannot account for these experimental results. 

The theory put forward by Frank and van der Merwe (1949) was that a prerequisite for 
the formation of oriented crystal overgrowths is the initial formation of a monolayer coherent 
in periodicity with the underlying lattice. This does not require that the configuration of 
the monolayer should be identical with that of an interior layer of the substrate crystal: indeed, 
the contrary must be anticipated in general, even for a crystalline substance giowing on 
itself, for the configuration of a surface layer of NaCl is certainly different from that of an 
interior layer of NaCl by reason of the Verwey displacements (see below). 

Frank and van der Merwe employed a simple model of similar particles on a substrate with 
a different natural spacing to find the order of magnitude of the greatest misfit in which 
coherence of periodicity between monolayer and substrate would be retained in stable or 
metastable configurations. When the parameters of the model were adjusted to correspond 
to a Lennard-Jones inert gas model, critical misfits of 9% for stability and 14°% for metast- 
ability were found ; with larger misfits coherence is lost by formation of dislocations 
between monolayer and substrate, i.e. by a general expansion or contraction of the monolayer, 
coupled with inhomogeneous deformations. 

Smollett and Blackman investigate another problem. They consider a model represent- 
ing an infinite monolayer of one alkali halide deformed to the configuration of an interior 
layer of a second one, and placed on the surface of a semi-infinite crystal of the latter. Their 
elegant methods allow them to discover the limits of stability of their model against small 
homogeneous deformations only. Non-uniform displacements normal to the surface, and 
dislocations, are thus both excluded a priori. It is found that if the positive or negative 
misfit is conceptually increased from zero, instability first appears in shearing deformations. 
These instabilities can be given a simple physical interpretation. When the monolayer is 
unduly compressed, the repulsive forces between its ions become dominant, requiring a 
change of two-dimensional crystal structure to a closer packing. When the monolayer is 
unduly extended the repulsive forces are inadequate for its stability, and the ions of opposite 
sign must tend to approach each other to form ion pairs. If restricted to homogeneous 
deformations, the monolayer can only reveal these tendencies by commencing one or other 
of its principal shearing deformations. If it were reasonable to retain the supposition that 

positive and negative ions of the monolayer lie in the same plane, while allowing inhomo- 
geneous two-dimensional deformations, we could infer that Smollett and Blackman have 
calculated an upper limit to the degree of misfit permitting stability in the proposed 
configuration. 

Moliére, Rathje and Stranski (1948, 1949) calculate that inhomogeneous displacements 
producing ion pairs occur even in the surface layers of some of the alkali halides by themselves. 
While the present writer believes that this particular result arises from defects of the model 
used, it is clear from their theory that the same ought to happen when an alkali halide 
monolayer is deposited on the surface of an alkali halide of larger lattice spacing. On the 
other hand, a plausible equilibrium configuration for the excessively compressed monolayer 
is one in which alternate rows of ions parallel to a cube edge are shifted in opposite directions 
along their length. However, before reaching any definite conclusion in either case one 
should pay attention also to inhomogeneous displacements normal to the crystal surface. 
According to Verwey (1946) these are important in the surface layer of any alkali halide; the 
positive and negative ions lie at different levels. The principal cause is ionic polarization, 
which becomes important in a surface layer. This is disregarded by Smollett and Blackman. 
It is considered by Moliére et al., who, however, disregard the difference between positive 
and negative ions, while neither considers inhomogeneous displacements normal to the 
surface. 
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Some day all these considerations may be taken into account in a way which has not yet — 
been attempted : i.e. a combination of the approaches of Verwey and Moliére et al., possibly © 


made more brief and elegant by borrowing from the mathematical methods of Smollett and 
Blackman. We shall then know, provided the crystal model employed is an adequate 


approximation, what is the equilibrium configuration of a surface monolayer of specified 


composition when constrained to have the periodicity of the underlying lattice. "Then we 


shall be left with the question: does it stably retain this coherence? The only model yet 


available with which to attempt an answer to this question, crude as it may be, is that of 


Frank and van der Merwe. ‘The writer supposes that it gives the correct order of magnitude 


for critical misfit, irrespective of the complications. 

Smollett and Blackman’s Table 2 requires no comment except with reference to the 
entries (2a)—cases in which parallel overgrowth occurs in spite of drastic differences in 
lattice parameter. Pending a searching experimental examination, the most natural 
supposition would be that these have a similar explanation to that given by Menzer (1938) for 
the case of silver on sodium chloride—the intervention of a layer of submicroscopic twins. 


The writer must express his thanks to Drs. Smollett and Blackman for allowing him to — 


see a copy of their manuscript in advance of publication. 


H. H. Wills Physical Laboratory, F. C. FRANK. 
University of Bristol. 
24th July 1951. 
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In his comments on our paper (Smollett and Blackman 1951) Dr. Frank has raised two 
points, one of a theoretical nature. In the cleavage plane of an alkali halide there will be a 
different position of the positive and negative ions in a direction perpendicular to the surface. 
It is suggested that this may affect our conclusions markedly. 

This effect, attributed by Frank to Verwey (1946), has been discussed by Madelung (1916), 
who made experiments, without success, to show the existence of the double layer which is 
expected to form when rocksalt is cleaved. 'The magnitude of the dipole moment, as cal- 
culated by Verwey for rocksalt, is, in the writer’s opinion, too large to account for the experi- 
mental observation of Tull (1951), that the mean inner potential of rocksalt was the same 
when flat (100), (110) and (111) surfaces were used. Nor is it clear how this type of distor- 
tion of the surface could stabilize a highly unstable configuration. 

This point is conceded, in part, by Frank in that he advances another explanation for the 
existence of apparently large misfits (20-50% : Table 2(b) of our paper); the contact plane 
of the deposit and substrate is supposed to be some plane other than the cube plane, and 
subsequent twinning leads to cube planes appearing parallel to the substrate surface. In 
support of this purely ad hoc assumption, Frank quotes the suggestion of Menzer (1938), in 
an analogous case where silver is deposited on to a cleavage face of rocksalt. This appears 
to be a (100) orientation with a misfit of 27°4,—again a ‘ drastic’ case. 

As this suggestion by Menzer has been used before (van der Merwe 1949) and seems to 
serve the Frank—van der Merwe theory as a saving clause, it is not irrelevant to point out 
that at least two independent investigations (Raether 1951, Tull 1951), using the electron 
diffraction technique, have shown that the contact plane of the silver is indeed the cube 
plane. The misfit of 27% must therefore be taken seriously. 

Frank has emphasized, as we have also done, the need for further experiments. Some 
are already in progress in the Physics Department of this College. Attention should also be 
directed to those experimental results already available. Uyeda (1942) has evaporated a 
layer of silver, of mean thickness 5 a., on to the cleavage (001) surface of molybdenite. The 


, 
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electron diffraction pattern shows the (111) plane of the silver to be parallel to the surface of 
- the substrate (there being two equivalent orientations), and the breadth of the silver ‘spots’ 
_ shows that the size of the silver crystals is only 15 a. The ‘side-spacing’ of the silver is 
_ normal and differs appreciably from that of the substrate, the diffraction pattern from which 
is also visible. There is here no question of a ‘ pseudomorphic’ layer. The misfit involved 
is 9%. 

Another experimental investigation (Lucas 1951) on a thin coherent film of ZnO, grown 
chemically on a cleavage face of zinc, also shows the normal ZnO electron diffraction pattern 
as well as the pattern from the underlying zinc. The thickness of the layer is less than 25 A.. 
Here again there is no trace of the initial ‘pseudomorphic’ layer demanded by the Frank— 
van der Merwe theory. 

I am indebted to Dr. D. W. Pashley for drawing my attention to the work of Uyeda, and 
to Dr. Frank for sending me his comments. 


Physics Department, M. BLACKMAN. 
Imperial College, London S.W. 7. 
8th August 1951. 
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The Structure of Thin Layers of Zinc Oxide Grown on a 
Zinc Single Crystal 


Finch and Quarrell (1933, 1934) reported some electron diffraction experiments which 
they interpreted as showing that the growth’of very thin layers of zinc oxide on zinc resulted 
in the crystal lattice of the oxide being deformed by the influence of the base metal to give 
a structure differing from that of zinc oxide in bulk. They concluded that the basal (001) 
planes of the oxide grew parallel to the corresponding planes of the zinc so as to fit exactly, 
there being a corresponding expansion along the c-axis of the overgrowth so that the 
volume of the unit cell remained unchanged. 

This conception of ‘basal plane pseudomorphism’, as it was termed by its authors, 
has recently been brought into prominence. It forms an essential part of the mechanism 
of oriented chemical growth as discussed by Cabrera and Mott (1949) in terms of the theory 
put forward by Frank and van der Merwe (Frank and van der Merwe 1949, van der Merwe 
1949). However, the idea has been by no means universally accepted, and in particular the 
interpretation of the above-mentioned zinc oxide experiments (the principal experimental 
evidence for its validity) is open to doubt, as has already been pointed out (Thomson and 
Cochrane 1949). In view of these considerations it seemed highly desirable to carry out 
some further experiments in the hope of producing more conclusive evidence of the nature 
of very thin films of zinc oxide on zinc. 

The best line of attack on the problem appeared to lie in the use of single crystals of 
zinc. ‘These crystals have the useful property of cleaving parallel to (001) planes and so 
giving, in a relatively simple manner, the plane of most interest in the present investigation. 
Moreover, this one plane can then be studied without the experiment being confused by 
effects from other crystal planes. 

The first cleavage faces of zinc that were studied by electron diffraction were found to 
be covered by an oriented layer of zinc oxide. This layer showed the normal structure but, 
as in the similar experiments reported by Raether (1950), no electron diffraction pattern 


61-2 


944 Letters to the Editor 


could be obtained from the underlying zinc crystal, and so it was only possible to estimate 
a lower limit for the thickness of the oxide. However, with improvement in the techniques 
of cleaving and mounting specimens it was found that the time of exposure of a fresh 
surface to air at atmospheric pressure could be reduced to less than one minute. When 
such a specimen was examined with a 55 kv. electron beam the pattern of Figure 1 was 
obtained. 

Figure 2 shows how the pattern of Figure 1 may be split up into its components. The 
dots represent the zinc pattern, and the circles that of the zinc oxide. Since the zinc oxide 
completely covers the surface of the zinc (no zinc pattern being obtained at 30 kv.) the 
diffracted beams emerging from the metal must pass through a layer of the zinc oxide before 
reaching the photographic plate. In doing so they may be expected to be diffracted a 
second time in the oxide, giving an oxide pattern around each metal spot as primary. 
In general these secondary diffracted beams would be difficult to detect on account of their 
weakness, but in the production of Figure 1 the crystal was set so that much of the diffracted 
intensity was thrown into the beam giving the large spot on the centre line, with the result 
that the secondary pattern arising from this beam as primary may be seen. This secondary 
pattern is indicated by strokes in Figure 2. 


¢ 
Figure 1. Reflection pattern from zinc Figure 2. Analysis of pattern Figure 3. 
cleavage surface covered by a very of Figure 1. 
thin layer of zinc oxide. [110] 
azimuth. 


The fact that an electron diffraction pattern can be obtained from the zinc mother 
crystal through the overlying layer of zinc oxide makes it possible to set an upper limit to 
the thickness of this layer. Thus if d (Figure 3) is the oxide thickness, / the total path 
length of the electron beam in the oxide, and 6 the angle of incidence of the beam, 2d=61. 
In its passage through the lattices the beam is of course weakened, and with sufficient path 
length it will become so diminished that no detectable pattern will be obtained from the 
metal substrate. On the basis of transmission experiments Raether (1950) estimates that 
the maximum path length in zinc oxide for electrons with the energies used in electron 
diffraction work is 500 to 1,000 a. That is to say, with 7 greater than this value no zinc 
pattern would be obtained through the overlying zinc oxide. Therefore, by taking the value 
of 1,000 a. for / and putting 6=0-05, which was the minimum value for which a zinc pattern 
was obtained, we get dmax—=25 A., which can be put as an upper limit to the thickness of the 
oxide film grown in this experiment. 

Thus the experiment shows that an extremely thin zinc oxide layer has been grown, 
no more than ten molecular layers thick. Such conditions are most favourable for the 
detection of the pseudomorphic layer suggested by Finch and Quarrell, and yet the oxide 
clearly shows its normal structure, with no trace of the pseudomorphic spacings. The 
experiment provides no support for the views of Cabrera and Mott (1949) and of Frank 
and van der Merwe (1949), which have also been objected to on theoretical grounds 
(Smollett and Blackman 1951). A fuller account of the work which is reported in part in 
this note is in preparation. 
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: Physics Department, L. N. D. Lucas. 
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9th August 1951. 
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On the Mechanical Action in Dielectrics 


In my account of this subject (Smith-White 1949) I have taken the mechanical force 


density F acting in a dielectric to be 
B= (PV) Es ee ee  eeonics (1) 


where P is the polarization of the dielectric and E is the electric force in the dielectric. 
Cade (1951) has objected that the formula (1) does not take account of the fact that a volume 
element, at the position of E, itself provides a contribution to the field of the same order 
as E, so that in finding the force on the element this contribution must be subtracted. 
This view of the matter, which is a common one, is due to misunderstanding. The point 
is a subtle one, and it has not been elucidated before. I regard (1) and the associated formula 


G=PCre, eee ee os (2) 


for the couple density as the fundamental formulae in the development of a mathematical 
theory of dielectric behaviour, based on a proper interpretation of the Poisson—Kelvin 
hypothesis. Cade does not object to (2). 

According to the hypothesis mentioned, the behaviour of dielectrics under electrical 
influence may be accounted for by assuming that each volume element dv of the dielectric 
acquires an electric moment P dv; this I take to mean that a dielectric body, considered 
from the macroscopic point of view, behaves, in electrostatics, as though it were simply a 
continuous volume distribution of electric moment with density P. On this view we are 
faced with the theoretical problem of determining the behaviour of continuous distributions 
of polarity; and until this question is decided the Poisson—Kelvin hypothesis as stated 
above does not furnish any theory of dielectrics at all. 

In electrostatic theory we ordinarily contemplate the following kinds of electrical 
distributions : (1) discrete distributions of point electric charges, g;,... , 3 (ii) continuous 
distributions of electricity defined by a charge density p; (iii) discrete distributions of point 
electric dipoles, p,, . . ., p,; (iv) continuous distributions of electric polarity defined by a 
moment density P. Now the primitive electrostatic law prescribes the mechanical action 
on each element in the type (i). But it does not apply directly to give the mechanical 
action on the elements in any of the other kinds of distribution. It is easy to infer the 
appropriate extensions of Coulomb’s law which enable us to specify this action on the 
elements of distributions (11) and (ii). This, however, involves an extension of the primitive 
hypothesis of the same kind as that emphasized by Kellogg (1929) in connection with the 
application of Newton’s law of gravitation to continuous distributions of matter. On the 
other hand, if we try to write down the mechanical force which acts on an element of the 
distribution (iv), by an integration over the distribution of a contribution from each of the 
other elements, we find that the integral which presents itself does not converge. 
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If a theory of electrical distributions of type (iv) is to be in line with the theories of the 
other three types of distribution it must be developed from a hypothesis concerning the 
mechanical action on the elements of the distribution. As just mentioned, the obvious 
procedure for finding. such a hypothesis fails. There is, however, another course open 
to us, and this can be regarded as a uniform procedure for all four types of distribution. 
We begin by noting that the mechanical action on the elements in (i), (ii) and (11) may be 
specified in an alternative, mathematically equivalent, way by first introducing the potential 
function ¢, then defining the electric field E=—Y¢, and then specifying, for (i) that the 
force acting on a charge q is gE, for (ii) that the force density acting at a place where the 
charge density is p is pE, for (iii) that the force acting ona dipole pis (p. V)E, and the couple — 
acting on a dipole pis p XE. ~ ; 

In line with these specifications for distributions (i), (ii), and (iii) we can specify the 
mechanical action on an element of the distribution (iv). Suppose the distribution occupy _ 
the volume v. Define the potential by 


1 
s-[ P.v(5) du; myn ae Ce (3) 


the electric field by E=—YV¢. We now assert that the force density acting on an element 
of the distribution is given by (1) and the couple density acting on it is given by (2). These 
statements amount to a hypothesis which I regard as the appropriate extension of Coulomb’s 
law for point charges, enabling us to construct a mathematical theory of distributions of 
electric polarity. The value of this hypothesis for the physics of dielectrics is to be judged — 
by a comparison of its consequences with the actual behaviour of bodies. As I have shown, 
one consequence of this theory is that it enables us to express the conservation of energy 
in dielectrics in a more satisfactory way than has been done before. 

Thus the formulae (1) and (2) constitute my interpretation of the analytical significance 
of the Poisson—Kelvin hypothesis. Formula (3) is a mere mathematical auxiliary to the — 
statements (1) and (2). In the literature of electrical theory (1) and (2) are generally 
neglected, while formula (3) is usually supposed to embody the whole content of the 
Poisson—Kelvin hypothesis. Transforming (3) by Green’s theorem, if f be the boundary 


of v, we get P “~P 
= { Sues i dW F Joe ane (4) 
fa v 


Yr 


where P,, is the component of P directed along the outward normal to f. This formula is 
usually interpreted by saying that the moment distribution is equivalent to a surface charge 
distribution of density o=Py, and a volume charge distribution of density p=—div P. 
These are Poisson’s equivalent distributions. On this interpretation we would be justified 
in translating the fundamental hypothesis to the form “ dielectrics under electrical influence 
behave as though they were distributions of electrical charge”’. Few physicists would 
accept the statement in this form, and the notion of equivalent charge distributions is 
in fact a fundamental mistake. Poisson’s distributions are equivalent only as regards the 
potential function $; they are not equivalent as regards the mechanical effects. 

Finally, I remark that if we develop in detail the theory of point dipole distributions 
and compare it with the theory of continuous distributions of polarity based on formulae (1) 
and (2) we see that the two theories are mathematically similar. We may note that in 
the. case of systems of hypothetical point dipoles there is no analogue to the Poisson’s 
equivalent charge distributions, and it is not possible to make the error which has confused 
electrical theory for so long. Moreover, the subtle difficulties in the formulation of a theory 
of continuous moment distributions are of no real consequence for physics. If we refer to 
such distributions as a matter of mathematical convenience we must only be sure that our 
theory of them is properly parallel to the theory of discrete distributions. 


Department of Mathematics, W. B. Sm1itTH-WHITE. 
University of Sydney. 
24th July 1951. 
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The Long-Range Tritons from the Reaction °Be(d, t)*Be 


The emission of tritons as a result of the bombardment of beryllium with 
| 1 mev. deuterons was identified chemically by O’Neal and Goldhaber (1940). Using an 
ionization chamber, Williams, Haxby and Shepherd (1937) measured, at 225 kev. deuteron 
' energy, the range of the tritons emitted from a thick (5 cm. air equivalent) beryllium foil. 
‘The energy of the reaction was thus estimated to be 4:32+0-06 mev., which is about 
| 7% less than that calculated from the masses given by Tollestrup, Fowler and Lauritsen 

(1950). Recently Resnick and Hanna (1951) have reported a value of 4:61 + 0-04 mev. for 
' this reaction energy which is in agreement with the value given from the masses. 
In the present experiment, the photographic plate method was used in detecting the 
_ particles emitted from a beryllium foil of thickness 0:4 cm. air equivalent due to bombard- 
| ment by a deuteron beam of 7-70 Mev. mean energy. For this purpose a camera described 
| by Rotblat, Burrows and Powell (1951) and Ilford C2 photographic plates with emulsion 
_ 200 p thick were employed. 

Search was made for the long-range tritons from the reaction where the residual 
nucleus ®Be was left in the ground state. Apart from the general appearance of a triton 
track by comparison with a proton track in the emulsion, the group of tritons was identified 
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by the agreement of the Q-values calculated for different angles of emission. Measurements 
were carried out for eleven different angles in the region from 20° to 60° with respect to the 
direction of the incident deuterons and representative histograms at angle 20° and 40° are 
shown in Figure 1. For each angle of observation, the group was fitted with a gaussian 
curve from which the mean range and half width were measured. Since the tracks were 
dipping in the emulsion, the true range was derived from the known mean angle of 
approach of these particles to the examined area of the photographic plate. Also, a 
correction to the energy of the emitted triton due to its passage through the thin beryllium 
foil was made for each angle of measurement. Using the range-energy curve (El-Bedewi 
1951) the mean energy of the triton group was found. ‘The Q-values were thus calculated 
for different angles of observation and the mean value for the present reaction can be 
stated as 4-67+0-03 mev., thus confirming the work of Resnick and Hanna (1951). 

Also, the angular distribution of this triton group was studied by counting the number 
of tracks belonging to each of the plotted histograms. By knowing the area examined 
in the photographic plate, together with its mean distance from the target, the relative 
intensity of the tritons emitted at each angle of observation could be calculated. The results 
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were then transferred to the centre-of-mass system and represented as shown in Figure 2. 
This distribution is based on measurements of a total of 2,000 tracks over eleven angles 
at intervals of 24° or 5°. 

The marked rise of the cross section near the forward direction is interesting, since it 
suggests that the mechanism of this reaction is a single stage process having a similarity 
to the mechanism of deuteron stripping (Butler 1950), in which a proton or a neutron is 


removed from the impinging deuteron without the formation of a compound nucleus, or 


rather to the inverse of this process. In other words, the deuteron was presumably to be 
regarded as sucking out a neutron during its passage past the target nucleus. If the 
deuteron were first absorbed to form a compound nucleus, a more or less symmetrical 
angular distribution would be expected. 


Nuclear Physics Research Laboratory, F. A. EL-BEDEWwI. 
The University of Liverpool. 
12th July 1951. 
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The Momentum Spectrum of Cosmic-Ray Protons 


Recently the momentum spectrum of cosmic-ray protons has been measured both at 


the top of the atmosphere and at sea level. At the top of the atmosphere Winckler et al. _ 


(1950) found a value for the exponent of the differential spectrum of —1:9+0-05 in the 
range 10° to 14 10° Mev. Van Allen and Singer (1950) give the value —2:1 for a similar 
range. In his theory of the nucleon cascade Messel (1951 a, b,c) has used a value of 
2:1 and has obtained good agreement with a wide range of experimental data. This 
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The graph gives the computed momentum spectrum of protons at sea level for an incident power 


law integral spectrum of exponent —1:1. The upper curve is obtained neglecting ionization 
loss 8. The lower curve includes ionization loss. 


theory, however, predicts that an incident power law should be reproduced throughout 
the atmosphere for all energies above the cut-off if ionization loss is neglected. If 
ionization loss is taken into account the spectrum is slightly flattened. We have computed 
the sea level spectrum for the energy range 10° to 10 « 10° Mev., using an incident power law 
spectrum of exponent —2-1 for both cases. The results are given in the Figure. It will be 
seen that in this range the points are well fitted by a power law of exponent 2-0 when 
ionization loss is included. 

The sea level spectrum has recently been measured by Mylroi and Wilson (1951). From 
their data using a least squares method we have computed the power law whose exponent 


f 
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gives the best fit. For the region 0-6 x 10? to 6 x 10° Mev. we obtain the result 2°5+0-2 
and for the region 10° to 6 x 10? Mev., 2:°8+ 0-3. The errors which give the 50% width have 
been calculated on the assumption that there is no error in the momentum determination 
and that the fluctuation in the numbers is gaussian. Both these assumptions are not 
_ strictly true and the actual errors therefore will be somewhat larger. However, these 
results are sufficiently in disagreement with the prediction of the theory to make it likely 
that the initial assumption of a power law at the top of the atmosphere is incorrect. To 
produce agreement we would have to assume that for energies higher than 15 x 10 Mev. the 
spectrum falls off more quickly than a power law of exponent —2-1. This would also be 
in agreement with the conclusions of Hilberry (1941), although his method of arriving at 
this result is open to criticism. 

We wish to thank Dr. T. E. Nevin and Mr. Franklin for their assistance, and 
Dr. J. G. Wilson for helpful discussion. 
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Residual Resistance in an Extremely Impure Semiconductor 


In a previous paper (Wright 1951), it has been shown that cadmium oxide can act 
as an ‘impurity metal’ type of semiconductor, i.e. the concentration of free electrons mg 
which is due mainly to excess Cd atoms in interstitial positions in the lattice, is independent 
of the temperature. It is considered that the impurity centres are so dense that they are 
always ionized in an impurity metal. The electrons are free in an impurity conduction 
band. Pearson and Bardeen (1949) have suggested that such a state of affairs also exists 
in very impure specimens of silicon. 

In the previous paper (Wright 1951) it was shown that theoretical treatment of the case 
of a constant concentration of free electrons in an ionic lattice indicated that the resistivity 
‘should be proportional to the absolute temperature, if the specimen were sufliciently 
degenerate. This was found to agree extremely well with the experimental results given 
for CdO. 

Investigation has continued, and interesting results have been observed with specimens 
containing an extremely large excess of Cd. The specimens were produced by pressing 
pure anhydrous CdO powder in moulds and sintering in air. The temperature of sintering 
appears to be very crucial. At temperatures of 550-650° c. there occurs in the specimen 
a colour change from brown to black which is reversible, but if the temperature is increased 
to 700-800° c. the change becomes irreversible and, upon cooling, the specimen is found 
to be dull black-purple in colour. This colour is probably due to the large excess of Cd 
in the lattice, and it is with these specimens that we are concerned here. The conductivity 
is extremely high, and the variation of the resistivity with temperature of such a specimen 
is shown in the Figure. 

It will be seen that the results fit accurately an equation of the form 


p=a7 Bs =5:39<10*-+4-76 x 10° 7. 


The term proportional to the temperature in the equation is likely to be the same as that 
obtained in the theory given previously. As the temperature falls below the Debye 
temperature ©, this part of the resistance due to the normal scattering by the lattice is no 
longer proportional to the temperature, and should tend towards a constant value of 
BQ@/2 as the temperature approaches 0° kK. However, the origin of thé constant A is of 
interest. It bears a direct analogy to the residual resistance in impure metals. A possible 
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explanation of A is the presence of intergranular boundaries, since the specimen originated 
from powder. It is likely that such a mechanism would be strongly dependent upon the 
temperature. Therefore, it seems probable that A is due to scattering by the impurity 
centres which are defects in the lattice, but which are also the sources of the electrons. 
The mean free path due to scattering by such ionized impurities has been investigated by 
Conwell and Weisskopf (1950), and was found to be dependent upon the square of the 
energy of the electron above the bottom of the conduction band. In CdO, since it is an 
impurity metal, all the impurities are ionized, leaving no non-ionized impurities which 
must be neglected on Conwell and Weisskopf’s theory. The considerations are simplified 
still further, since the electron gas is extremely degenerate, so that only the electrons 
possessing energies near the Fermi energy E* are effective in conduction, and E* varies 
but slowly with temperature. The effective energy, and thus the mean free path J; for 
scattering by impurities, are therefore independent of temperature to a first approximation. 

Thus 1/l;=ngQ; where Q; is the effective area of cross section of scattering, and 
ne is the concentration of both free electrons and impurity centres. Qj; is dependent upon 
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the dielectric constant, the effective mass of the electrons, E* and mg, but is independent 
of temperature. Since the degeneracy is high, this leads us to the simple expression for 
A=he~(3n,/877)'/?Q;. Assuming a probable value for np~5%x10?° we deduce that 
Q;~4:4 x 10-™, or that the effective radius of the impurity centre is (Qj/7)1/?~3-8 x 10-8 cm. 
The average distance apart of such centres is of the order of 12 x 10-° cm. 

The transition from semiconductor to impurity metal, as observed by Pearson and 
Bardeen (1949) in silicon has yet to be satisfactorily explained. A certain amount of 
importance may therefore be attached to Qj, since it is a measure of the ‘sphere of 
influence’ of the impurity centre when in an impurity metal. It is hoped that from 
further investigation of the conductivity of samples of CdO with different amounts of 
impurity content, and also of the Hall effect or thermoelectric power which will give the 
concentration of free electrons ng, accurate values of Qj; may be obtained, which may 
enable light to be thrown upon the problem of transition. 


Physics Department, R. W. WRIGHT. 
University College, 
Ibadan, Nigeria. 
25th April 1951. 


ConweELL, E., and Wetssxopr, V. F., 1950, Phys. Rev., 77, 388. 
Parson, G. L., and BarDEEN, J., 1949, Phys. Rev., 75, 865. 
Wricht, R. W., 1951, Proc. Phys. Soc. A, 64, 350. 
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REVIEWS OF BOOKS 


Einfithrung in das Studium der Physik, by W. FINKELNBURG. Pp. 128. (Heidel- 
berg: Carl Winter, 1950.) DM 4.50; bound in cloth DM 6.30. 


This interesting booklet is meant to be a guide for prospective German physicists. 
It contains a little of everything : the aims and problems in physics, the history of physical 
thought, the methods of research, the University syllabus, the main sections in physics and 
the type of person likely to be an ideal physicist. There is a chapter on the physicist in 
various professions which is, alas, the weakest chapter of all. It contains contributions by 
several sub-authors about the physicist as academic and school teacher, industrial researcher, 
employee in the electrical and chemical industries and about women physicists. A chapter 
is included on physical literature listing textbooks, handbooks and periodicals, and giving an 
introduction into the international decimal classification. Finally there are some notes about 
the German learned societies and professional bodies. 

In Germany physicists are trained at universities and technical academies. During the 
first two years (=4 semesters) 14 to 18 hours per week are spent on lectures on compulsory 
subjects and about three hours on optional subjects apart from the time for practical work 
andseminars. At the end of that period a preliminary diploma examination is taken. During 
the following two years the number of lectures is reduced and the time for practical and 
seminar work increased. At the end of this period the diploma examination (or teaching 
examination) is taken; this requires a miniature thesis but not original research. A further 
two years spent on original research are necessary to complete a dissertation and after an 
examination in the form of our viva, the degree of Dr.Phil. or Dr.Ing. is awarded. 

One cannot help feeling that in those chapters in which the author ‘ advises ’ his future 
colleagues, the dogmatic element becomes sometimes so strongly pronounced that it makes 
rather amusing reading; for example he says, “‘ anyone who wants to stick to fixed working 
hours and wants to be left alone in the evenings or spend his spare time on a hobby, should 
abandon the idea of ever becoming a good physicist’. The author’s only support for this 
opinion is that of a well-known mathematician who once said that ‘‘ modern physics is 
actually much too difficult for physicists ’’. 

May I suggest that the author mentions in a future edition the value of personal contact 


by attending meetings at home and abroad, of a knowledge of foreign languages, and last 


but not least, the need of a revival of the travelling scholar. A. VON ENGEL. 


Selected Topics in X-Ray Crystallography, edited by J. Bouman. (Monographs 
on Theoretical and Applied Physics, III.) Pp. xv+375. (Amsterdam: 
North Holland Publishing Company, 1951.) 38s. 


The preface of this volume accurately describes it as a series of partly interrelated mono- 
graphs, bearing on several applications of x-ray diffraction. "The monographs have been 
prepared by members of the Delft X-Ray Institute, and have been edited by Professor J. 
Bouman. 

The work is divided into eight sections, of which the first, Section A, deals with general 
methods, and includes chapters by Professor Bouman on the reciprocal lattice and the use of 
reciprocal space in discussing problems in x-ray optics, and a chapter by Dr. van Reijen on 
errors due to the termination of Fourier series. Section B deals with distortion in crystals, 
including under that heading not only the effects of the cold-working of metals, but also of 
thermal agitation. Section C deals with twinning, and the rate of nucleation during the 
recrystallization of aluminium, Section D with the x-ray investigation of photographic 
emulsions, and Section E, by Professor J. A. Prins, with the amorphous states of certain 
elements. Section F deals with quantitative work with powders, with special reference to 
problems of absorption and focusing, Section G with the crystallization of natural rubber, 
and Section H with certain biological investigations. 

The method of presentation as a series of separate monographs, while sometimes 
allowing technical details to be discussed more fully than would be possible in a general 
textbook of comparable length, makes the logical development of the arguments more 
difficult, and leads inevitably to a certain amount of overlapping, although this is not alwaysa 
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disadvantage. The introductory theoretical section is interesting, and on the whole lucidly 
set out; but it is not clear why the sharp distinction is made between reciprocal space and 
Fourier space in Chapter I. By reciprocal space, the author of the chapter appears to mean 
the reciprocal lattice, whereas by Fourier space he means the reciprocal lattice with its 
points weighted by the appropriate values of the structure factors, that is to say the values of 
the Fourier transform of the unit cell, at those points of reciprocal space which correspond to 
the reciprocal lattice points. The distinction is perhaps a little arbitrary, for the reciprocal 
lattice is the Fourier transform in reciprocal space of the point lattice upon which the crystal 
is built up, and the weighted reciprocal lattice the Fourier transform of a periodic continuous 
distribution of matter based upon that lattice when the number of unit cells approaches 
infinity. 

The book gives a very interesting account of the various lines of research carried out by 
this most active school of x-ray crystallography. It is very well produced, and the printing 
and diagrams are excellent. R. W. J. 


Crystal Growth, by H. E. Bucxtey. Pp. xv+571. (New York: Wiley and 
Sons; London: Chapman and Hall, 1951.) 72s. 


This book deserves a place in the library wherever crystals are grown or studied. It has 
useful chapters on the various methods of artificial preparation of crystals and on habit 
modification by impurities—the author’s most particular interest—the latter supported by 
an appendix of detailed information on this subject. A miscellany of significant and 
sometimes little known observations about crystals are scattered through the work, or 
collected in the later chapters. The extensive bibliographies, despite important omissions, 
are a particularly valuable feature. 

A considerable amount of space is devoted to historical accounts of the earlier theories of 
crystal growth and relevant crystal properties, moderately complete up to about 1930. 
These accounts nearly always lack precision, and are frequently quite wrong : the reader will 
do well to use them as guides to the literature only. 

The author’s information on discoveries made since 1930 is quite astonishingly scrappy, 
whether it has to do with progress in the theory of the solid state, who invented the electron 
microscope, or (what one might have expected to be within his sphere) the plain facts about 
the supercooling of water. In the circumstances, the reviewer must be gratified that his own 
contribution to the theory of crystal growth receives mention at all. F. C. FRANK. 


Theory of Electrons, by LL. ROSENFELD. Pp. xv+119. (Amsterdam: North 
Holland Publishing Company, 1951.) 15s. 


For half a century successive generations of honours physics students have experienced 
the delight of reading for the first time H. A. Lorentz’s classical work The Theory of Electrons. 
Professor Rosenfeld’s’ book of a similar title worthily sustains the high standard set by its 
celebrated predecessor. 

After a short historical introduction the second chapter discusses how the electric and 
magnetic polarization densities of Maxwell’s equations for material media are related to the 
structural properties of atoms as expressed by their electric and magnetic dipole moments. 
The modification of atomic charge and current distributions under the influence either of 
external fields or of the mutual interaction between atomic systems is then discussed, and the 
ground thus cleared for a very stimulating treatment of magnetism, electrical polarization, 
and dispersion, in subsequent chapters. 

The treatment is classical throughout but care has been taken to point out where this 
approach is inadequate and the modifications required by quantum theory are mentioned. 
It aims not so much to give a detailed account of the electric, magnetic and optical properties 
of matter, but rather to examine and clarify the underlying principles on which the inter- 
pretation of these properties is based. In this aim the author has been eminently successful. 
The approach is physical rather than mathematical but when mathematical formalism is 
introduced it is distinguished by the elegance and clarity we have learned to associate with 
Professor Rosenfeid’s writings. 

The book is based on a course of lectures given at Manchester and one feels envious of the 
students privileged to attend such a course. It will certainly constitute essential reading for 
advanced students and research workers in physics for many years tocome. _E. H. S. BURHOP. 
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Heat and Thermodynamics, by J. K. Rogerrs, revised by A. R. MILLER. 
Pp. xx+593. 4th Edition. (London: Blackie and Son, 1951.) 35s. 


The book by Roberts has become a classic in its field, used by several generations of 
students since it first appeared in 1928. It contained adequate information on practically 
any branch of its subject for all but specialists, and gave very clear and useful treatment of 
theory. 

This fourth edition has been revised by Dr. Miller of Cambridge. Without sacrificing 
the characteristic features of the original work, he has brought it up to date, especially in the 
field of low-temperature physics, where, after all, the main advances have taken place. 

It gives a very clear account of the science of heat, and can be highly recommended to all 
‘students prepared to work hard enough to master it. Va kina 


The Elements of Thermodynamics, by CHARLES Fasry, translated by G. F. 
Herrenden Harker. Pp. x+229. (London: Frederick Muller, 1951.) 
15s. 


This is a very lucid introduction to the basis of thermodynamics, classically considered 
and presented, with some account of its application to physics and chemistry, as distinct from 
engineering. It does not include any account of the statistical mechanical approach, nor of 
Caratheodory’s presentation of the classical theory, but within these limits is very useful. 
‘The translation appears to be very good and certainly does not betray itself as translation, 
but, if it is faithful, then the original is much less limpid than we hope for and often find in 
French writing. 

Others may not agree with me, but for myself I see no reason why time need be spent on 
the translation of scientific work from French to English at all, especially when the result 
has to be sold at such a relatively high price. Every science student ought to be able to 
read this in French without real difficulty. J. H. A. 


Theory and Application of Infinite Series, by KoNRaD Knopp. 2nd English 
Edition. Pp. xii+563. (London: Blackie and Son, 1951.) 35s. 


This second edition is a translation of the second German edition revised in accordance 
with the fourth German edition. As the only essential difference between the second and 
fourth German editions is a new chapter on Euler’s summation formula and asymptotic 
expansions already included in the first English issue, the two English editions differ only 
in some improvements and simplification of proofs. The references to literature have been 
brought up to date. 

Professor Knopp’s work is a well-established standard textbook on infinite series with 
two special features. One is the simplicity of the proofs in spite of the highest mathematical 
rigour maintained throughout. The second feature is an attempt to develop the subject in a 
systematic way for which we should be grateful to the author since our knowledge about 
infinite series mainly consists of separate bits of information. In particular, the various 
tests of convergence and divergence are presented in a masterly way. The price paid for 
rigour and coherence is an elaborate presentation of the idea of real number and of sequences 
of such numbers in the first two chapters (over 100 pages), including a detailed account of 
the notion of limit. The discussion is kept on a high logical level with good commonplace 
description of the notions involved. Irrationals are defined through Dedekind cuts. Mathe- 
miaticians will find this first part rather indispensable and physicists may turn to those pages 
to learn the why and wherefore of many proofs. 

The second part, Chapters III-VIII, gives an excellent account of the elementary 
properties of infinite series dealing with the ratio and root tests, power series, especially those 
of elementary functions, infinite products, and the transformation of series for numerical 
evaluation. The text is full of well-chosen illustrative examples followed by carefully 
grouped and graded exercises. 

The third part, Chapters IX—XIV, is the development of the theory, first for positive 
terms, complete and beautifully presented, then for arbitrary and variable terms including 
uniform convergence and Fourier series, and for complex terms leading to the definition of 
elementary functions for complex variable and to Dirichlet, faculty and Lambert series. 
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In Chapter XIII, on divergent series, the special processes of summation known by the 
names Cesaro, Hélder, Abel and Euler are extensively dealt with but the systematization of 
such processes by infinite matrices is not included. The last chapter deals with Euler’s 
summation formula and with asymptotic expansions. 

Bibliography and general index complete this excellent treatise in Miss R. C. H. Young’s. 
pleasant translation. A word of warning: ‘ Application’ in the title does not mean 
physical or geometrical application. It really means the application of general results to: 
important particular cases like Stirling’s formula or Dirichlet series. The production of the 
book is first rate. ; PAUL DIENES. 
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ABSTRACTS FOR SECTION B 


Some Phenomena of the Upper Atmosphere, by SYDNEY CHAPMAN. (Presidential 
Address.) 


ABSTRACT. The height distribution of the atmosphere depends on the temperature and 
the mean molecular weight of the air; recent progress in our knowledge of these two 
quantities will be reviewed, the main sources of information being outlined. ‘The motions 
(winds and tides) in the upper atmosphere will similarly be considered. 

The main gross uncertainty regarding composition at high levels relates to the proportion 
of atomic nitrogen present. Besides the major constituents that chiefly determine the mean 
molecular weight, there are important rare constituents, such as ozone, sodium and hydroxyl, 
which reveal their presence by their absorption or emission spectra; and there are others 
whose presence can be inferred from theirs. 

The composition of the atmosphere changes by escape from the top (hydrogen and 
helium) and, by addition, from below (e.g. helium and carbon dioxide) and, as shown by the 
auroral spectrum, from above (hydrogen); cosmic rays also both add to the atmosphere and 
produce changes of composition by nuclear reactions. 

Near the magnetic equator there are important ionospheric phenomena not yet fully 
explored and explained. Among these one of the most interesting is the abnormal intensifi- 
cation of electric current flow, specially notable over Huancayo in Peru. 
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A General Formula for the Conductivity of a Gas containing Free Electrons, by 
L4G. H. Huxtey. 


ABSTRACT. This paper contains a general discussion of electron drift in gases according 
to the method of free paths in which a fallacy contained in some earlier treatments is avoided. 
A general formula is derived for drift velocity in an electric field which is an arbitrary 
function of the time, and from it formulae for drift velocities in the following special cases 
of importance are deduced: drift in a constant electric field, in an alternating sinusoidal 
field, in a constant electric field with an applied magnetic field, in an alternating electric field 
with an applied magnetic field. Formulae for the electrical conductivity of an electron-laden 
gas follow from those for the drift velocities. Applications of the formulae to ionospheric 
and other studies are mentioned. 


The Refractive Indices and Dielectric Constants of Air and its Principal Constituents 
at 24,000 Mc/s., by L. Essen and K. D. FRoome. 


ABSTRACT. The refractive indices of air and its principal constituents have been 
measured at a frequency of 24,000 Mc/s. with a precision comparable with that obtained in 
the optical range. The method is based on the measurement of the resonant frequency of a 
cavity resonator first when it is filled with the gas and then when it is evacuated. The source 
of oscillations used is a Pound-stabilized velocity-modulated oscillator, and its frequency is 
measured by reference to a high harmonic of a quartz standard, with a precision of 1 part in 
108. The cavity is provided with a tuning plunger, which changes the resonant frequency 
through a range of about 10 Mc/s., and is calibrated to an accuracy of within 1 kc/s. 

Most of the measurements were made by the frequency-change method and the plunger 
was used in a narrow region only for a precise setting to resonance. It is possible by a 
larger movement of the plunger to compensate for the whole frequency change due to the 
removal of the gas and thus to work at a fixed frequency. 

The following results were obtained for (n—1)10® where n is the refractive index at 
0° c., 760mm. Hg: dry CO,-free air, 288:15+0-1; nitrogen, 294:1+0:1; oxygen, 
266°4+ 0-2; argon, 277°8+0:2; carbon dioxide, 494+1; and the value for water vapour 
at 20° c., 10 mm. Hg pressure was 60:7+ 0:1. 

The dielectric constants can be calculated from the relationship we=n?, where p is the 
magnetic permeability and « the dielectric constant, the values of (u—1)10° being taken as 
0-4 for air, 1-9 for oxygen and zero for the other gases. 

Accurate formulae are given for obtaining the refractive index of moist air at different 
atmospheric conditions, and are reduced to the following simple formula which is applicable 
for normal atmospheric conditions : 


103-49 177-4 86:26 5748 
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(14, »—1)10°= Pe Pah ap oo 

where fj, po, 3 are the partial pressures of dry air, carbon dioxide and water vapour, ¢ is the 
temperature in degrees c., and T=273+t is the absolute temperature. A value of 
1-839 + 0-002 x 10-48 E.s.U. was derived for the dipole moment of the water vapour molecule. 


Alpha-Particle Assay and the Measurement of the Thorium-Uranium Ratio in 
Radioactive Ores, by D. H. PErRSON. 


ABSTRACT. ‘The measurement of the alpha-particle activity of a ‘ thick’ radioactive 
source containing the natural radioactive series is discussed critically. 

Two methods of determining the thorium—uranium ratio are considered for use in ore or 
mineral assays. The energy discrimination method is simple in procedure and can take 
account of lack of secular equilibrium between members of a series in the radioactive 
sample. The beta—alpha coincidence method, making use of the difference in the disinte- 
gration constants of radium C’ and thorium C’, is more complicated. 

Initial experimental verification of the theory has been obtained. The estimated 
statistical error in determining the thorium—uranium ratio by the energy discrimination 
method is +5% in the region of Th/U=1, with a 4% sample (U;O, equivalent concen- 
tration) in 4 minutes counting. A similar accuracy is obtainable in the coincidence method 
with a counting périod 10 times greater. 
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The Trajectories of Heavy, Solid Particles in a Two-Dimensional fet of Ideal Fluid 
Impinging Normally upon a Plate, by C. N. Davies and Mary AyLwarp, 


ABSTRACT. An expression is derived for the flow field in a two-dimensional jet of ideal 

fluid impinging normally upon a plate situated at an arbitrary distance from a parallel-sided 

orifice. Numerical data are given when the orifice is at an infinite distance and in three other 

cases. In each example the trajectories of solid particles carried by the fluid are calculated 

by a stepwise method, and the conditions determining whether a particle will strike the plate. 
or not are found. 5 


Fluctuations in the Energy-Loss of Fast Electrons in a Proportional Counter, by 
P. H. M. RoTHWELL. 


ABSTRACT. ‘The fluctuations in the energy lost by fast electrons over a fixed small 
fraction of their range have been studied by measuring the pulses the electrons produce in 
passing through a proportional counter. The energy loss distribution calculated theoreti- 
cally by Landau is approximately gaussian with a high-energy tail. The experimental 
distribution is in agreement with the theoretical on the high-energy side of the peak, but 
is wider than Landau predicts on the low-energy side; this discrepancy is greater in krypton 
than in argon. The experimental values of the most probable energy loss agree with the 
theoretical values. 


The Mechanism of the Low-Frequency Electrodeless Discharge in Chlorine and the 
Influence of Irradiation, by W. L. Harriss and A. von ENGEL. 


ABSTRACT. The current flowing through a 50 c/s. electrodeless discharge in a short 
cylindrical glass vessel with plane electrodes at the ends filled with chlorine at 5-50 mm. Hg 
was investigated with an oscillograph. It consists essentially of between one and fifty or 
more distinct pulses per half cycle, the number increasing with the voltage. At sufficiently 
large voltage these pulses cover the greater part of a half cycle including the instant of zero 
voltage. By irradiating the discharge with light from an incandescent lamp the average 
pulse height decreases, radiation below 4,800 a. being strongly active. By irradiating only 
the central portion of the vessel the pulse height does not change. Irradiating one wall 
reduces the height of the pulses in every second half cycle, namely those pulses whose 
electron avalanches start from that wall. 

Once the discharge starts, electric charges collect on the inner surface of the vessel. Each 
current pulse consists of a series of avalanches which develop between small areas of the 
inner glass walls. One such area may discharge within a half cycle in several steps. The 
time variation of the applied field, the wall-charge field and the field in the gas are discussed 
and a new mechanism is developed which appears to apply to discharges between solid 
dielectrics. 

Atroom temperatures, chlorine molecules form an adsorbed layer, several molecules thick, 
on the glass walls. Irradiation produces photo-dissociation of Cl,, the atoms being effective 
in capturing electrons in this layer where slow electrons and many-body collisions are likely. 
This reduction of the number of secondary electrons reduces the pulse height. Adsorbed 
molecular gas layers also explain observations at higher temperatures. At higher voltages 
the interval between two successive pulses becomes so small that one pulse reduces the height 
of the following one. This is probably due to metastable Cl atoms and resonance radiation 
diffusing to the wall where they cause dissociation of molecules in the adsorbed layers, the 
chlorine atoms capturing electrons as before. It also explains why at high voltages the effect 
of irradiation becomes insignificant. 
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